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STABILITY AND PERFORMANCE OF MANNED CONTROL SYSTEMS 
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The human- operator transform is usually specified as a function of the 
controlled- element characteristics and type of input power spectrum. Thus, any 
single description for the human operator really is valid only for a set of limited 
conditions. This study demonstrates the human-operator transfer function ap- 
proach may be generalized; from a variety of published experimental data, "capa- 
bility bounds" upon the transfer function parameters are formed. With such ranges 
defined, these parameters, forming a type of variable structure, can be used as a 
mapping function to display in the complex plane, the boundaries of human- adaptive 
capacity. These boundaries contain a collection of "critical points", similar to the 
-1 + jO point on the usual Nyquist diagram; when the plant Nyquist contour is plot- 
ted in conjunction with these bounded areas (at particular frequencies), interpreta- 
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tion of system stability characteristics is possible. Several plots of such bounda- 
ries are presented for important frequencies in the human range. A man-moni- 
tored system configuration also appears to be susceptible to such stability analy- 
sis. 

Experimental verification of the theory is obtained by analyzing certain 
systems which have been studied empirically in the past and for which results 
have been published in the literature. To this end, transparent overlays are pro- 
vided which facilitate stability determination. The theoretical predictions parallel 
the experimental results remarkably well; one of the sets of data relates to a sys- 
tem which is open-loop unstable. 

A discussion of scalar performance indices and some of their limitations 
is followed by a laconic survey of methods for solving the ’’best’’ performance 
problem, especially for varying-parameter plants. The notion of a ’’vector” per- 
formance index is introduced and a technique for establishing one response as 
’’better” than another is presented. This technique rests intimately upon the 
formulation of an ’’artificial” function which, when reduced to zero, identically 
satisfies the upper bounds on system requirements. These requirements are set 
forth in the ’’vector” index. The design for a ’’best- overall” compensator to meet 
specifications upon a varying-parameter plant is obtained by using a ’’direct 
search” method (especially suited to the electronic computer). This method in- 
corporates a ’’pattern strategy” to aid in finding the next set of trial values for the 
compensator; the strategy affords rapid convergence to a solution, as is demon- 
strated by an example. 

Finally, the example itself displays the versatility of the variable-para- 
meter mappings presented previously. Setting the human-operator transfer func- 
tion equal to unity shows that the mappings have use for ’’inanimate” system design. 
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‘'Instead, we shall regard all that has been presented so far in this volume, 
and all that will follow, as merely a set of conceptual and mathematical tricks, a 
bit of scientific legerdemain, which, for reasons unknown, but for which we are 
profoundly grateful, prove to be startlingly more successful in many applications 
than anyone who examines the basic assumptions would have any right to expect. " 

Richard Bellman 
Adaptive Control Processes: 

A Guided Tour 
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CHAPTER ONE 

STABILITY AND PERFORMANCE OF MANNED CONTROL SYSTEMS 

1. 1 Introduction and Historical Background 

The unique capabilities of man, while seeming to pale in the face of 
increasingly sophisticated and automatic equipment, still offer characteristics 
either wholly unattainable by machines or, at the least, constrained in their 
attainment by weight, reliability, or economic considerations. Ultimately, the 
purpose of any good system, e. g. , economic, political or control, is to ease 
the burdens of man , either directly (automatically serving his convenience, 
transporting him to some preferred location), or indirectly (aiding him in 
attaining his desires, whether the production of goods or the retention of vehicle 
control). It is not surprising, therefore, to find that a good deal of study and 
experiment Lave been devoted to various aspects of human behavior. These 
range from investigations of a psychological variety to those of a less ambitious 
nature involving the characterization of human dynamics. Even in the latter, 
less extensive domain, any attempt at describing the rich diversity of human 
dexterity is an ambitious one. To permit mathematical description of human 
response to stimuli, the scope of the problem must be reduced (to a tracking 
task with visual inputs) and the input-signal characteristics constrained to cer- 
tain classes. Using such techniques, Tustin, working on the problem of man- 
operated-gun-fire -control equipment, first attempted to describe the human 
controller in servomechanism terms. About that time, Sobczyk^^^ and 



Tustin, A. , An Investigation of the Operator’s Response in Manual Control of 
a Power -Dr iven Gun , C.S. Memo 169, Metr opolitan- Vickers Co., Ltd., 
Sheffield, England, 1944. 

( 2 ) 

Sobczyk, A., Aided Tracking, MIT Rad. Lab. Rpt. 430, Sept., 1943. 
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Phillips performed similar work on "aided” gun-laying systems . It was 
Tustin, however, who first arbitrarily separated the input-output characteristics 
of man into a linear operator term and a "remnant" term, to account for that 
part of the human-output response not attributable to the linear operation on the 
input , Fig. 1-1. It is clear al so from Fig. 1-1 that the system under considera- 
tion restricts the human controller to processing the error signal; this configura- 
tion is called "compensatory tracking", a term assigned by the engineering 
psychologists. Certainly, not all human control can be classified as tracking, 
nor does the tracking task imply human control. A sufficiently large group of 
situations is covered by the compensatory tracking mode, however, to make 
its study very worthwhile; even in those situations where, due to other complica- 
tions, this mode does not provide an accurate description of system conditions, 
an understanding of such manual-tracking capabilities can provide insight into the 
more complex situation. 

It is evident that a mechanism so adaptable that it can "shoot marbles" 

or pilot high performance aircraft does not lend itself to description by a single 

transfer function. After World War II, several now-classic studies of the 

(2)(3) 

human controller were performed. From these and the comprehensive 

work of McRuer and Krendel, there began to evolve a general form of human- 

(4) 

operator -transfer -character is tic in which some confidence could be placed. 

It became apparent that, in general, the particular transfer description was a 
function of both the controlled element and the input power spectrum. Thus, 

In aided tracking, the human controller’s output signal directly affects both 
the position and rate of the gun; thus, it is possible to simultaneously correct 
errors in gun-mount displacement and speed. 

^^^Phillips, R.S., H. M. James, andN.B. Nichols, "Theory of Servomechanisms ", 
McGraw-Hill, New York, 1947. 

( 2 ) 

Russell, L. , Characteristics of the Human as a Linear Servo Element, 

M.S. Thesis, MIT, May, 1951. 

(3) 

Elkind, J.I. , Characteristics of Simple Manual Control Systems, Ph. D. 

Thesis, MIT, April, 1956; also issued as MIT Lincoln Lab. Rep. No. 111. 

(4) 

McRuer, D. T. and E.S. Krendel, Dynamic Response of Human Operators, 

WADC Technical Rept. 56-524, Oct., 1957. 
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Fig. 1-1 : The Human Controller in the Compensatory Mode 
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combinations of useful sets of plant dynamics and input power spectra were 
presented in a table with certain "adjustment rules" which allowed the designer 
to "set" the parameters in the human- operator transfer function in accordance 
with empirical evidence. As a result of the environmental dependence of the 
human-operator transfer function and the occasional lack of point-to-point 
correspondence between the model and actual tracking records, the term 
"quasi-linear " is used to qualify the transfer function description; thus the 
model portrays a smoothed, time -averaged type of "linearity". Since the 
publication of the McRuer compendium, additional studies have demonstrated 
further the utility of the transfer function approach, ^ The assumption of 

human operator stationarity is inherent in most of the studies involving the 
transfer function approach, although non-steady behavior is, to some extent, 
measured by the remnant term; work has been conducted to determine time- 
varying descriptions of the human controller. Although many elegant 

models have been presented in the last few years, none has been tested and 
established to the extent of the model contained in the studies of McRuer and 
Krendel. This quasi-linear mathematical description, in its most general form, 
appears as 

K e'-’-''^(l+jrT ) 

(1-1) 

° (l+jxT )(l^jxT ) 



That is, during any given measurement run, one assumes that human response 
is not time dependent. 

^^^Hall, I.A.M. , Effects of Controlled Element on the Human Pilot, Princeton 
Univ. , WADC Tech. Rept. No. 57-509, Aug., 1958. 

( 2 ) 

Smith, R. H. , On the Limits of Manual Control, IEEE Trans, on Human 
Factors in Electronics, Sept. , 1963. 



(3) 

Sheridan, T. B. , Studies of Adaptive Characteristics of the Human Controller 
Tech. Rept. No. ESD-TDR 62-351, L.G. Hanscom Field, Bedford, Mass., 
Dec . , 1 962 . 

(4) 

Wierville, W. W. , G. A. Gagne, A Theory for the Optimal Deterministic 
Characterization of the Time-Varying Dynamics of the Human Operator, 
Cornell Aero. Lab. , NASA Contractor Rept. CR-170, Feb. , 1965. 
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It is noted that the description contains a gain term, a pure lag 
(frequently associated with human reaction time), a pole with time constant 
(taken to represent neuromuscular lag), and an equalization characteristic with 
time constants T and T . This equalization capability represents the human 
operator’s capacity to "tune” his response in much the same manner a control 
engineer compensates a system by introducing lead or lag networks. Considerable 
experiment has shown that descriptions more complicated than (1-1) offer little 
added accuracy in "goodness of fit" to the data; frequently, however, the form 
of (1-1) can be simplified by eliminating some of the time constants whenever 
they are not required to obtain a good-fit to actual human- operator response 
data. In part, this represents the analytical analog of the human -operator 
adaptive capacity. Equation (1-1) is obtained by analysis of power spectra 
records; the expression is derived under the assumptions of time -invariance 
and "trained" operators, i.e,, the human-controller transfer characteristic 
represents post -learning period response -capability. This assumed familiarity 
with the task introduces sufficient response -consistency to permit mathematical 
description of the data. 



1,2 Human-Operator Adaptive Capacity 

The transfer function approach continues to offer the most promise to 
engineers involved with human-c ontr oiled systems. Studies designed to increase 

(l)-(^) 

the utility of this technique are still being pursued. This thesis presents 



( 1 ) 

( 2 ) 



Muckier, F.A. and R. W, Obermayer, The Use of Man in Booster Guidance and 
Control, NASACR-81, July, 1964. 

Muckier, F.A. and R. W, Obermayer, Control System Lags and Man-Machine 
System Performance, NASA CR-83, July, 1964, 



( 3 ) 

Summers, E.G. , and K. Ziedman, A Study of Manual Control Methodology with 
Annotated Bibliography, NASACR-125, Nov., 1964. 

(4) 

Bergeron, H. P. , J.K. Kincaid, andJ.J. Adams, Measured Human Transfer 
Functions in Simulated Single-Degree-of-Freedom Nonlinear Control Systems, 
NASA TN D-2569, Jan., 1965. 

(5) 

Obermayer, R. W. and F.A, Muckier, On the Inverse Optimal Control Problem 
in Manual Control Systems, NASA CR-208, April, 1965 
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a different approach to meeting the stability and performance requirements of 

certain man-machine systems. In the many studies conducted, there appear 

many different values for the parameters ,T, T , T and T of (1-1); these 

h Li I N 

distinct values indicate the ’’quasi -linear " character of the human-operator 
transfer function; also they simultaneously imply the adaptive capacity of the 
human controller. From these empirically derived values, ’’capability bounds” 
upon the human operator's adaptive range can be formed. Certain mathematical 
manipulations and a suitable method for displaying such bounds, when used in 
conjunction with plots of the contr oiled-element dynamics, permit stability 
estimates to be made for the combined closed-loop system. The assumption 
that the human operator regards system stability as a primary criterion is 
intrinsically contained in these estimates; it has been established in experiments 
that this is a valid supposition. Experiments also have shown that the human 
controller is most satisfied with the dynamics of a system when he is able to 
participate (and maintain control) as a simple amplifier.^ This property 

implies that any "margin of human capacity” remaining can serve as a measure 
of human satisfaction with the system ("handling qualities”); these margins 
appear as a part of the "capability bounds” display. The bounds formulated and 
used in this thesis cannot be regarded as definitive parameter ranges for large 
populations since such do not now exist. Indeed, if past experiments share any 
common feature, it is the limited number of subjects used to provide data. 
Usually, the small range of variation among well-trained subjects is cited as 
justification for limited subject sample size. The cross-sectional aspect of this 
thesis, i. e. , its use of the data from many experiments, to some extent 



^ ^ ^Birmingham, H. P. and F. V. Taylor, A Human Engineering Approach to the 
Design of Man-Operated Continuous Control Systems, NRE Kept. 43333, 
April, 1954. 

( 2 ) 

McRuer, D. T. and E.S. Krendel, The Human Operator as a Servo System 
Element, J. Franklin Ins t. , Vol. 267, Nos. 5 and 6, May and June, 1959. 
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attenuates the effect of reduced sample size. The ranges established for the 
parameters of (1-1) are as follows: 





: 0.6 - 250 


(l-2a) 


T : 


o 

1 

o 


(l-2b) 




: 0-5.3 


(l-2c) 


h ^ 


: 0 - 25 


(l-2d) 




0 

1 

o 


(l-2e) 



The units for depend upon the units of the input- signal to the human 
controller, the others have units of seconds. It should be noted that the ranges 
on and t do not restrict the analysis in any way; as is seen below, these two 
parameters enter into the method in such a way that any range desired can be 
accommodated quickly. 

1.3 The ^'Difference Vector" Plot 

For the system of Fig. 1-1 the closed-loop system function is 



V 

R 



H G 
Q P . 

1+H G 
o p 



(1-3) 



The system is stable if the denominator contains no zeros in the right-half plane. 
It is noteworthy that in a monitoring configuration. Fig. 1-2, the closed-loop 
function is 



H G +G G 
jL - op c p 

R 1+H G +G G 
op c p 



(1-4) 



Equation (1-3) does not consider the effect of n(t). The stability analysis which 
follows does not depend upon this exclusion. 

In a monitoring configuration, the human operator observes system error and 

enters the loop as a "signal processor" only after he decides his contributions 

can improve system response. To include the "no entry" condition the lower 

bound on K, (l-2a) is reduced to zero, 
h 




'm 
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Fig. 1-2 : The Human Controller as a Monitor 
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Again, the system is stable only if the denominator contains no right-half plane 
zeros. The denominator may be written 



H G 

(l+G G ) ll + TT " / 
c p \ l+G G 

c p, 



If 



G - ^ 

p l+G G 



c p 



(1-5) 

( 1 - 6 ) 



then (1-5) becomes 



(l+G G )(1+H g') 

C P op 



(1-7) 



In (1-7), the term containing the variable -parameter portion of the system is 



The term 



1 +H G' 
o p 



(1-8) 



1 +G G (1-9) 

c p 

in (1-7) presents no serious problem in stability analysis since it has a fixed- 
parameter structure. It should be noted that (1-8) has the form of the denomi- 
nator of (1-3); of course, any effects from the poles of (1-6) remain to be 
considered. Expression (1-8), set equal to zero, is the characteristic equation 
for the system; for stability this equation must have no roots in the right-half 
plane . 

1 + H G' = 0 (1-10) 

o p 

Substituting (1-1) into (1-10) and manipulating, one obtains 



(l+jxT^)(l+jxTj) 






-JTT 



(1+jiTj^) 



G' - 0 
P 



( 1 - 11 ) 
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In Eq. (1-11), if the first term is called, Q, and the second, L , then (1-11) 

n 



bee omes 

(l+joiT )|1+J®T ) 

ntL = ‘ tK = 0 

” (1+jojT^) ^ P 


(1-12) 







By an extension of the Nyquist criterion, the encirclements that -L makes 

n 

of the Q region provide stability information about the system. 

Usually, in a fixed-parameter problem the system structure depends for 

its stability upon the nature of encirclements of the ’’minus one” point by its 

Nyquist contour; with a varying parameter structure, this "critical point” 

generalizes into a ’’region of critical points”, each of which is obtained from a 

specific set of values assigned to the triple (T , T , T ). Unfortunately in order 

BIN 

to display this "region of variation” for the ranges on T , T , T , the term Q of 

J— » i N 

(1-12) must be plotted for particular frequencies uu = uu. . Thus, to interpret 

"encirclements” and infer stability characteristics of (1-12), several sets of 

such "regions of variation” must be plotted for pertinent frequencies. Fortunately, 

in the human controller problem only a very few frequencies are required as a 

consequence of the rapid deterioration, at frequencies above uu?^20 radians per 

second, of both the human operator and the accuracy of the mathematical model. 

It should be noted that (1-12) is written such that Q represents the VARYING 

portion of the system and L. the FIXED portion. In plotting these to determine 

n 

stability (discussed below), the effects of varying -human- operator gain is 
readily observed on the Nyquist curve for -L ; the human reaction time delay is 
handled by selecting pertinent values for t and viewing their overall effects. Any 
number of such T values can be included quickly and easily, since an added phase 



( 1 ) 



Satche, M. , Discussion, J. Appl. Mech. , p. 419, Dec., 1949. 



il 



shift to *>L< is all that is required, 
n 

Rather than plot the sum (1~12) and observe any rotations about the 

origin, the following technique is used. The negative of L» is plotted, i.e. , 

-L , as well as the region of variation Q; Fig. l-3(a) displays a hypothetical 
n 

situation. If a vector is drawn from the -lj^(juu) locus to the Q ’’locus”, each end 
of the vector being at the same value of uu, that ’’difference vector” is actually 
the function Q + L< and the rotations it experiences give stability information for 
that function. Thus, stability is determined from the nature of the origin 
encirclements, if any, that occur in the ’’difference vector” plot, Fig. l-3(b). 
The vectors shown in Fig. l-3(a) are for specific values of the varying para- 
meters as well as for particular frequencies. Thus, the locus of (1-12) maybe 
visualized without carrying out the summation. The replacement of Q by unity 
results in an ordinary Nyquist plot with the head end of the vector fixed at 

+ 1 +j0 and the other tracing out the -L< (juu) locus. Thus, as a dividend the plot 

n 

of -L allows immediate observation of the plant-stability characteristics in 
n 

those cases when G =0; the critical point from which stability is determined in 
c 

such instances is at + 1 + jO. 

1.4 The Mapping Function and the Stability Criterion 

With the human operator’s gain- varying capability and delay time lumped 
into the fixed-parameter structure, the varying portion of the characteristic 
equation is given by 

(l+jx'T )(l+jxT ) 

(1-13) 

(1+jxTJ 

It is shown in Chapter Four that for fixed m*, T , and T , variations in T map 

In 1 1-1 

into circular arc segments, while with the other parameters fixed, variations in 
T^ or T^ map into straight lines. In this way, the bounded region of variation 
for the ranges selected (1-2, c through e). Fig. 1-4, may be plotted, A typical 
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plot appears in Fig. 1-5; this plot is a mapping through Q of the boundaries (end 
points) of Fig. 1-4. Frequently, and in the case of (1-13), the "end point" 
mapping is not "complete". The property of completeness signifies that every 
point of the mapping Q, as the space of points in Fig. 1-4 is scanned, is also 
a point of the "end point" mapping. This property is required if rigorous 
statements regarding stability criteria are to be made; i.e., completeness 
permits the establishment of necessary and sufficient conditions for stability. 

For human-controller studies, it may be shown that the mapping of Fig. 1-5 is 
"almost complete" and its application to engineering design studies is acceptable. 
Plots of Q for the frequencies o) = 0.5, 1.0, 2.0 and 3.0 are given in Figs. 1-6 
through 1-9; these plots are "universal" in the sense that they may be used with 
any system of interest. 

With the critical points governing stability defined by a region 0, 
estimates of human-operator handling capability are made in the following way: 
if there is some area of Q unenclosed by the contour, -L , then presumably the 
human operator can assume such parameter values and stabilize the system. 

The parameter values necessary to accomplish this appear to give insight into 
the amount of effort the human operator must expend to retain control. If the 
entire 0 region is enclosed, the human operator cannot control under any cir- 
cumstances (for the capability ranges assumed and plotted). If no portion of 
the n region is enclosed, there is stability no matter what the actions of the 
human operator (within the sense of the transfer function description). 



The ultimate justification for this approach is provided by predicting 
whether the human operator can control (stabilize) certain systems. The 
systems used in this thesis were analytically equivalent to those reported on in 



McRuer, D, R. , I. L. Ashkenas and C. L. Guerre, A Systems Analysis View of 
Longitudinal Flying Qualities, WADD 60-43, Wright Air Development Division, 
Jan. , 1 960 . 
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Fig. 1-5 : Typical Human Operator Variable Parameter Mapping (o)=l) 
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Fig. 1-6 : Variable Parameter Mapping 
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Fip. 1-8 : Variable Parameter Mapping 
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Fig. 1-9 : Variable Parameter Mapping 
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the literature. The published results were determined experimentally; their 
classifications were obtained from human operator evaluations (pilot opinion 
of handling qualities). The a priori predictions parallel the empirical results 
very well, including some for systems which are open-loop unstable. 

Qualitative information concerning handling qualities and engineering insight into 
how the systems may be modified to improve their stability and the human 
operator's satisfaction with their characteristics is also provided. 



1. 5 Applications of the Theory 

In a simulator study by Ian Hall at Princeton University the effects upon 

a pilot of various c ontr oiled- element longitudinal dynamics were investigated.^ ^ 

McRuer et al performed a detailed analysis on the results of the Hall report in 

( 2 ) 

an attempt at developing a "flying qualities" theory. The McRuer approach 

follows the "adjustment rules" procedure described in Section 1.1. Another 

experimental study conducted at the University of Cincinnati by Ralph Smith 

investigated the maximum capabilities of a human in controlling certain open- 

(3) 

loop-unstable systems. The results of both these studies furnish experimental 
justification for the theoretical predictions in this thesis. 

As an indication of the analysis technique, a single example from the Hall 
study is offered in this chapter, other applications appear in Chapter Five. The 
experiment was conducted in the compensatory arrangement as shown in Fig. 1-1 



Hall, I.A.M. , Effects of Controlled Element on the Human Pilot, Princeton 
University, WADC Tech. Rept. No. 57-509, Aug., 1958. 

( 2 ) 

McRuer, D. T. , I. L. Ashkenas and C. L. Guerre, A Systems Analysis View of 
Longitudinal Flying Qualities, WADD 60-43, Wright Air Development Division, 
Jan. , i 960 . 

( 3 ) 

Smith, R.H., An Experimental Determination of the Limits of Pilot Controll- 
ability for Unstable, Oscillatory Second Order Systems, Thesis for degree of 
Aeronautical Engineer, Univ. of Cincinnati, I 96 I (also appeared as a paper. 

On the Limits of Manual Control, IEEE Trans, on Human Factors in 
Elec tr onic s , Sept., 1963). 
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with the plant (aircraft) dynamics having the generic form (short-period 
appr oximati on ) : 



G = 
P 



5(1+0. 6s)uu 

^ 

2 2 

s(s +2Cuu s+uu ) 
n n 



(1-14) 



Hall used twenty combinations of Q and uu^ to test pilot response. The plot for 

-L. = -G from (1-14) for <^=0.2, uu =1.57 is s hown in Fig. 1-10; this c onfigura- 
n p n o o 

tion was labelled C21 in the Hall report and is so labelled in Fig. 1-10. 

Vellum (or transparent) overlays of figures such as Figs. 1-6 through 
1-9 can be made and used to good advantage with system plots as presented in 
Fig. 1-10. Such overlays facilitate envisionment of whether the "difference 
vector" between -L and Q undergoes some type of rotation. In the case pre- 
sented, a glance at Fig, 1-10 and the general trend of the Q mappings for the 
frequencies given (oj = 0. 5, 1 . 0, 2, 0 and 3. 0), indicate that system stability is 

determined by the relation between the -L plot and the mappings for uu = 2.0 and 

n 

uu = 3.0. This situation results since the -L plot leaves the Q region of 

n ’ 

criticality for uu«=?2; this implies that a condition of instability may pertain. It is 
not yet clear how to take the difference vector between -L (j2) and Q(j2), i, e. , 
conceivably, it could be connected from -L (j2) to any point within or on Q(j2); 
for instability, it must be demonstrated that a clockwise encirclement occurs as 
uu oo irrespective of the selected orientation of this vector (all the Hall configura- 
tions were open-loop stable). 

The time invariance hypothesis implies, however, that once a settled 
condition of operation has been adopted by the human controller, any "selected 
orientation" of the difference vector defines the sub-region in space to be used for 
the remainder of the stability analysis. For example, if the human operator 
assumes a set of parameters near point of Fig. 1-8 (corresponding to a 
neighborhood of max T^, max T^, max T^), there results, from Fig. 1-10, an 
encirclement by the difference vector and instability. If the point F^ is assumed. 
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the human can "hold out" longer before enclosure, but eventually an encirclement 
again occurs. Unfortunately, the prospects of being at are not encouraging 
since, although the point corresponds to max T (lead) which is obviously needed, 
the location also implies which, if the physiological correlations are 

credible, corresponds to no neuromuscular lag. If there is any hope for system 
stabilization or controllability, the human operator must assume parameter 
values somewhere between arc and arc A'd' and in the vicinity of the 

origin of Fig. 1-8, also Fig. 1-5 (detail). Particularly, is this the case when 
the human time-delay is imposed upon the system. Thus, the characteristic of 
C21 projecting beyond the Q region is further aggravated by the introduction of 
the time lag, i.e., the point -L^(J2) located at approximately 5.0 on the real 
axis rotates, with the distance from the origin to 5+jO as radius arm, a number 
of degrees proportional to the time delay. 

Since C21 represents a set of dynamics for which empirical evidence of 
human dynamic response exists, it is possible to check the inferences above. A 
"best fit" human operator description was calculated by Hall and a plot of this 
function can be made. In Section 5. 3 it is shown that the human controller did, 
indeed, assume those characteristics previously implied by the analysis; further- 
more, it is obvious that his capacity is being strained severely just to maintain 
control of the C21 system. That is to say, few points in the variable mapping 
space Q are not enclosed by the plot of Fig. I-IO, particularly after introduction 
of time-delay effects. 

It is interesting that the "subjective evaluation" of this system by two 
experienced Navy pilots was a rating of "unacc eptable -quite unsuitable. " Good 
correlation between the analytical predictions and the pilot ratings also occur for 
other sets of dynamics (Chapter Five). As implied already, correspondences 
between boundary lines of the mapping and the physical capacity of a human 
operator seem to be possible. 
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K6 ^'Best Performance and *'Best- Overall" Compensation of Varying - 
Parameter Plants 

As in all engineering design, several performance measures have been 
formulated for man-machine systems to provide a so-called figure of merit or 
index of quality; typically, deficiencies have been found to exist with all these 
scalar indices. The difficulty in attempting to quantify the performance of a 
complex system by a single number is that satisfaction of a requirement on this 
unique scalar may still result in a system with other unsatisfactory character- 
istics, irrespective of formulation and weighting techniques. A 'Vector” 
performance index might offer significant advantages. Such an approach to the 
performance problem is pursued in this thesis. Hopefully, it may be suggestive 
of a standard technique, or algorithm, for handling far more severely complex 
situations than can be exhibited here. Also, one feels intuitively that design 
compromise, on an adequate compensator should be possible so as to provide 
"best” performance even over some range of plant- parameter variation. A 
simplified example of a method for satisfying such demands is given below. 

The goal is a "best-overall” compensator for any given system with a 
human controller. Fig. 1-11. The approach may be used for a plant with variable 
parameters; the solution to be presented, however, depends intimately upon an 
initial stability study conducted in accordance with the procedure previously out- 
lined. Thus, in Fig. 1-11, the result desired is values for K , a and 8 which 

c 

provide good performance for all possible values of Q, This portion of the 

thesis departs from stochastic -input-signal considerations and studies 

deterministic -input- signal response. Thus, performance maybe measured by 

the classical indices t , o and t (response rise-time, overshoot and settling- 

r s s 

time resulting from a unit step input); these follow exactly from some arbitrary 
definition and may be regarded as a vector quantity in three space. As a result 
of a stability analysis as outlined above, approximate values for K , a, and 8 
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Fig. 1-11 : Human Control with ’’Best" Compensation 
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conducive to stability are determined. In addition, the designer attempts to use 
the freedom in selecting a and 8 (for this imposed class of compensator) so that 
the plant's Nyquist curve passes through that part of the 0 region co.-**r es ponding 
to reduced task requirements on the human operator; i.e, , one establishes as 
large a "margin of human capacity" as possible. With the ranges of a and 8 that 
provide acceptable stability available, nominal values for the human controller 
can be read from the plot using the time -invarianc e assumption. Simultaneously, 
these values must retain the "high-frequency" end of the compensated Nyquist 
contour in a location that continues to satisfy the stability criterion, as indicated 
by the corresponding 0 region. Thus, the triple (T , T , T ) must have value 
such that it lies to the left of the plant's Nyquist curve (traversed in the direction 
of increasing fy) and still be within the range of human adaptivity. Also, 

Tf^, Tj and should be as small as possible (the human prefers to operate as a 
simple amplifier), particularly T , since this parameter is almost a direct 
indication of "comfortable" adaptation. 

The above analysis, conducted for a nominal C» is repeated for several 
values covering the range of l. The greatest common range of a and g, obtained 
from the various ranges resulting from the repetitive analysis, then provide 
satisfactory stability over the domain of varying Q. From these final ranges of 
acceptable a and p, the selection is made of those values of a and g which produce 
the best performance. There is no difficulty in determining 

given set of K , a and 0; the problem is how to select subsequent trials of a and 3 
c 

from the final range established, so that each iteration is at least as good as the 
last. 

The function 

f - (t -T )^ + (o -O )^ + (t -T (1-15) 

r r s s s s 

This statement assumes the human operator uses an equivalent performance- 
criterion. 
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is formed, where 

0 < t < T 
r — r 

0 < o < O 
s — s 

0 < t < T 

S — s 



That is, T ,0 and T are the maximum acceptable bounds on t ,o and t . The 
r s s r s s 

function is an artifice created specifically to facilitate finding some point on the 

surface of acceptable performance; at f = 0, the specifications are met identically 

at their upper bounds. A minimization process conducted on the function (which 

contains time-domain parameters) produces useful information regarding 

acceptable K , a and B (which represent frequency-domain parameters), 
c 

Appropriate' reductions in T^, and T at timely points in the search procedure 

provide resolution of corresponding boundaries of acceptable -performance in 

both (t , o , t ) and (K , a, B) space, 
r s s c 

A very straightforward search technique is used to find a compensator 
satisfying the requirements. ^ ^ This method of "direct search" with "pattern 
strategy" is especially suitable to minimization of a sum of squares, but its 
overriding virtues are its straightforward simplicity and compatibility with the 
electronic computer. "Direct search" means a step-by-step examination of trial 
solutions and a method of comparing each to determine which is "best" in some 
sense. Additionally, a method is afforded for choosing the next trial point. This 
last aspect forms the strategy part of the search procedure; a "pattern" approach 
is used and from this is derived the strategy's name. 

The simplified example used to demonstrate this search technique is 



( 1 ) 



Hooke, R. and T. A. Jeeves, Direct Search Solution of Numerical and 
Statistical Problems, J. Assoc. Comp. Mach. , April, 1961. 



•Xj 




es 



shown in Fig. 1-12; the human- operator description has been set at unity since 
its presence, while having no effect upon the search procedure, complicates the 
explanation. It is assumed that the ranges of a and g have been established which 
offer good stability; specific values within that range are now desired which meet 
the step-response requirements on rise time, overshoot and settling time. For 
economy of notation is redesignated K. The initial "base point" of the search 
^^1**^1* ^1^ is chosen arbitrarily or the designer may exercise a "hunch". Again, 
the function, f, (1-15) is formed; the point (^ 2 * ^2’ ^2^ said to be "better" than 
(K^, 8^) if 

f(K2, ^ f(K^,a^, 8^) (I-I 6 ) 

By f(K.,a. , B.) is meant the value of (1-15) found from the substitution of the t , 
111 r 

o^, and t^ associated with the system step-response when the compensator has 
gain and time constants, K. , a. , and B., respectively. If (1-1 6) holds, a "nothing 
succeeds like success" philosophy dictates a continuation of the search in the 
same "direction"; therefore, exploration is continued by perturbing the inde- 
pendent variables with increments whose proportions are based on the frequency 
of success achieved along each dimension, K, a, B, in the recent past. 

In the system of Fig. 1-12, the specifications are chosen to be 

T =2 seconds, O = 15 percent, T ~ 4 seconds (1-17) 

r ^ g 

also, the plant variations possible are given to be 

0. 01 < C < 1. 0 (1-18) 

Exploration is conducted using step sizes as follows: 

For K : take step 6 =0.2 



(U19a) 
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Fig. 1-12 : Compensation of a Varying- Parameter System 
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For a : take step 6^ = 0. 1 (l-19b) 

For g : take step 6^ = 0.05 (l-19c) 



The initial base point is chosen as 

K, = 1 (l-20a) 

1 

= 2 (l-20b) 

6^ = 0. 5 (l-20c) 

In (l-19a,b,c) the subscripts 1, 2, 3 are assigned to represent step 
increments along the K^,aand g dimensional axes, respectively. Usually, the 
problem and design experience suggest step sizes for each independent variable. 
During the termination procedure, the nominal step sizes are reduced until 
impractical as a consequence of solution indifference or equipment limitations. 

The search is started at the varying plant's "sluggish end" = 1.0), The 
initial base point for the search is called, b^, (l-20a,b,c); the horizontal bar 
signifies the vector (multi-dimensional) nature of the point, i.e. , 



b^ = (1,2, 0.5) 



(1-Zl) 



The step response of the system with a compensator which possesses such 
characteristics is 

(t ,o ,t ) = (4. 8, 0,9) (1-22) 

r s s 



From (1-15) and (1-17), 



f(1.2,0.5) = [4. 8-2]^+ [0-15]^ + [ 9 - 4 ]^ = 257. 81 



(1-23) 



This value for f must now be compared with the value at the next exploration 
point. The next observation is performed at the point where 

b^ +6j = (1.2, 2, 0.5) 



(1-24) 
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The response is found to have performance measures 

(t , o , t ) = (3. 5, 0, 6. 8) 
r s s 



(1-25) 



Therefore 

f(1.2,2, 0. 5) = [3.5-2]^+[0-15]^+[6.8-4]^ = 235.05 (1-26) 



Thus , 



f(b^ +6^) < f(b^) 



(1-27) 



and becomes the new temporary base point, from which exploratory 

moves along the a axis are conducted. 

The process is continued in this way, the best of three exploratory points 

(only two are shown above, generally, f is also evaluated at 3,lso) always 

leading to a new temporary base point. The procedure is performed along each of 

the axes (K , a, 3) and yields a new base point, b , after the completion of such a 
c 2 _ 

search cycle. With a successful ’’direction'* of search established (b^-b^), the 
next temporary base point used for exploration purposes is found by doubling the 
magnitude proceeding in the same direction. This technique forms the 

’’pattern strategy” portion of the search procedure since future points of 
exploration are determined as a function of the frequency of success achieved 
along each dimension in the perturbation steps such as (1-24). Great strides 
can be made in pattern moves after very few iterations. This example was solved, 
effectively, after six iterations; the geometric aspects of the first two iterations 
are shown in Fig. 1-13. It is found that the triple 



(K ,a,B) - (5, 1.0,0.05) (1-28) 

c 

satisfies the requirements 

T = 1. 0, 0 =: 10, T ^2 

r s s 



(1-29) 
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Fig.l'l3 : "Direct Search" with "Pattern Strategy” 
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over the range of 0. 01 < ^ < 1 . 0; this is superior performance to that originally 
specified. 

Empirical justification of the above approach (with a human operator in 

the loop) can be offered only after additional research is conducted for deter- 

ministic inputs driving controlled elements other than pure gains. It is 

recognized that the bulk of the evidence justifying the form of human-operator 

transfer function used in this study has been accumulated for systems driven by 

stochastic input signals. There is no reason to suppose significant changes in 

this form for deterministic inputs; the initial range selections for K ,t, T , T 

h El 

and were quite broad, but modifications, if needed, are easily incorporated. 

Other search techniques might also be used to identify those values of 

K , a and B yielding "best-overall” compensation. Any method employed, 
c 

however, must be suitable to finding the optimum of functions which are unknown 

to the observer; i.e. , the analytic form of the function is unknown or too 

complicated to be used. This was the case in the previous example; the direct 

functional relationship between (K , a, 8) and (t , o , t ) is not easily specified. 

c r s s 

First, the transient solution for given K^, a, must be obtained by performing an 

inverse Laplace transformation. From the particular transient response resulting 

the values of t , o , t can then be determined. No direct relation between the 
r s s 

time-domain performance -indices t , o , t and the frequency -domain compensator 

t s s 

parameters K , a, S can be written. The end result is that experimental 

. . (1) . 

optimization procedures become most attractive. Wilde describes several 
search techniques that could be used alternatively, e.g. , Method of Steepest 
Ascent, Contour Tangents, Parallel Tangents or Method of Rotating Coordinates, 

> 1 ' 

Such results could not be located by the author. 

^^^Wilde, D.J., "Optimum Seeking Methods," Prentice -Hall, Englewood, 

Cliffs, New Jersey, 1964. 
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1 , 7 Summary 

Chapter One provides an overview of the main points considered in this 
study. The goals, accomplishments and limitations of the investigation are 
sketched briefly and a concise description of several of the major techniques 
used is presented. Further detailed discussion may be found in the main body 
of the thesis . 
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CHAPTER TWO 

THE HUMAN OPERATOR IN CONTROL SYSTEMS 



2. 1 Introduction 

The possibility of describing human operator ’’behavior" in terms of 

mathematical relationships began in earnest sometime during World War II. 

Confronted with the problem of designing complex systems such as gun turrets 

capable of tracking fast and maneuverable aircraft, engineers found themselves 

needing considerably more quantitative information about the human operator's 

transfer characteristics than was previously required in weapons design. Whole 

teams of professionals in mathematics, engineering, psychology and physiology 

set to work to determine more adequately the bounds of human performance and 

the capabilities that could be expected. In addition to Tustin's work, efforts 

in the United States were being conducted along similar lines by Phillips; some 

(2) 

mention of this is made in the now classic text of which he is a co-author, 

(3) 

The interested reader may find a good historical sketch in Sheridans’s paper 

and the background interplay between the engineering and psychological philoso- 

(4) 

phies given by McRuer. The resurgence of interest in the human controller 
probably may be marked by Russell's paper^^^ and the definitive work by 

Tustin, A. , An Investigation of the Operator's Response in Manual Control 
of a Power Driven Gun; C. S. Memorandum No. 169, Metropolitan- Vickers 
Electrical Co. Ltd. , Attercliffe Common Works, Sheffield, England, Aug. , 
1944. 

James, H. M. , N. B. Nichols and R. S. Phillips, "Theory of Servomechanisms’,' 
McGraw-Hill, 1947. 

Sheridan, T. B, , The Human Operator in Control Instrumentation, "Progress 
in Control Engineering", Vol. 1, Academic Press, New York, 1962. 

McRuer, D. T. and E. S. Krendel, Dynamic Response of Human Operators, 
WADC Tech. Rpt. 56-524, Oct., 1957. 

Russell, L. , Characteristics of the Human as a Linear Servo Element, (M. S. 
Thesis, M. I. T. , May, 1951). 
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Elkind. ^ Military high performance aircraft, space capsules, deep sub- 

mergence vehicles and other systems whose dynamics are under the control or 
partial control of a human operator have served to intensify the interest in the 
human as a controlling element; additionally, the need for greater understanding 
of the man-machine interface has spurred research in this area. 

Almost all the human operator studies conducted to date, and this thesis 
as well, address themselves to the problem of a human controller in a tracking 
situation; in particular, these investigations have been done in the "compensatory” 
mode, i. e. , with the controller operating off the system error, Fig. 2-1. This 
topology lends itself most readily to analysis techniques; in the face of consider- 
able ignorance about the human mechanism and the fantastic description-defying 
complexity of that mechanism some simplifying assumptions are justifiable. 
Particularly is this the case when the simplifying assumptions still permit useful 
conclusions to be derived which succeed in predicting empirical results. 

Several characteristics peculiar to human-operator control are not treated 
in this study. Psychological aspects of human behavior (that is, stress, com- 
petitive situations, rewards, other similar task influences and their effects upon 
tracking) are not considered as such. Physiological loading influences (such as 
fatigue or alcohol) also have not been incorporated. This is not to discount the 
importance of these factors upon perceptual-motor behavior and their resultant 
effect on any mathematical model of the tracking situation; again, as many 
superfluous variables as possible must be eliminated if useful analytical models 
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of the human controller are to be formed. The completely successful explanation 
of human response in this situation of reduced complexity (a milestone still not 
yet attained) places the future investigator closer to the possibility of including 
these additional parameters in a more accurate description of geneial human 
behavior. 

There have been recent studies attempting a time-variant description of 

humans^ and some experimental work relating to the adaptation process. 

Almost no work has been published on man as a monitor although the problem is 

apparently under active study. Much of the published material on the monitoring 

(5) 

situation has been essaylike , although there has been some hardware generated 
by the Federal Aviation Agency as part of the supersonic transport program. This 
system, known as "Force Wheel Steering", allows the pilot to override the auto- 
matic elements; flight tests have shown considerable promise. A recent paper 
(7) 

on monitoring reports a tentative theory to explain some experimental results 
but the theory is quite preliminary and, although some model fitting is attempted, 
nothing to compare with the demonstrated validity of the quasi-linear transfer 
function is produced. 



Operating as a monitor, the human controller observes the character of the error 
signal, but instead of remaining in the loop and processing the error signal con- 
tinuously, he enters only as frequently as he believes he is required to augment 
normally "automatic" control. 
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The following definition of quasi-linearization appears in Elkind. 
"Quasi-linearization is the process of representing the characteristics of a non- 
linear device by linear transfer functions whose parameters depend upon the 
environment of the device. These transfer functions are called quasi-linear 
transfer functions because they behave linearly when the set environmental para- 
meters remain fixed and nonlinearly when these parameters are allowed to change. " 
"Quasi-linearization is particularly suited to the kinds of nonlinearities 
observed in the human operator. His most significant departures from linearity 
occur when the input statistics or system dynamics change. Since its parameters 
depend upon the input and the system, a quasi-linear transfer function can handle 
this kind of nonlinearity. " 

It is interesting to quote the following comment from a recent NASA report. 

"The potentials of quasi-linear models have not been fully investigated and further 

( 2 ) 

studies along this line are warranted. " 

In this paper a technique is developed to display graphically a range of 
human-operator capability which has been derived from a variety of such quasi- 
linearization studies; these plots are then used, in conjunction with standard 
Nyquist contours which represent the controlled elements, to predict human 
control capacity. It is shown that these predictions parallel previously published 
experimental results. Following the stability studies, there is a chapter devoted 
to performance considerations, Chapter Six. The classical measures of step 
input response characteristics, rise time, overshoot and settling time are used 
to demonstrate that a "vector" performance index can be formed and satisfied. 

It is shown that the required compensator design for the satisfaction of such 
specifications can be determined using computer search methods in combination 
with an "artificial" function which, when minimized to equal zero, identically 
meets the stipulated response requirements. 
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Z, Z Controllability and Adaptation 

The utility of the quasi-linear transfer function approach to engineering 
system design has been established many times, over recent years; even in the 
face of its many acknowledged limitations it represents a technique which provides 
useful results. This thesis uses as its starting point the quasi-linear ’’best-fit” 
functions and the empirical results that have been published in the literature, A 
fresh look at the analytical tools extant and some novel manipulation of system 
transfer functions afford a new method for incorporating the entire adaptive range 
of the human operator and his pure time delay in a stability analysis. To date 
servo analysis techniques used in man-machine studies have utilized a nominal 
transfer function where the specific values for the time constants are selected to 
reflect the gross form of the process, i, e, , human-operator transfer functions 
for particular dynamic elements (plants or processes) have been experimentally 
determined and presented in a ’’table” which is then used by the systems designer. 
The table may be regarded as offering ’’adjustment rules” indicating how to ’’set” 
the parameters in the human-operator transfer function for varied control tasks 
and input power spectra. It should be noted that these studies normally employ 
simple transfer function approximations, usually of a bilinear form, for the pure 
time lag introduced by the human operator. 

Some of the ’’inappropriate” and lightly justified aspects of the quasi-linear 
modeling techniques are subdued in this ’’new” theoretical approach to human- 
operator studies since a broad range of human-operator adaptation and control 
attributes are inherently included in the analysis technique. This broad range 
tends to reduce the sensitivity of analyses to inaccuracies in the nominal values 
for the human-operator parameters, which previously had to be specifically 
designated. Since no new experimental work is presented there is no attempt to 
put forth a new model or model matching technique; rather, the investigation 
should be regarded as providing another means of viewing the human operator in 
engineering terms with the physical justification of the theory obtained from pre- 
viously published empirical results. 

The mathematical model required to adequately describe the human operator 
for different tasks necessarily varies in correspondence with the human’s actual 
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adaptive capabilities. To date, this variation has been characterized by selecting 
specific values for the quasi-linear description as a function of the particular 
task. The method offered in this thesis is sufficiently encompassing of the broad 
range a mathematical model of the human operator must take, that it compensates 
for some of the inherent rigidity this selection procedure introduces. The method 
accomplishes this compensation by its complete display of human- operator capacity. 
Simultaneously, the method provides both a means of analysis and a channel for 
developing better intuitive understanding of human-parameter effects upon a con- 
trol system and system effects upon the human. 

For man operated systems the most important questions to be posed are: 
"Can the human ’handle' the system?" and "What can be done to ease his control 
burdens?" A comprehensive analysis technique for certain man-machine systems 
is constructed by capitalizing on studies already conducted and results previously 
published in the literature. It is, perhaps, a first attempt at a theory for coming 
to grips not merely with the human controlled system, since this, indeed, has 
been tried before, but with that system under the added complication of the 
human- ope rator- adaptive range. Thus, to answer the question, "Can a human 
operator stabilize this system?" would be useful. The method to be offered will 
allow estimation of whether a human operator can stabilize any given system 
assuming he can adapt his "parameters" over certain ranges; in general, the 
investigator need not select specific values for these parameters. Insight into 
the relative human effort required to stabilize is also provided. No approximation 
to the pure delay introduced by all human operators need be made. The ultimate 
justification for the new approach is then provided by predicting whether the 
human operator can stabilize certain systems. The systems used were analytically 
equivalent to those reported on in the literature, along with their associated human 
operator evaluations; the published classifications were determined experimentally. 
The predictions parallel the empirical results remarkably well even for systems 
with open-loop instability. Qualitative information concerning handling qualities 
and engineering insight into how the systems may be modified to improve their 
stability and the human operator's satisfaction with their characteristics is also 
provided. The theory offers a convenient and illuminating display of the effects 



41 



upon the system of the various parameters which form the human- operator trans- 
fer function. Thus, it becomes easier to observe which ones (and their associated 
physiological characteristics, Section 2.5) are causing difficulty or limiting con- 
trollability. 

2. 3 Thesis Organization 

The dissertation is effectively in two parts. The first (Chapters Two-Five) 
deals with system stability and controllability by the human operator and has been 
described qualitatively in the above paragraphs. The second part is concerned 
with system performance and introduces a type of "vector'* performance measure 
which could have considerable usefulness in future human operator studies. This 
portion of the thesis departs from the stochastic input signal and studies the rami- 
fications of deterministic forcing functions. Much experimental work in the psy- 
chological literature pertains to step input tracking and the techniques described 
in Chapter Six may have direct application there. Too often, however, the psy- 
chologist’ s reports are qualitative and include varied methodologies thus making 
cross-experimental comparisons with the engineer’s work almost impossible. 

With all its difficulties the engineer's approach at least provides a consistent 
language. Chapter Six is offered not as a solution to the performance problem but 
as a stimulant to further study of performance measures. In conjunction with this 
"vector" type performance criterion, a "direct search" method is used to obtain 
a specific compensation of a given class which provides "good overall" performanc 
in the presence of a plant (process) having variable dynamics. The results are 
applicable to automatic control systems, in general, and systems with a human 
operator, in particular. 

Additionally, in this last chapter, the displays of human operator capacity, 
offered in the first part of the thesis are shown to have utility in the search for a 
solution to the performance problem; indeed, it is demonstrated in an example 
that these displays can be used as a tool in the design of compensators for 
"inanimate" systems. 

2, 4 The Quasi-Linear Model 



In general, human-operator response characteristics are strong functions 
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of the input signal form. For signals such as single sine waves the operator tends 
to "synchronize” his response with the input; such situations where the human may 
exercise his capability for prediction have been called "precognitive". Precogni- 
tive response can be considered a form of open-loop control rather than tracking, 
although certain tracking tasks do permit this type of behavior. If both the opera- 
tors corrective responses and the input are displayed separately, the form of 
tracking is called "pursuit"; the classic example of such tracking is aiming a 
rifle at a moving target through open sights. In "compensatory" tracking the 
operator cannot distinguish between the effects of his response actions and changes 
in the input, except under zero input conditions, i. e. , the operator views an error 
display only. The analogous example for "compensatory" tracking is aiming a 
rifle at a moving target through a telescopic sight with a narrow field of vision. 

As mentioned above, this study restricts itself to the compensatory case. 

A very brief derivation of how the quasi-linear human- operator transfer 
function is obtained, follows. From Fig. 2-1, 
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where N( j uu ) is the transform of n(t), a term arbitrarily added at the operator's 
output to account for that portion of the output not described by the linear transfer 
characteristic; this term is called "the remnant. " 

Also, 
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where the asterisk signifies the conjugate and the subscript T indicates the 
Fourier Transform has been taken on the associated truncated time function, e, g. , 
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Fig. Z-1: The Human Controller in the Compensatory Mode 
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under the assumption that the signal and noise are uncorrelated. Therefore, 
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i. e. , the cros s- spectrum between input and operator output divided by that 
between input and error yields a set of data in complex form which represents the 
human operator block of Fig. 2-1 (solid lines). 

These data can be fitted closely by a linear analytic relation of the form: 

K. 

■ (l+j«JTj^)(l+juuTj) 

Jackson has shown that negligible increased accuracy results if (2-1) is 

made any more complicated. ^ Equation (2-1) provides a “best-linear-fit" in a 

minimum mean- square error sense assuming random inputs of a Gaussian type, 

( 3 ) 

i. e. , normal amplitude distribution. The human operator adjusts his para- 
meters to some unknown internally-generated criterion; comparisions with the 

mean- square-error criterion seem to support its validity as a substitute analytic 
(4) 

index. The Gaussian assumption has proved to be reasonable and realistic in 
the light of the low-pass nature of most control systems, i. e. , low-pass char- 
acteristics in feedback systems tend to make the amplitude distributions Gaussian 
(central limit theorem)!^^^^^ Additionally, the Gaussian signal lends itself to 



This is not so restrictive as might appear since any criterion which is a 
monotonically increasing, even error function is also minimized by the resulting 
( 2 ) 

linear fit. 
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experimental work; also, Booton has shown that "the best mean- square 

error approximation to any nonlinear devices having a Gaussian input will 

( 1 ) 

be realizable and is determined by (2-la)*’- 
F rom 



§ 

uu 




is found 



where 



uu 



H 

o 



1 + H G 
o p 



$ + 
rr 









1 



1 + H 



G 

P 



2 



nn 



which may be considered the "closed-loop remnant" spectral density. Let 
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Greatest confidence in the description (2-1) occurs when it represents 
the human operator "adequately", i. e. , when n(t), the remnant term, is 
small. A figure of merit termed "linear correlation" has been generated 
which relates the amount of linearly correlated operator output to total 
operator output. From (2-2) is found the fraction of the linearly correlated 
operator output power spectrum to total operator output power spectrum 
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If the square of the linear correlation is written 



then 



h 1 ^ t 



p = 1 - 
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Of course, the closer p is to unity the greater the confidence in (Z-1) as a 
description of human- operator response characteristics. Frequently, the 
linear correlation is good (p ^ 0, 8), thereby justifying use of the quasi- 
linear model for low-frequency conditions. Further discussion of remnant 
characteristics and possible sources for the remnant term can be found in 
McRuer. 



Z, 5 Adaptive Ranges of the Human Operator 

From various published data in the literature^^^^^^ the values assumed 

by ^ L,* assembled and from the largest and smallest 

of these values ranges for each parameter were established. These data 

were associated with controlled elements (system dynamics) having forms of 

pure gain, lags, integrators and second-order systems. Even simululated 

approximations to human operator control of the short period longitudinal 

dynamics of an aircraft and the dynamics of an F-80A aircraft engaged in a 

tail chase have been studied. In general, the most complicated human-operator 

transfer functions that could be justified by these data have the form (2-1). 

(4) 

Some double leads and lags were used in the linear fits by Hall but these 
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were isolated cases and they are of no great concern here. Again, the 
human operator adjusts his parameters to conform to each control task; the 
tasks themselves are changed by modifying the controlled element or the in- 
put power spectrum. 

Since a considerable amount of work has been expended in obtaining 
these data and they are highly indicative of human-operator capabilities they 
should be sufficient to permit upper and lower bounds to be drawn upon the 
human adaptation process. It is not maintained that the ranges given below 
definitively bound human-operator capabilities, but they are suggestive of 
limits obtained for a variety of forcing functions and controlled elements. 

The question as to whether a human operator can freely adjust these para- 
meters independently to optimize is not answered although some studies seem 
to imply he does. ^ The "adaptive behavior" of the human operator refers to 
responses to tasks AFTER the learning or training period has been accom- 
plished. Inherent in the final adaptive- range potential is the assumption that 
the operator has gained some familiarization with the system. Studies are 
currently being conducted on the learning or adaptive process; some results 
have already been reported in the literature. ^ 

Thus, the human-operator transfer function with the range of variation 
indicated below should provide a very broad, useful set of results until fur- 
ther definition of bounds is empirically determined. The validity of the ap- 
proach is not conditioned upon the absolute accuracy of these ranges which 
always can be modified when desired. 
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The '*best fit” (in the minimum mean-square sense) human- operator 
transfer function based upon Gaussian input spectral analysis takes the gen- 
eral form (2-1) , 

(1 + + jojTj) 

The ranges on the various parameters are 
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where the units for depend upon the units of the input to the human-opera- 

tor block and the units for T, T^, T^, and are in seconds. 

.To date, little work has been done on the monitoring problem, and few 

of the techniques of control system analysis have been applied to either multi- 

( 2 ) 

axis or monitoring tasks involving the human operator. The stability 
analyses presented below are also applicable to the monitoring situation. The 
analysis is conducted to establish stability conditions only for the situation 
where the human enters the system and remains as a controller. Further 



Some work has been done on sampled data models of the human operator 
based upon the philosophy that human performance is of an intermittent 
nature rather than continuous. This concept is not used in this study; 

details are available in the liter ature! ^ ^ Bekey’s work also contains a 
delightful critique on the psychological literature. 
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work on periodic entrance and exit could be useful in application to time- 
sharing a human operator, i. e. , the human is responsible for N systems 
which he may ’’trim up" on a "multiplexing" basis. As an amplification of 
the sense of the word "monitoring" as used in this study, the following quote 
is offered. "Force wheel steering operationally maintains the human 
pilot as an active control element with an automatic system for the purpose 
of having the human pilot monitor, supervise, and act as a redundant element 
for those contingincies which require pilot action. . . THE AUTOPILOT ACTS 
ON THE SAME GUIDANCE INFORMATION THAT IS PRESENTED TO THE 
PILOT; therefore, the pilot can easily monitor autopilot performance. For 
those approaches when the autopilot is operating satisfactorily, the pilot 
makes a minimum number of inputs with the Force Wheel. But the pilot can, 
upon assessment of the total situation, modify the autopilot performance using 
the same information as the autopilot by making normal control actions. . . 
pilot inputs can be made WITHOUT disconnecting the autopilot. " A possible 
block diagram display of such a monitoring system is shown in Fig. 2-2. In 
order to incorporate the "hands off" mode in monitoring, i. e. , when the 
human controller effectively eliminates himself from the loop as during the 
"off duty" period in true monitoring, the lower bound on has been reduced 

to zero. To facilitate certain manipulations, and ease the algebra in some of 
the mapping theorems, the lower bound on T^ also is considered as zero; 
this in no way restricts the generality of the results. The broadened range 
does ensure that when the stability conditions are met, greater confidence in 
the physical realization of stability follows. On the other hand, the number of 
cases which are stable is reduced; the modification required on T^ is not too 
significant however, and is justified by the simplification provided. 

There has been an attempt to associate some of the parameters in the 
human- operator transfer function with physiological characteristics; the con- 
sensus is generally as follows: 
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Fig. Z-2: The Human Controller as a Monitor 
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human- controlle r gain; normally adjusted for overall 
system stability and low-frequency performance 
pure delay (characterizes human reaction time) 
equalization; these provide the human operator's 
adjustable compensation similar to lead or lag networks 
in "inanimate” systems 

neuromuscular lag ; appears to be partially adjustable 
for the task 



This study does not enter the debate on the validity of these assignments. 



2, 6 Summary 

Normally, the human-operator transfer function is specified as a 
function of the controlled-element form and type of input power spectrum. 
Thus, any single description for the human operator really is valid only for 
a set of limited conditions. This study demonstrates the human- operator 
transfer function approach may be generalized by using a variety of published 
experimental data to form "capability bounds" upon the transfer function para- 
meters. With such ranges defined, the human-operator capacity for adapta- 
tion can be displayed and used in conjunction with controlled-element system 
dynamics to determine certain "man-machine" closed-loop stability char- 
acteristics. This technique, to some extent, depends upon the assumption 
that the human operator uses as a basic criterion the maintenance of system 
stability; this assumption has been verified as a reasonable one after many 
experiments. A man-monitored system configuration also appears to be sus- 
ceptible to such stability analysis. The difficulty experienced by a human 
operator in attempting to maintain control of the system is implied by the 
"margin of human capacity" remaining; this capacity appears graphically as 
a part of the analysis technique. 
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CHAPTER THREE 

HUMAN OPERATOR ADAPTIVE RANGES AND STABILITY THEORY 
3. 1 Introduction 

In this chapter the system transfer function mathematics are presented 
and the relations are manipulated so the adaptive properties of the human operator 
are isolated. The concept of a nominal transfer function is introduced to facilitate 
the separation of the closed-loop- system function denominator into a fixed and a 
varying part. Finally, a brief discussion of stability interpretation for this 
modified Nyquist display is given. 

Many of the concepts in this study may be found in different form through- 
out the engineering literature. Several such published papers led to considering 
the possibility of a study on the effects of variations in the parameters of a trans- 
fer function the FORM of which was known. Also, techniques similar 

in certain respects to those presented here are offered in Horowitz's text. Of 

considerable aid in contemplating the theoretical ramifications of a pure time 

(?\ 

delay was Choksy's work . The utility of the method described lies in its inclu- 
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sion of a range on the parameters of interest so a specific expression of a trans- 
fer function need not be selected. The case of parameter variation in a transfer 
function of given form is precisely the state of affairs in certain man-machine 
systems where the human operator provides an adaptive transmission link. 



3. Z System Transfer Function Manipulation 

The closed-loop system function for Fig. 2-2 is found to be 



H G +G G 
V _ op c p 

R"1+HG+GG 

op c p 



(3-1) 



It is noteworthy that in solving the stability problem for this system several other 
configurations are handled simultaneously. For example, removing G^ in Fig. 
2-2 and placing it in feedback around G^ leads to the same characteristic equa- 
tion as in (3-1). The stability of (3-1) depends upon whether any of the zeros of 
its denominator lie in the right-half plane. For stability, there must be no zeros 
of 

1 + L=1 + HG+GG (3-2) 

o p c p f 

in the right-half plane . 

Equation (3-2) can be written 

H G 

l+L = (1 + G^G )(l+ ^ ) 

^ c p 

or 



l + L = (1 + L )(1 + L, ) 

P n 

where L is the open-loop transfer function corresponding to the compensator- 
plant and is the open-loop transfer function corresponding to the human con 

troller-plant path where the plant is somewhat more complicated than and 






That is, no roots ofl+HG+GG =0in the right-half plane. 

o p c p 
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is represented by 



1 + G G 
c p 



Of course, 1 presents no serious problem in stability analysis; the tech- 

niques for stability checks and performance studies in this situation are well 
known since the structure is FIXED PARAMETER. 

The term containing the VARYING PARAMETERS is given by, 1 



where 



H G 

L = Q P 

h 1 + G G 
c p 



= H G' 
o p 



that is, 



1 + G G 
c p 



The symbol, H , represents that part of the expression which has 

variable parameters and G^ that part with parameters nominally fixed. 

Since; 

1 + L, = 1 + H G' (3-3) 

h op 

this can be written as 



= 1 +G H 

P o 

* r 



H 



H 



where H represents some nominal value of which is obtained by picking 

n 

any set of values of the varying parameters in the given range of their variation. 
Now, the zeros of 



are the zeros of 



H 



H 



1 + G H 

P 



— + G' H 



H 1 

o 



H 



nj 



(3-4) 



H 



P o 

r- 



H 



+ L 



(3-5) 



The examples in this study are for G(,= 0 in order to simplify the presentation, 
permit an intuitive grasp of the method, and allow comparisons with experimental 
results reported in the literature, since few empirical data have been reported 
for the monitoring configuration. No loss in generality occurs; the condition of 
an Lj^ with right-half plane poles, a situation representing the worst effect of a 

non- zero G , will be discussed in an example, 
c 



i 
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where 



L = G' H 

n P o 



(3-6) 

is the nominal loop-transmission with FIXED PARAMETERS, Fig. 3-1. The 

stability problem is now one of shaping the function with non-varying parameters 

H H 

o o 

L such that in the face of the parameter variation in — the zeros of -«— + L 
n hi H n 

o o 

are all in the left-half plane; the criteria ensuring stability will be discussed 
below. 

3. 3> The Nominal Function 

Of course, H should be chosen to facilitate the analysis. Equation 
o 

n 

(2-1), representing the human operator is 

(l+sTJ 

" (l + sTj^)(l + sTj) ’ 

where the final ranges considered in this study are 



(3-7) 



0 < 




< 


250 


1 < 


T 


< 


0. 4 


0 < 




< 


5. 3 


0 < 




< 


25 


0 < 




< 


0. 7 



As an example, the nominal expression for H is taken as 

n 



H = e 
o 

n 



-. Is 



(3-8) 



where it is clear that (3-8) is (3-7) with 



K|, = 1. ’ = 0. 1, = Tj = = 0 



Thus, 
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(b) 



Fig. i’l : Separation of Fixed and Varying-Parameter Elements 

(a) Original System, (b) System after introduction of nominal function, H 

n 






I 
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H - 1 s 

o 

• (l + sTj^)(l + sTj) 
^°n (Tj^s+ l)(TjS+ 1) 

-Tif = -Tyr^iTTT— 



(3-9) 



It is noted that the human-operator gain K can be lumped with the fixed elements 

n 

of (3-5), This results in a design for a single overall open-loop gain, K^, 

from which the human-operator gain needed, K , or the plant gain needed, K , is 

n p 

determined upon specification of one or the other, i, e. , K^=K K , In a human- 

1 h p 

operator study the plant gain would be fixed at some nominal value since the in- 
vestigation i'S focused upon the human operator aspects of the analysis. In general, 
this approach is used throughout, thus from (3-4) 



H 



1-fG H 
P o 



n 



H 



1 + G H 
P o 

~ r 



+sTj_^) 



(l + sTj^)(l + sTj) * 



1 + K K G H 

h P P o 



(1 + s T^) 



(1 + sTj^)(l+ sTj) 



0 (3-10) 



The study regards only fixed values of T and covers the entire range of T 
variation by repetition over a sequence of values uniformly distributed, or 
otherwise, across 0. 1 < T < 0. 4. Since the plots are in the Nyquist plane the 
additional burden is negligible and consists only of a phase shift in the 

- i 'UT 

magnitude of e being unity. 
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where the plant gain, , has been removed from G' and H has been set 

equal to e , Equation ( 3- 1 0) , put into the form (3-5), is written 



H 



n 



H 



+ L = 
n 



(l + sT^)(l + sT,) 



r 



(1 + s T, ) 



-Ts 



Throughout the remainder of this work the expression for 



= 0 

H 

o 

n 

H 



(3-11) 



is written as 



H 



Q 



H 



(l + s Tj^)(l+ sTj) 
(l + sTj_^) 



(3-12) 



and the human-operator gain is lumped into as in (3-11). It should be noted 

that the term of (3-11) has a FIXED PARAMETER structure, Fig. 3-1; each 

analysis considers only fixed values of T across the range of interest, since the 

plots are in the complex plane the effect upon L consists only of an additional 

n 

phase shift. 



3. 4. Stability Criterion 

The standard Nyquist criterion for fixed parameter systems gives informa- 
tion about the characteristic equation of the closed-loop system. Thus, if the sys- 
tem function is 

®o _ KP(s) 

" 1 +KFP(s) 

then system stability is determined by whether the characteristic equation 

l + KFP(s) = 0 

has roots in the right-half plane. If the Nyquist plot of the open-loop transfer 
function encloses the -1 point, the number of such encirclements and their sense 
gives the number of zeros minus poles that l + KFP(s) has in the right-half plane. 

If the system is open-loop stable, i. e. , no poles of KFP(s) in the right-half 
plane, the number of encirclements gives directly the number of zeros in the right- 
half plane. Thus, in this instance, any clockwise encirclement is sufficient to 
guarantee instability. If the system is open-loop unstable there are poles in the 



I 




I 
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right-half plane, and these must be accounted for in interpreting whether or not 
zeros are in the right-half plane. Clockwise encirclements are regarded as posi- 
tive and counterclockwise encirclements as negative. For example, if 



KFP(s) 



K 

s(“ 1 + ST 



the system is open-loop unstable (there is a right-half-plane pole). However, the 
Nyquist plot yields one clockwise encirclement, i. e. , 



number of zeros - number of poles =Z-P = +1 



but P is known to be unity thereby yielding Z=2. In this case, the system is 
closed-loop unstable, the characteristic equation has two roots in the right-half 
plane. 

Can the Nyquist criterion be applied to systems containing time lags? 

Such a condition is always present with a human operator in the loop! The Nyquist 
criterion indicates the number of zeros of a function in the right half of the s 
plane by the number of revolutions the mapping of the function makes around a 
critical point of another complex plane for all values of s bounding the original 
right-half plane. The primary condition ensuring the validity of the Nyquist 
criterion is the analyticity of the function; since the exponential term representing 
the time lag effect is analytic in the region of interest the criterion is directly 
applicable without modification . 

The stability criterion used throughout this study follows from (3-11) 



Q + L 



(l + sT^)(l+sT^) 



(1 + sTj 



+ K„G' = 0 

T p 



(3-11) 



More detailed justification can be found elsewhere. 



(1)(2)(3) 



Satche, M. , Discussion, J. Appl. Mech . ,p. 419, Dec., 1949. 

Dzung, L. H. , The Stability Criterion, "Automatic and Manual Control", 
Tustin, A. , Cranfield Conference papers, Butterworths (publishers), 1952. 

Choksy, N. H. , Time Lag Systems, "Progress in Control Engineering", Vol. 
1, Academic Press, N. Y. , 1962. 



( 3 ) 



61 



The stability of the system depends upon the number and type of encirclements 
the plot of (3-11) makes of the origin of its complex plane. Variation in 
and causes 0 to take on the form of an area, however, and rotation of the 
’Vector” representing (3-11) as the values of s traverse the bounds of the right- 
half plane must be interpreted with care. All points within the region of variation 
which result in a rotation of this ’Vector” such that instability is established 
represent combinations of values for the triple T^, T^) which yield an un- 

stable system. If no points exist which provide a stable result this reflects a 
situation where the human operator cannot stabilize (or control) the given system. 
The concept will be examined in detail later by building from simple examples. 

Thus, the region (3-12) is plotted for all possible values of the triple (T , T , T ); 

H LIN 

o 

the boundaries of the regions Q = — (juj) (for various fixed uu each is called 

o 

the "hull”) are universely applicable for the selected range of variation on the 
time constants T^^, Tj, T^. This range was obtained from empirical results pub- 
lished in -the literature and was established by selecting the largest and smallest 
values that T. , T and T ever assumed, as indicated from the many different 
experiments conducted to determine them. The universal plots may be used in 
conjunction with any system of interest. The first generalized version of this 
technique seems to have been given by Jones^^^; regions of variation were not 
considered. 

Rather than plot the sum (3-11) and observe any rotations about the origin, 
the following technique is used. The negative of is plotted, i. e. , as well 

as the region of variation 0 ; Fig. 3-2(a) displays a hypothetical situation for 
explanatory purposes. If a vector is drawn from the -L^(joj) locus to the 0 



The origin here plays the role of the "critical point”. 

Jones, P. , Stability of Feedback Systems Using Dual Nyquist Diagram, Trans. 
IRE Prof. Group on Circuit Theory , Mar. , 1954. 
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**locus", each end of the vector being at the same value of uu , that ’’difference 
vector” is actually the function Q+ and the rotations it experiences give sta- 
bility information for that function. That is, stability is determined from the 
nature of the origin encirclements, if any, that occur in the ’’difference vector” 
plot, Fig. 3-2(b). The vectors shown in Fig, 3-2(a) are for specific values of 
the varying parameters as well as for particular frequencies. Thus, the locus 
of (3-11) may be visualized without carrying out the summation. The replace- 
ment of Q by unity results in an ordinary Nyquist plot with the head end of the 
vector fixed at + 1 + jO and the other tracing out the -L^(juu) locus. Thus, as a 
dividend the plot of allows immediate observation of the plant- stability 

characteristics in those cases when G^=0; the critical point from which stability 
is determined in such instances is at + 1 + jO. Examples will be offered in the 
next chapter. 

It is not necessary, of course, to introduce the concept of the nominal 

human- operator transfer function H , since (3-3) 

n 

1 +H G' = 0 (3-3) 

o p 

could have been plotted in the form of either 

H = (3-13) 

o G 

P 

or 




and the "difference vector" between the two loci interpreted similarly to the 
explanation above. In (3-13) the design engineer using the techniques of this study 
would be required to plot the inverse frequency locus and in either (3-13) or (3-14) 
the general contours to be developed which define the region of variation would be 
complicated in their calculation by a linear phase shift due to the delay time. The 
use of the nominal expression for the human operator removes the delay from the 
varying portion of the plot and places the burden for its consideration on the FIXED 
PARAMETER portion, i. e. , the user of the study always plots this locus since it 
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represents his particular system. His plot, however, need be modified by the 
addition of a linear phase shift to a single contour only instead of to a bounded 
region. 

3. 5 Summary 

Under the assumption that the human operator can be described by a 
transfer function composed of five variable parameters, the system mathematics 
are manipulated so that three of these parameters (time constants) are displayed 
together. The remaining two parameters, gain and pure delay time, are lumped 
with the rest of the system; their effects upon the closed-loop stability character- 
istics are easily determined, since their influence upon the Nyquist contour of 
the plant is readily observed. It is shown that stability analysis of certain man- 
monitored system configurations also is possible. Stability estimates are rea- 
sonably accurate as long as the human-operator description used continues to 
represent the human controller adequately; this condition seems to hold for 
angular frequencies below 18 to 20 radians per second. 
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CHAPTER FOUR 

THE STABILITY CRITERION AND THE VARIABLE PARAMETER MAPPING 
4. 1 Introduction 

One intent of this study is to provide a technique for viewing human con- 
troller-system response characteristics which is more intuitively satisfying than 
previous approaches. As a result of some interesting system transfer function 
manipulations, certain man-machine systems can be studied for stability on the 
familiar Nyquist plane. The time delay inherent in the human operator is easily 
handled and the effects of many of the human operator's physiological character- 
istics are plainly observed. The goal is to introduce a new approach to handling 
the quasi-linear human-operator transfer function without excluding the enriching 
properties of human adaptive range and pure time delay. 

In this chapter a complete analysis of a very simple single-varying-para- 
meter system is conducted to provide familiarity with how the theory discussed in 
Chapter Three actually is used. Following this presentation, the multidimensional 
variable parameter mapping is discussed in detail. The mapping function is that 
portion of the human-operator transfer function which previously was isolated as a 
separate term (the three variable time constants). Certain properties of the map- 
ping which relate to necessary and sufficient conditions on system stability are 
considered. Finally, plots of the mapping are presented (also on transparent over- 
lays); it is these plots which are used in Chapter Five to predict human control 
capabilities. 

4. Z Analysis of a Simple System 

To facilitate an understanding of the modified Nyquist diagram used in this 
study consider the system depicted in Fig. 4-1. The root loci shown in Figs. 4-2 
and 4-3 are representative of this system, i. e., for the same open-loop gain K = K^ 
it is possible to have values for the variable time constant x such that a stable 
closed loop occurs, Fig. 4-2; or an unstable closed loop prevails, Fig. 4-3. The 
encircled crosses depict closed-loop poles while the others represent open-loop 
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Fig. 4-1 



: A System With One Variable Parameter 
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poles. Such a system, with only one varying parameter, can be studied by way of 
example to gain an understanding of the dual Nyquist plotting technique. Generali- 
zation to a three varying parameter mapping follows easily. 

A specific example of Fig. 4-1 where b=l and ^ follows. The closed- 
loop function has a denominator of the form 



K 



D(s) 



1 + 



= 1 +. 



2K /x 
P/ 



(sx+ l )(-2 + l)(s + 1) 



(s +-)(s + l)(s+ 2) 



If X = then 



8K 



K 



(s+ 4)(s+ l)(s+ 2) (s + 4)(s + l)(s+ 2) 



K = 8K 



Therefore, 

D = 1 + F= s^+ 7s^+ 14s + 8+ K = 0 



is the characteristic equation. From the Routh array, the stability limit for the 
root locus gain is found 



3 

s 

2 

s 

1 

s 

0 

s 



1 14 

7 (8 + K) 

-(8 + K)4 98 
7 



to be K = 90, K — ^ = 

V o 

the ju) axis is 



11. 25; the point at which the closed-loop poles cross 
7s^ = -98 



s = ±j3. 74 

. 1.T - 1 

This situation is shown in Fig. 4-4 by the solid line loci. Now suppose, x ^ 



I 




IjH 

n 

^'1 
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instead of — then 
4 

6Kv _ ^ 

(s+3)(s + l)(s+2) “ (s + 3)(s+ l)(s+2) 

and 

D=1 +F = s^+6s^+ lls+6 + K = 0 

is the characteristic equation. Again, using the Routh array, the stability limit 
is found 



1 11 

6 (6+K) 

-(6+ K)+66 

S 

to be K = 60, = 10 ; the crossing of the j uu axis by the closed-loop 

poles occurs at 

6s^ = -66 

s =-11 
/. s = ±j3. 32 



3 

s 

2 

s 

1 

s 



0 

s 



The loci for this system configuration are shown dotted in Fig. 4-4 as is the pole 




X 



If X is allowed to vary over the arbitrarily selected range 0 < x < 5 the 
system will be characterized by both stable and unstable response properties 
depending upon the value x assumes, provided the gain is "properly” fixed. Thus, 

for K = 90 the system is at its stability limit for — = 4 , and is surely unstable 

1 ^ 
when — = 3 for the same gain, since the gain driving the closed-loop poles onto 

^ 1 

the imaginary axis for — = 3 is only K=60. Thus, at K = 90, there is system 
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stability for i > 4 and instability for — < 4, see Fig. 4-5. 

The plot of these situations using the dual Nyquist approach is now con- 
sidered. Again, letting 



D(s) 



1 + 



11. 25 

(sx+ l)(s + l)(-“ + 1) 



and with the previous knowledge that 

K = 1 1 . Z5 and x = “ 
V 4 

corresponds to the stability limit, while 



K 

V 



1 1 . Z5 and x = -:r 



corresponds to instability, the denominator of the closed loop function, D(s) is 
rewritten. 

D(s) = l+V(s)R(s) 



The V(s) portion will contain the varying parameters and the R(s) portion the 
fixed ("rigid'*) parameters. 

Thus, 



and 



so 



R = 



11. 25 



[s+1)(4 + D 



D = 1+RV 



(4-2) 



(4-3) 



The nominal function, V 

n 

Then the roots of 

D 



are those of 



is taken as unity, i. e. , x = 0, in the expression for V, 



1 + RV 



'V 

n 

V 



n 



0 



V 



0 



(4-4) 
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or 



{sx + 1) + 



11. 25 

(s+ 1){|+ 1) 



= 0 



(4-5) 



In equation (4-5), the expression (4-2) plays the role of in (3-5). Thus, 



and similarly. 



plays the role of 



H 



H 



L = 



V 

n 

V 



= Q 



11. 25 



(s+l){| + l) 



= { sx + 1 ) 



(4-6) 



(4-7) 



-L 



In analogy to the method described in Chapter Three, the expressions for, 

n 

representing the fixed parameter contour and, , the region of criticality 



(here a straight line) representing the varying parameter structure, are plotted. It 
is noted that this varying parameter structure is characterized by the mapping 
over the range of x and for specific frequencies uu . The remainder of the exam- 
ple clarifies this notion. 

The expression (4-7) plots as a straight line for any given uu with its bounds 
defined by the end points on the range of x. For example, 



At x = 0 



At x = 5 



let 


S = = j 


i. e. uj = 1 


then 


II 

•>— 4 


jx + 1 




il 

•>— 4 
> 1 


1 




II 


j5+ 1 
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At X = 0 



At X = 5 



For ^ = 0. 5 : 
then 

V 

^{jO. 5) = j(0. 5) x+1 



V 

(jO. 5) = 1 



V 

(jO. 5) = j2. 5+ 1 



At X = 0 



At X = 5 



For uu = 2: 
then 

V 

(jZ) = j2x+ 1 



V 

(j2) = 1 



V 

^ (j2) = jlO + 1 



The results are plotted in Fig. 4-6 and are represented by the vertical line running 
from 1 +j0 to 1 as UJ-»oo, Actually shown on Fig. 4-6^are the situations 

reflecting ^ = 0.5, uj = 1 and uj = 2. Thus, when 1, the plot may assume 

any point along the vertical line segment from 1 + jO to 1 +j5 depending upon the 
particular value of x. This segment is defined in Fig. 4-6 by the vertical distance 
dimensioned ^=1, and the others similarly. 

, A plot for the FIXED PARAMETER structure, is obtained from the 

negative of (4-6). It is easy to see that 

-L^(jO) = -11, 25 = 11. 25 /l80 = 1 1. 25 /-1 80 



and that the 
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Fig. 4-6 : Dual Nyquist Plot of a Varying-Parameter System 
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lim L-L (jU))] = 0/0 

OJ «5 ^ 



thus providing the end points of the contour. Other points may be found in the 
normal, tedious fashion or by transference from a Bode^ plot. This contour is 
also shown in Fig. 4-6. 

It is (^fficult to see at this point that the "difference vector" drawn from, 
-L^, to the contour (head end at the latter), and connecting identical values of 

u) on each contour, does itself undergo a clockwise rotation if . Without 

some experience viewing the plot, initial interpretation may be formidable; for 
this reason, consider x assuming only two specific values in its range. When 
X = i , the system is stable, and when x = , the system is unstable; these two 
conditions are now studied. Thus, for 



and 



F^(s) 

3 



11. Z5 

(|+l)(s+ 1)(|+ 1) 



Fj(s) 

T 



11. 25 

(|+l)(s + l)(|+l) 



(4-8) 



(4-9) 



there is stability with (4-8) and instability with (4-9). This can be checked quickly by 
counting sign changes in a Routh test.^^^ 

It should be recalled that 



F^(s) 

4 



11. 25 

(-| + l)(s + l)(|+l) 



(4-10) 



represents a system precisely at its stability limit. 

Separation of the varying (in this instance, "fixed**, but arbitrary) and 
fixed parts of the characteristic equation results in (4-5), repeated below. 



nr 



Gardner, M. F. and J. L. Barnes, 
Wiley, New York, 1942. 
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whe re 




Thus, 



( sx + 1) + 



or 



5 * 



for 




ll. 25 

(s+ 1)(| + 1) 



L 

n 



0 



(4-5) 



(f + 1) 



11. 25 

(s+l)(|+l) 



0 



The -L^ plot is the same as in the Fig. 4-6 plot and is reproduced in Fig. 4-7; the 
(*|'+1) plot is again a straight line (this time for varying ^ and fixed x) upon 
which are marked, to the right of the line, the corresponding frequency points. 

It is noted in Fig. 4-7, that the difference vector (shown for uu = 1) never under- 
goes a complete encirclement (360^ or more). At ttJ = 0, the difference vector 
points toward^O^ (due right), while at slightly less than tu = 2 (first intersection 
of -L^ and --^)» the difference vector points toward -90^ (straight down). Some- 
where in the neighborhood of ^ = 3, this vector experiences its largest clockwise 
excursion (approximately -120^), after which it begins a counterclockwise return 
rotation until at = 4 it has almost returned to 0^. As it undergoes a 

counter-clockwise rotation approaching 90^ in the limit. Clearly, a complete 
encirclement does not occur. 

For X = i : 

s . 

The plot is, of course, unchanged. The (-j+1) plot superposes over the line 

s 

describing (-^- + 1) but with different frequency correspondences. These are shown 
by the markings which appear to the left of the vertical line in Fig. 4-7. It can be 
seen from Fig. 4-7 that as runs from ^ = 3 to = 4 the difference vector 
between the two contours undergoes an encirclement. It is clear that whether or 
not encirclement occurs depends upon the relative values of frequency on the -L 

g 

and (—+1) contours. 

A magnification of the situation is contained in the enlarged sketches of 

Fig. 4-8. In Fig. 4-8{a) is shown the hypothetical situation where ^ ^ (at inter- 

V ^ 

-^) = 3. 7 and y (at intersection with “ 4. 7. The effects of 

n 



section with 
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Fig. 4-7 : Detailed Dual Nyquist Plot 
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Fig. 4-8 : Hypotht'tiral Crossing Situations (in the small) 
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encircling the right half of the s-plane upon the mapping of the difference vector 
between -L^ and (sx+ 1) is shown in Fig. 4-9(a) and (b). In a qualitative plot 
of the difference-vector's angular rotation, Fig. 4-10, there is no encirclement 
of the origin and the system therefore has no right-half-plane roots in its char- 
acteristic equation; since the open-loop system is stable this conclusion is per- 
missible. Thus, the system is closed-loop stable. 

A different state of affairs is depicted in Fig.4*ii(b) where ^ (inter- 

n 

section) = 3. 7 and (intersection) = 2. 8. The effect upon the mapping of the 

n 

V 

difference vector is such that. Fig. 4-1 1(a), (b) and the qualitative plot of Fig. 4-12, 
two encirclements occur indicating a system with two closed-loop poles in the 
right-half plane and unstable characteristics. 

In this simple case, it can be inferred that if the value of the frequency on 
the contour is greater than that on the locus of the varying structure, at the 

critical intersection point, then instability occurs. In short, if the locus 

"encloses" the varying locus there is an unstable condition; "encloses" is inter- 
preted to mean: if the frequency at the intersection of the two loci is on 

and if in approaching the point -L (j^.) as the -L contour is traversed in the 

n 1 ^ n 

direction of increasing frequency, the point - y.— ( j ^ J is passed to the right of the 
-L.^ contour (called enclosure), then the system is unstable. Again, there is 
inherent here the asumption of open-loop stability. 

Before this simple example is left, the precise values of the intrinsic para- 
meter ^ (at the points of intersection) are calculated for each contour. This is 
done for completeness of the example. It is not required for the stability test 
under ordinary circumstances since, with practice, encirclement can be envisioned 
from the relative locations of discrete frequency points on the dual Nyquist plot. 

For the stable configuration: 




-L 

n 



-11. 25 
(s+ 1)(| + 1) 



Intersection occurs in this example when 
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(b) 
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Fig. 4-9 



Hypothetical Crossing Situations (in the large) 
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Fig. -1-10 : Exc ursion of iho Difference Vc'ctor 
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Fig. ^-11 : Hypothetical Cross.ing Situation (in the large) 
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Eiccursion of the Difference Vector 
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V 

Re = Re[-LJjuJ)] = 1 

where ”Re" stands for "the real part of". 

Thus, from 




it is found that uu^ = 1. 87 and 0)^ = 3. 75, where the subscripts indicate the fre- 
quencies of first and second intersection respectively. 

It is now clear (since in this simple case the locus of the varying portion is 

the same straight line but with different frequency markings) that, when x = i » 

V 1 ^ 

-L (j3. 75) does not "enclose" — —(}3, 75); when x = , -L (j3. 75) does 

n y V 3 n 

"enclose" (3. 75). Of course, this statement can be checked analytically. Thus, 




(4-11) 



where "Im" signifies "the imaginary part of". For uu = 3. 75, (4-11) evaluates to 
0. 936. But 



Im 



at the point of intersection. Therefore, for x = 



V 





fv ] 


1 

r 





Im = ^ = 0.936 



or 



uu = 4. 68, on ju) at intersection 



V 

and, clearly, ~^(j3. 75) is to the left of the contour indicating a stable 

system. 

For the case when x = -j , similarly, 

V 



Im^(juu) = y = 0.936 
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OJ = 2. 81, on at intersection 

V 

and, clearly, y — (j3> 75) falls to the right of the contour -L^(jo)), indicating an 
unstable system. 

A thumb rule for interpreting the stability criterion in this simple situation 
has now been established. In the original example, Fig. 4-7, an attempt at inter- 
pretation can be made using the additional insight just gained. It is known that the 
intersection of interest occurs at the point -L^(j3. 75), and instability obtains for 
any value of x such that 

uUy (intersection) < (intersection). 

n " n 

V 



In the particular case where the locus of the varying portion is a straight line 
there is the condition that 

Im(sx+1) = Xy. X = 0.936 
n 

But for stability, 

>3.75 

n 

“V 

therefore, since 

. 936 

X = 

^'V 

n 

“V 



for stability 



X 



. 936 
3. 75 



. 250 



Thus, for stability, x < . 25 which reflects precisely the conditions seen in the 
root locus plots of Fig. 4-4 and 4-5. 



4. 3 Multidimensional Variable Parameter Mapping 



The system examined in Section 4. 2 had but a single varying parameter. 



I 





85 



It is for this reason that the locus of the critical point governing stability moved 
along a single curve, in that example, a straight line. When more than one para- 
meter in a system varies (multidimensional parameter variation), the locus of 
critical points assumes the form of an area. In a well-behaved situation, for 
example, a mapping of the boundaries of varying a and 3, Fig. 4-13(a), appears 
as in Fig. 4- 13(b). In this case, every element XY which scans the constrained 
■ range of cc 3 space is inside the mapping's periphery 0^1 2 ^3 ^ called the "hull". 

In such a case the mapping is "complete" and with a suitable generalization of 
Section 4. 2 the stability criterion is necessary and sufficient. 

In the general case, however, end point mapping of the region of variation 

produces results as appear in Figs. 4-14(a) and (b). Here, for a given distance 

along the CX-axis the element X Y maps into a line contained inside the hull; 

a a 

but, further along it is clear that has points which are not contained in the 

hull, e. g. , X^ Y^ , Fig. 4-14(b). Thus, although the end points of the — ^ line 
segments are the images of the end points of the rectangular region in cl 3 space 
the mapping may not be "complete". It is clear, however, that all elements X^Y^ 
are inside the triangle 0^1 ^3^; this "polygon" is called the "convex hull". 
Obviously, the convex hull contains points which are not a result of mapping the 
a 3 variation. It is for this reason that exclusive use of the convex hull provides 
only a sufficient condition for stability, i. e. , for no "encirclement" a system is 
stable but, since an "encirclement"might be of a point or points not part of the 
mapping, an "encirclement" does not necessarily indicate "instability". This 
becomes clear in the examples. It should be noted, however, that a complete 
mapping can always be obtained. One need only plot a few points such as X^ Y^ 
of Fig. 4- 14(b) and the true hull can be sketched quickly. It is shown later that a 
closed form expression (a parabola) can be obtained for the envelope of the lines 
X^Y^ in Fig. 4-14(b); this will have direct application in the human operator 
situation. 

Frequently, the requirement for mapping additional points is precluded by 
the relative location of the fixed and varying contours. For example, assuming the 
situation shown in Fig. 4-15 holds in a similar way for all liJ of pertinence in 
addition to it is obvious that whether 0 12 3 or 0 1 2 2 3 is used 
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Fig. 4-13 : Boundary Mapping Resulting in Complete Hull 
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Fig. 4-14 : Boundaiy Mapping with Hull wliich is Not Complete 

(also shows Convex Hull) 
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Fig. 4 -15 : Hypothetical Variable Parameter Situation 



REAL 



88 



makes little difference. The shaping of is carried out with respect to a 

region common to both. This property does not hold dogmatically when the human 
operator is considered, since being an adaptive mechanism, he will strive to 
assume characteristics in that part of his 'Varying region" which afford stability. 
It is through this understanding, however, that insight may be gained concerning 
the amount of effort the human operator has to expend to maintain control. At 
times it is possible to demonstrate that a mapping of the end points in a 8 space 
yields almost the true hull and is negligibly different from the complete mapping, 
so that practical use can be made of the end point mapping. This is the situation 
with the human operator and the mapping is called "almost^complete". 



4. 4. A Map of the Human- Operator Adaptive Range 

With the human-operator's gain-varying capability and delay time lumped 
into the fixed parameter structure, the varying portion of the characteristic equa- 
tion is givenby (3-12), repeated here 



(1+ sTj^)(l + sTj) 
17“ ^ (1 + sT, ) 

O 



(3-12) 



In order to study the human operator pure delay effect, T, over its entire range of 
variation, several plots for different values of T must be made. Usually, due to 
the limited range of t and the prospect of interpolating middle range effects from 
the end bounds on T, one or two plots will suffice. The plots (of (3-11)) are 

all identical except for the phase lag introduced by the particular value of T used. 

The region of variation for T^^, T^ and T^ is shown in Fig. 4-16, i. e. , 
equation (3-12) is used to map this three-dimenstional volume into the complex 
plane plot of Q(ja)). The region covered by 0 as the triple assumes 

all values on the bounds of the rectangular par allelipiped is shown in Fig. 4-17. 
(Later the completeness of this "end point mapping" is considered and it is shown 
that the end point hull is almost complete as described in Section 4. 3. ) For motion 



along the boundary AB of Fig. 4-16 there is obtained 
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Fig. 4-17 : Mapping of Human Operator Adaptive Range 
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(Tjs+ 1)(0. 7s+ 1) 

^ = fST Ss+l) (4-12) 

but for a fixed u) this is the equation of a straight line in T^. For example, let 
uu = 1, then 






(Tjj+ 1)(0. 7j+ 1) 
(5. 3j+ 1) 



and after some manipulation 



n(ji) 



Tj( 4. 6 + j4. 71) + (4. 71-J4. 6) 
29 



which is recognized as the general vector equation of a straight line 

R = A + tB 

Since Tj has a range of interest, its end points, T^= 0 and Tj=25, define a line 
segment in the space of the mapping ; this line segment has end points given by 
the images of Tj= 0 and Tj=25, i. e. , 
at Tj=0 

^ 6 = 0, 162-jO. 159 

while at Tj= 25 

^ _ 25(4.6 + j4.71) + (4. 71-j4.6) 

29 



n = 4. 13 + j3. 91 



This line is plotted as in Fig. 4-17. It is obvious that lines GF, OF and 

DC map in a similar manner into straight lines in the image (or mapping) space. 
Also, variations in T^, for fixed and T^, map into straight lines; thus, 

OD, GA, FB and EC can be mapped easily. The results of this mapping yield the 
end points in the image space for uu = 1 as shown in Table 4-1. 

These lines are mapped accordingly, Fig. 4-17; their images are shown by 
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TABLE 4-1 



Line Mapped 


^L 




T 

N 


End Points of Image Line 




5. 3 


0 


0. 7 


0. l62-j0. 159 


AB 




5. 3 


25 


0. 7 


4. 13 + j3. 91 




5. 3 


0 


0 


0, 0345-jO. 183 


GF 




5. 3 


25 


0 


4, 61 + jO. 679 




0 


0 


0 


1 


OE 




0 


25 


0 


1 +j25 




0 


0 


0. 7 


1 + jO. 7 


DC 




0 


25 


0. 7 


-16. 5 + j25. 7 




0 


0 


0 


1 


OD 




0 


0 


0. 7 


1 + jO. 7 




5. 3 


0 


0 


0. 0345-j0. 183 


GA 




5. 3 


0 


0. 7 


0. 1625-jO. 1586 




5. 3 


25 


0 


4, 61 + jO. 679 


FB 




5. 3 


25 


0. 7 


4. 12 + j3. 91 




0 


25 


0 


l+j25 


EC 




0 


2^ 


0. 7 


-16. 5 + j25. 7 
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the same letter designations but primed. Thus, faces ODCEO and GABFG of 
Fig, 4-16 map into the triangles and line segment appendages of 
and G^A^B^F^G^ of Fig. 4-17, respectively. 

The mapping of lines AD, GO, FE and BC into Q space can now be con- 
sidered. The most general case, i. e. , T^, T^ 0, is represented by line BC. 
Again, for uu = 1 , 

_ (j25+l)(j0. 7 + 1) 

^ 



-16.5 + j25. 
l + jT^ 



(4-13) 



Equation (4-13) is recognized as the equation of a circle; a proof is given in 
Appendix A. Since three points determine a circle, the arc corresponding to 
(4-13) can be found from the two end points on the range of T and the property 
that as ® the circle passes through the origin, i. e. , 

Q(juu)“*0 as 



Thus, on BC the three points are: 



H 



H 



-(ji) 



Tl=0 



H 






0 



-16. 5 + J25. 7 



Ti^=5.3 



H 

-^(jl) = 4. 13 + J3.9 
o 



TTT 



Churchill, R. V, ."Complex Variable and Applications”, 2nd Ed. , McGraw- 
Hill, New York, I960. 
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Alternatively, the equation of the circle is given by (A-8) where in this case, 
a="l6, 5, b = 25. 7 yielding 

(u+8. 25)^+(v- 12. 85)^ = (15.28)^ (4-14) 

Of course, only the arc of the circle corresponding to the range 0 < £ 5. 3 is 

of interest. In a similar manner, the mappings of lines GO, AD and FE are 
obtained; the information, for uu= 1, is assembled in Table 4-2. 

TABLE 4-2 



Line Mapped 






T 

N 


End Points of 
Image Arc 


Equation of 
Corresponding Circle 




0 


25 


0. 7 


-16. 5 + j25. 7 


(u+ 8. 25)^+ (v- 12. 85)^ 


BC 


5. 3 


25 


0. 7 


4, 13 + j3. 9 


=(15. 28)^ 




0 


25 


0 


1 + j25 


(u-i)^+(v- 12.5)^ 


FE 


5. 3 


25 


0 


4, 61 + jO. 679 


=(12.53)^ 




0 


0 


0 


1 


(u--) +v 


GO 


5. 3 


0 


0 


0, 0345-jO. 183 






0 


0 


0. 7 


1+jO. 7 


(u --|)^ + (v -0. 35)^ 


AD 


5. 3 


0 


0. 7 


0. 1623-jO. 1586 


=(0. 6105)^ 



It should be observed that, interestingly enough, the equations of the 
circles which represent variation in T^ are not themselves functions of T^, 
although the end points of the arc segments reflecting the range of T^ surely 
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are. Hence, as long as the ranges on T^. and remain unchanged so do the 

equations of the circles. This property proves to be extremely useful in this 
study later on, and will provide an opportunity for obtaining insight and intuitive 
understanding of the effect of human operator adaptation upon any given system. 

The end result of Tables 4-1 and 4-2 is shown plotted in Fig. 4-18 and 
represents the mapping (3-12) plotted for the single frequency uu = 1. A sufficient 
number of such plots must be made for selected frequencies to permit determina- 
tion of whether the locus of (3-11), the difference vector between -L and 
H 

o 

— = Q, undergoes encirclements of its origin or not. From this, stability 
o 

information may be inferred. Since the human-operator transfer function is valid 
only over a limited frequency range, the plots are not infinite in number; addition- 
ally, their tendencies are obvious and a few judiciously chosen discrete frequencies 
are seen to suffice for the purposes of this study. Several points are worth noting 
at this time: 

1. The human operator's pole variations map into straight line segments 
while his zero variations map into arc segments. 

2. The formation of a largest polygon (with some arc boundaries) around 
the bounds of the region in Fig. 4-18, i. e. , formation of the convex hull, defines 
an area of possible criticality re stability. 

3. It must be determined if the end-point mapping results in a complete 
mapping, i. e. , every point of the mapping Q pnust be a point of the end-point 
mapping if the mapping is to be complete 

4. If the hull is not complete, rather than use a convex hull, can the ’’end 
point” hull be shown to be almost complete ? 

5. The situation has changed progressively through the usual Nyquist cbn- 
dition where the critical point is -1 + jO, to a one parameter variation where the 
locus of critical points is a curve, to a final condition where the critical points 
are defined by an area Q . 

It is shown that Point 4. pertains, but that the hull of the complete mapping 
can always be found and may be described analytically; however, the end point 
mapping of Fig. 4-18 differs so little from the complete hull (indeed, very little 
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Fig. 4-18 : Typical Human Operator Variable Parameter Mapping (o)= 1 ) 
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from the convex hull also) that for all practical purposes it may be used. 

The interpretation of Point 5. is: if there is some area of Q unenclosed 



meter values and stabilize the system. The parameter values necessary to ac- 
complish this give insight into the amount of effort the human operator must 
expend to retain control. If the entire Q region is enclosed, the human operator 
cannot control under any circumstances (for the capability ranges assumed and 
plotted). If no portion of the Q region is enclosed, there is stability no matter 
what the actions of the human operator (within the sense of the transfer function 
description). Later, examples and discussion clarify these notions. 

For the study to proceed, the regions of varying parameters, fi, must be 
plotted for selected frequencies. Generally speaking, the human-operator trans- 
fer function has validity only over a range of frequencies, uu < ZO. It can be seen 
that beyond an tju = 3, the plots offer little additional information and their trend 
is obvious. The region Q is plotted therefore, for (ju = 0.5, 1.0, Z. 0, 3.0 and 
also for an uu = 10; in order to do this it is necessary to derive the mapping 
relations for the various bounds of Fig. 4-16 as functions of uu . Two such 
derivations are given below; the first for a human-operator pole and the next for 
a human-operator zero, the remainder are obtained in a similar way. 



The human- operator pole variation, Tj , along AB of Fig. 4-16, is found 



by the contour then presumably the human operator can assume such para- 



from 



Q 



(Tjjju+ 1)(0. 7jJU+ 1) 



(5. 3j UU+ l) 



After some manipulation. 



n 



(Tjju)+ D(1 + 3. 71 /- 4. 6 jju) 



(4-15) 




Thus, at Tj= 0, 



Q 



1 + 3. 71 u)^ - 4. 6 jut) 
1 + 28ou^ 



(4-16) 
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and, at Tj=25, 



^ _ 1 + 118. 71uu^+ j uu(92, 8uu^+ 20. 4) 
1 + 28uu^ 



(4-17) 



The mapping of AB is a straight line between the points obtained from (4-16) and 
(4-17). For variation in the human- operator zero, T^, along BC in Fig. 4-16, 



n _ (25juu+ 1)(0. 7juu+ 1) 

(TJUU + 1) 



(-17. 5uu 4 25 juu4 0. 7 jou4 1) 
1 + T^j (ju 



If the complex number n(juu), is represented by u4jv, then 

jv u - JV 



1 + T juu = ^ ^ ~ 

Lj^ U + IV * U - IV 



Z Z 

U-17.5UU u + 25. 7juuu-jv+17. 5juuv+25. 7 uuv 

2 , 2 

U + V 



The real parts may be set equal to each other (the imaginary parts permit recovery 
of Tj^ for given uu, Appendix A): 

, u- 17. 5uu^u + 25. 7 uuv 



If the square is completed, one obtains 




(4-18) 



Relation (4-18) is the equation of a circle, an arc segment of which, BC maps into, 
as a function of frequency. 
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At T^=0: 



= 1- 17. 5 1 ^5. 7j 



(4-19) 



and at 



I + 1 18. 6:*^+ j ; (9^. 8 /+ ZO. 4) 
1 + Z8: 



(4-ZO) 



In a similar manner, corresponding relations are obtained for the remaining 
bounds of Fig. 4-16; these are displayed in Talkie 4-3. 

It is now possible to plot the mappings of the parameter variations for 
particular frequencies; these plots are shown for ; =0.5, 1.0, Z. 0, 3.0 and 10.0 
in Figs. 4-19, 4-ZO, 4-Zl, 4-^iZ and 4-Z3, re specti\ely. The question of the com- 
pleteness of these end point mappings may now be considered. As is clear from 
Fig. 4-Z4 the end point mapping is not significantly different from the convex hull. 
This is particularly true in the light of system plots of ” contours for typical 
physical plants, as can be seen in later examples. The difference between the end 
point hull and the true hull (complett* m«ipping), howev^er, is even less than between 
the end point hull and the convex hull, in the complete mapping, a small parabolic 
arc extends between segment and D as well as between G A and 

For the purposes of this study, this situation may be ignored. Details and 
a proof arc‘ contained in Appendix B. 

One further point of some signifi Ctince relating to the mapping should be 
noted. As stated above, face ODCEO of Fig. 4-16 maps into the triangle plus 
line segment O^D^C^E'O", shown in Fig. 4-17; similarly, face GABFG maps 
into G^A^B^F^G\ The foriner is associated with a while the latter 

corresponds to a 5. 3; this value represents the upper bound on the range for 

used in this study. It is interesting to observe, and useful later in the physical 
interpretations, that as T^ varies from 0 to 5. 3, face ODCEO "rotates" 
into coincidence with face G'a'b'f'g'; the "rotation" occurs by virtue of con- 
strained motion of the significant apex points of the "triangle" along the circular 
arc segments of Fig. 4-18. 

It may be preferable to reflect upon the effects of a two parameter mapping 
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TABLE 4-3 





End Point Value 


General Equations 


Line Mapped 


of Varying Parameter 


of the curve 






0 


1 + 3. 71uu^-4. 6j uu 








1 + ZSi/ 


AB 




25 


1 + 118. 71u)^+juu(92. 8uu^+20. 4) 






1 +28u)^ 






0 


1 - 5. 3juu 
1 +28u)^ 


GF 


= 


25 


1+132. 5tu^+ 19. 7j(i) 






1 + 28uu^ 






0 


1 


OE 




25 


1 +25juu 






0 


l + .7ju) 


DC 












25 


1 - 17. 50J +25. 7juu 






0 


1 + 132. 5uu^+ 19. 7ju) 
1 +28uu^ 


FB 


Tn = 


0. 7 


1 + 118. 71uu^+ juu(92. 8o)^+20. 4) 
2 








1 + 28o) 






0 


1 


OD 


Tn = 


0. 7 


1 + . 7ju) 
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TABLE 4-3 (cont'd. ) 



Line Mapped 


End Point Value 
of Varying Parameter 


Uenerai Equations 
of the curve 






0 


1-5. 3juu 


GA 


T = 




1 +28(1? 




0. 7 


1 +3. 71(1)^- 4. 6ju) 








1 + 28UU 


EC 


II 


0 


1 + 25j uo 




0. 7 


1-17. 5u)^+ 25. 7ju) 






0 


1-17. 5u?+ 25. 7ju) 


BC 


^L = 




1 Cu-(. 5-8. 75(11 )] +Cv-12. 85 (u3 
I =(. 5-8. 75(U^f+(12. 85(1))^ 



, 1 + 118. 7l/+ j(ju(92. 8(ju^+ 20. 4) 

1 + zqJ' 

T5 1 




5. 3 
~ 



AD 



T 



L 



0 



1 - 5. 3jU) 

1 +28w^ 

1 + . 7ju) 

Cu-|]^[v-.35uu]^=(if+ (.350,)^' 

1 + 3. 71gj^-4. 6ju) 

1 + 28uj^ 

1 + 25ju) 



FE 



[u--i]^+Cv-12. 5u)]^ = 



+ (12. 5uu) 



2 



5. 3 



1 + 132. 5(ju + 19. 7j(ju 

1 

1 +28(ju 
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Fig. 4-19 : Variable' Paramclcr Mapping 
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Fip.4-Z0 : Variable Parameter Mapping 
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Fig. 4-21 : Variable Parame*ter Mapping 
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Fig. 4-22: Variable Parameter Mapping 
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Fig. 4-Z3 : Variable Parameter Mapping 
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Fig. 4-24 : The Convex Hull 
(heavy lines) 
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before continuing on to the practical results, derived from the three parameter 
mapping, which are described in the next chapter. To this end, and also to offer 
further stability interpretation practice, consideration of Appendix C may be 
appropriate at this time. 

4. 5 Summary 

After the analysis of a simple system containing only one variable para- 
meter, the situation of multi-parameter variation is considered. With the famili- 
arity gained from the single parameter example the generalization to a multiple 
parameter problem follows more easily. In particular, the three-varying para- 
meter problem (one zero and two poles) is considered in detail, since it relates 
directly to the three time constants in the human- operator transfer function that 
are within the adaptive capacity of the human controller. It is shown that these 
parameters, forming a type of variable structure, can be used as a mapping func- 
tion in the complex plane to display the boundaries of human- adaptive capacity. 
These boundaries contain a collection of ’’critical points”, akin to the -1 + jO 
point on the usual Nyquist diagram; when the plant Nyquist contour is plotted in 
conjunction with these bounded areas (at particular frequencies) interpretation of 
system stability characteristics is possible (Chapter Five). Several plots of such 
boundaries are presented for important frequencies in the human range. 
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CHAPTER FIVE 

APPLICATIONS OF THE THEORY 



5. 1 Introduction 

The groundwork which is laid in Chapters Two, Three, and Four now 
makes it possible to consider actual physical systems which are operated by 
human controllers. If the empirical results can be predicted (or implied) by the 
approach outlined in the previous chapters, then this technique for viewing the 
response characteristics of systems with human operators in the loop is justified, 
and difficulties with human controllability of such systems can be anticipated via 
the theory. 

Of the several experimental studies performed with human controllers, . 

two are particularly well suited to the purposes of this investigation. The Hall 

study^^^ is notable for its comprehensive nature and for the use of a non-trivial 

controlled element. The predictions obtained from the theory presented in the 

previous chapters can be checked for their physical validity using the results of 

(Z) 

the Hall experiment. The Smith data is derived from human attempts to con- 
trol an unstable system; with these empirical results, the effect of open-loop 
instability upon the theory also may be compared with actual physical accomplish- 
ment. 

5. 2 The Hall and Smith Studies 

In 1956 and 1957 a simulator study of the human pilot was conducted by 
Ian Hall at Princeton University. The intent was to investigate the effects upon 
the pilot of various controlled-element-longitudinal dynamics. A very excellent 
analysis by McRuer et al, uses the data of that report as a basis for developing 

Hall, I. A. M. , Effects of Controlled Element on the Human Pilot, Princeton 
Univ. , WADC Tech. Rept. No. 57-509, Aug., 1958. 

Smith, R. H. , An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
On the Limits of Manual Control, IEEE Trans, on Hum. Factors in Electronics, 
Sept. , 1963). 
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a "flying qualities" theory. Briefly, the McRuer approach follows the "adjust- 

ment rules" procedure outlined in Section 2. 1 of this study. Another experimental 

study conducted at the University of Cincinnati by Ralph Smith probed for the 

( 2 ) 

human limits in controlling certain open-loop-unstable systems. ' The results of 
both these studies are used to significant advantage here; indeed, the empirical 
evidence of Smith^^^ and Hall^^^ furnished the experimental justification for the 
results of this study. 

5.3 The Hall Study 

The experiments were conducted in a mockup cockpit simulating a Navion 
aircraft. The subjects were two Navy pilots, one with considerable jet time 
(Pilot I) and the other with a preponderance of propeller aircraft experience 
(Pilot II) . The seemingly subjective pilot rating scale of the Hall study used the 

(3) 

following definitions (quoted): 

'» GQQD (G): aircraft response very suitable for the given tracking task. 

ACCEPTABLE (A): possessed of satisfactory flying qualities, but less 

easy to control than a "good" configuration; aircraft response either 
faster or slower than that which the pilot visualizes as most desirable. 
POOR (P): possessing characteristics which require considerable con- 

centration to accomplish the given tracking task; configurations 
are either quite touchy, unresponsive, or poorly damped. 
UNACCEPTABLE (U): quite unsuitable. 



The author discussed the study with Pilot II in November, 1964. 

McRuer, D. T. , I. L. Ashkenas, C. L, Guerre, A Systems Analysis View of 
Longitudinal Flying Qualities, WADD Tech. Rept. 60-43, Wright Air Dvlpmt. 
Div. , J an. , I960. 

Smith, R. H. , An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper. 
On the Limits of Manual Control, IEEE Trans, on Hum. Factors in Electronics, 
Sept., 1963.) 

Hall, I. A. M. , Effects of Controlled Element on the Human Pilot, Princeton 
Univ., WADC Tech. Rept. No. 57-509, Aug., 1958. 
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The subjects differed only between P and U. Pilot II tended to classify aircraft 
which were unsuitable for the given tracking task, but could conceivably be flown, 
as P, reserving U for configurations which he considered as unsafe to fly under 
any circumstances, while Pilot I tended to classify all configurations with which 
tracking was very difficult as U.” Due to the use of plus , "somewhat more than", 
and minus , "somewhat less than" modifiers on the G, A and P ratings the resul- 
tant number of classifications was ten, i. e. , A + -*4, G “*2, G + 1, etc. Those 

familiar with the military aircraft rating scheme are aware that the standard 
Cooper scale also has ten gradations; they are not exactly the same. It may 
appear that the entire rating system is seriously suspect by virtue of its apparent 

subjectivism. The consistency of pilot opinion, however, is remarkably good and 

( 2 ) . 

the variation between pilots astonishingly small. The inference, of course, is 
that objective factors are coming to play and it is there that the prospect of a pos- 
sible quantitative theory for handling qualities arises. 

The experiment was conducted in the compensatory arrangement as shown 
in Fig. 2-1 with the plant (aircraft) dynamics having the generic form (short-period 
approximation): 



5(1 +0. 6s)0J 



G = 



n 



^ s(s^+2Cuu s +(JU^) 



(5-1) 



n n 

Hall used twenty combinations of Q and to test pilot response. Of 
interest here is a representative sample covering the spectrum of pilot opinion. 
The configurations designated by Hall as C30, C25, C32, C31, and C21 are to be 
studied, Table 5-1. 

The following pertinent performance comment is a quote from the Hall 
report. "As a measure of tracking efficiency. . . an average absolute value 



Cooper, G. E,, Understanding and Interpreting Pilot Opinion, Aero. Eng. Rev. 
Vol. 16, No. 3., Mar., 1957. 

Ashkenas, I. L. and D. T. McRuer, A Theory of Handling Qualities Derived 
from Pilot- Vehicle System Considerations, Aerospace Engineering , Feb., 1962. 
(3) 

Hall, I. A. M. , Effects of Controlled Element on the Human Pilot, Princeton 
Univ. , WADC Tech. Rept. No. 57-509, Aug., 1958, 
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TABLE 5-1 



Configuration 


Opinion Rating 


Mean Absolute Tracking 
Error [c] (degrees) 


c 

number 


(JU 

n 


c 


Pilot 

I 


Pilot 

II 


Average 


Pilot 

I 


Pilot 

II 


Average 


C30 


2. 51 


1. 0 


3 


1 


2 


. 85 


.64 


. 75 


C25 


1. 57 


1.0 


6 


3 


4. 5 


.90 


. 77 


. 84 


C32 


3. 77 


0. 35 


4 


6 


5. 0 


. 83 


. 77 


. 80 


C31 


3. 77 


0, 2 


5. 5 


8 


6, 75 


. 84 


, 82 


. 83 


C21 


1. 57 


0. 2 


10 


10 


10 


1. 30 


1. 04 


1. 17 



criterion was preferred to a root mean square error, since it was believed that 
the occasional large errors, which essentially decide the value of RMS error, 
more often were due to lapses in concentration rather than conditions associated 
with longitudinal dynamic characteristics, . .In discussing average error it is 
interesting to keep in mind that without any pilot output this average would be 1. 36 
degrees. . . due to longitudinal. . . excitation signals. . . Average longitudinal errors 
of less than 1, 36 degrees indicate that the pilot is doing some good— something 
which must not be taken for granted, " 

Additionally, the system excitation was random appearing with Gaussian 
amplitude probability distribution which resulted in a realistic target course, 
specifically, the input was white noise cut off at 18db/oct after 1 rad/sec. It is 
noteworthy that "attempts to increase the cutoff frequency to 2 rad/sec. and 4 
rad/sec. had resulted in loss of the pilot's faith in his ability to ever catch the 
target. 

In the case of C21, a Routh test shows this configuration to be closed- loop 
unstable (without a human controller in the loop). 



( ILallj I. A. M. , Effects of Controlled Element on the Human Pilot, Princeton 
Univ. , WADC Tech. Rept. No. 57-509, Aug., 1958, 
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5(d^{ 1 + 0. 6s) 

1 +F(s) = 1 + 5 -: (C = 0. 2 ,( 1 ) =1.57) 

s{s +2 s + uu ) 

' n n 



D(s) = s^+ 0. 628 s^+9. 84s + 12. 30 = 0 



and 



3 

s 

2 

s 

1 

s 



1 

. 628 

-(12.30-6. 18) _ , 
7^28 



9. 84 
12 . 30 



0 

s 



-( 6 . 12 ){ 12 . 3 ) 
n — 



+ 72 



the two changes of sign in the first column signify two positive roots and therefore 
closed-loop- system instability. The human operator* s task is clear if any attempt 
is made to control this system. 

Next, the closed-loop-stability characteristics of C 32 in the Hall report 

are determined. 2 

5 uu ( 1 + 0 . 6 s) 

1 +F(s) = 1 + y (C = 0 . 35, 0 ) = 3.77) 

s(s + 2 C(i) s + 11 ) ) 

n n 

.•.D(s) = s^+2.64s^+56.8s+71 = 0 

The Routh array is given by 

1 56.8 

2.64 71 

.(.89){71) - + 



3 

s 

2 

s 

1 

s 

0 

s 



since there are no sign changes, there are no positive roots, and the system is 
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closed-loop stable (without the human operator). It should be expected that less 
human effort is needed to control this system. 

In plotting these systems in conjunction with the variable parameter map- 
ping representing the human adaptive range, one hopes to see the -L (j uu) contour 

n 

for CZl cross the Q region in a less acceptable manner than the plot of -L» ( j uu) 
for C32 crosses it. The term "a less acceptable manner" may not be adequately 
defined at this point, but the two extremes of stability are clear, i. e. , the case 
of stability regardless of the pilot* s presence, or the other extreme of instability 
irrespective of the pilots attempts. 

The configuration C21 and C3Z are plotted in Figs. 5-1 and 5-2, respec- 
tively. The solid line through the specified frequency points is merely a suggested 
locus and is offered mostly for visual ease. As indicated in Chapter Four, the 
important properties of each system are interpreted from the characteristics of 
the difference vector drawn between all frequencies of the system locus, 
and corresponding frequencies in the variable parameter mapping, Q. The dif- 
ficulty, of course, is just what sub-region, set of points, or point is used in the 
large expanse of Q ? It is recognized that if a point in the space of the Q map- 
ping is selected this is equivalent to selecting a set of values for the triple 
(Tji^, Tj, Tj^). It might seem that this amounts to using a specific set of parameter 
values and working with a nominal transfer function, so the human adaptive range 
is not really used in the study. This is rather a pedantic point of view. It is seen 
that all the possible capabilities available to the human (by way of parameter 
variation) are clearly displayed in the plots and the net effect of choosing any set 
of them (a triple) is equally visible. The usefulness of such information is un- 
deniable, particularly when experience and experiment allow correlation of cer- 
tain of the sub- regions of the Q space with human effort. Such implication is 
suggested even by the few examples viewed here. In any case, selecting a set of 
parameter values corresponding to a point in Q space, or even a cluster of para- 
meter values (Tji^, T^, T^^), each over a tiny range whose image yields a sub- 
region in Q space, is really no more than the human operator is forced to do. 

The human inability, or reluctance, to change parameter character (break fre- 
quencies) discontinuously, corresponds to a constraint upon jumping about in Q 
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Fig. 5-1: Plot ol Conl'iguration C21 



116 



Kcj) 




Fig. 5-2 : Plot of Configuration C32 



space once some "neighborhood” or sub- region of that space has been chosen. As 
to how such a neighborhood or sub-region is chosen, some insight regarding this 
devolves from the examples along with the inherent assumption that the human 
operator will always try to do his "best". It is not unreasonable to associate 
"doing one*s best" with the maintenajice of stability, which is the prime require- 
ment of any system. Once in this sub- region it seems entirely proper that signifi- 
cant departures from it do not occur and interpretations as to whether any portions 
of this sub-region are "enclosed" or not by -L^ are possible. 

From (5-1) is obtained the FIXED portion of the system, the negative of 



5(1 + . 6s)uu^ 
n 

s( s^+ 2 Cou s +uu^) 
n n 



(5-2) 



(5-2) is plotted to obtain -L (180° phase angle is added). Since the delay time, 

e"^^ , and human gain, K , are lumped with (5-2) it should be borne in mind that 

^ - T S . 

Fig. 5-1 actually represents a situation where ® i. e. , t = 0, K^=l. 

The significance of these conditions comes to focus shortly. Equation (5-2) may be 

written 2 

-5(1 + 0. 6juu)uu 

-L = (5-3) 

^ j uu((D - UU + 2 Cut) j uu) 

n n*' 



Equation (5-3) is rationalized to give 

10 uu^(l +0. 6j uu)+ j5 uu^ (D(tiJ^- uu^)( 1 + 0. 6j (d) 

, _ n n n (5-4) 

n n 

The imaginary part of (5-4) is set equal to zero 

uu^(6Ca> -5)+5uu^ = 0 
n n 

which for C21 (C = 0. 2, u)^= 1. 57) gives a real axis crossing at (i) = ±1.99. For C32 
the solutions are imaginary, reflecting no real-axis intercepts. The conditions 
portrayed in Fig. 5-1 and 5-2, respectively. 

Transparent overlays of Figs. 4-19 through 4-23 are enclosed in an enve- 



are 
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lope fastened to the inside back cover of this study. These overlays, which are 
universally applicable, for the given scale, facilitate considerably the determina- 
tion of system stability characteristics. Detailed consideration of configuration 
C21 follows, 

A glance at Fig. 5-1 and the general trend of the Q mappings for the fre- 
quencies given (uu = 0. 5, 1. 0, 2. 0, 3. 0, 10. 0), indicates that system stability is 
determined by the relation between the plot and the overlays for uu= 2. 0 and 

3, 0, This is obvious, since the plot leaves the Q region of criticality for 

thus implying that a condition of instability may pertain. It is not yet clear 
how to take the difference vector between -L^(j2) and 0(j2), i. e. , it could con- 
ceivably be connected from -L^(j2) to any point within or on Q(j2), and for in- 
stability it must be demonstrated that a clockwise encirclement must occur as 

, irrespective of this vector’s selected orientation (all the Hall configurations 
were open-loop stable). Due to the time invariance assumptions, however, once 
a settled condition of operation has been adopted by the human controller, any 
’’selected orientation” of the difference vector defines the sub- region in space to 
be used for the remainder of the stability analysis. One notes, for example, that 
if the human operator assumes a set of parameters near point of Fig. 4-21 
(corresponding to a neighborhood of max T^, max T^, max T^), there results, 
from Fig, 5-1, an encirclement by the difference vector and instability. If the 
point F^ is assumed, the human can ’’hold out” longer before enclosure, but 
eventually an encirclement again occurs. Unfortunately, the prospects of being 
at F^ are not encouraging since, although the point corresponds to max (lead) 
which is obviously needed, it also implies T^= 0 which, giving some credence to 
the physiological correlations, corresponds to no neuromuscular lag. It is soon 
discovered that if there is any hope for system stabilization or controllability, the 
human operator must assume parameter values somewhere between arc 
and arc and in the vicinity of the origin of Fig. 4-21, also Fig. 4-18 (detail). 

Particularly is this the case when the human time-delay is imposed upon the sys- 
tem. Thus, the characteristic of C21 projecting beyond the Q region is further 
aggravated by the introduction of the time lag, i, e. , the point -L^(j2), at approxi- 
mately 5, 0 on the real axis, rotates with that distance as a radius arm, a number 
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of degrees proportional to the time delay. For example, if T=.15 is assumed, 

, 1 r -j ^ T / • 

the angle of e (its magnitude, of course, is unity) for uu= 2 is (2)(. 15) = 

. 3 radian = 17. 2 , taken in a clockwise direction (lag), and similarly for other 
pertinent frequencies; clearly, conditions are not made better. It might be added 
at this point that the effect of the e ^ term should not be considered as adding 
ever increasing phase, thus further complicating the plot; since the human- 
operator transfer function is not valid beyond a maiximum frequency of uu ^==5^ 20, 
the corresponding phase angle contributed by this frequency can be regarded as an 
upper bound. 

Since C21 is more than a hypothetical situation, it is possible to check the 
implications of the above analysis. A ’'best fit" human-operator description was 
calculated by Hall from the empirical data and a plot of this function can be made. 
The Hall fitted transfer function was (Pilot II): 



H 

o 



0. i(i + |)^ 



15s 



(5-5) 



Several properties of (5-5) should be evident. First, this is one of the cases where 
Hall required a double lead; this characteristic is not accommodated by the gen- 
eralized human-operator transfer function assumed for this investigation. Equation 
(2-1). Nevertheless, the description is analyzed in order that the curve repre- 
senting the Hall fitted transfer function can be observed relative to the regions of 
criticality which have been used in this study. To be noted are the value of human 
gain assumed, i. e. , = 0. 1, and time delay tolerated, T =. 15; the system 

dynamics, -E^, are plotted for ^ although here it is KNOWN that actually 

was 0. 1. The effect of this gain reduction must be imposed upon the system 
dynamics as indicated previously; in an actual study, where the approach is from 
the direction of synthesis, the plots imply that reduced human gain is a necessity 
for this man-machine system if there is to be any chance of controllability. 

Equation (5-5) represents those characteristics assumed by a real life 
human operator and when mapped properly, indicates those points in Q space 
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actually assumed under the task conditions imposed by C21. Table 5-2 displays 
the pertinent information for 

Q = (5-6) 

(1 + 0, 5juu)^ 



TABLE 5-2 (C21) 



uu 


0 


. 5 


1 


2 


3 


7 


10 


n(juu) 


ik 


. 94 /-I 8 . 8 


8zl-34. 8 


. 55 /- 54 . 6 


. 39 /- 6 I. 6 


. 17 /- 52 . 6 


. 01 /- 4 I. 8 















The last row in Table 5-2 displays the values of Hall's "best fit" transfer 
function for several pertinent frequencies. A plot is not required to see that the 
human controller did, indeed, assume those characteristics previously implied 
by the analysis; furthermore, it is obvious that his capacity is being strained 
severely just to maintain control of the C21 system. That is to say, precious few 
points in the variable mapping space 0 are not enclosed by the plot of Fig. 5-1, 
particularly after introduction of the time-delay effects. It is also interesting to 
observe that the double lead does not impose any gross effect upon the experi- 
mentally determined 0 that is not already incorporated in the generalized map- 
ping (3-12). Additionally, it is noted that reduced gain. Fig. 5-1, would ease the 
instability tendencies considerably. 

Since the dynamics of C21 were rated "unsatisfactory" by the pilots in the 
Hall study, it is interesting to look in detail at an analysis similar to that above 
for C32 , to which a rating of "acceptable" was assigned. In Fig. 5-2, the plot 
for C32 of the Hall study, it is clear that controllability should not be as difficult 
a task as presented previously by C21. It should be obvious what the trend is for 
the plots of Figs. 4-19 through 4-22 (the transparent overlays) as the frequency 
continues to be increased beyond UJ = 3, Fig. 4-23. The plot for UJ - 3, Fig. 4-22, 
is already so large as to preclude displaying the entire hull upon a page. For this 
reason, only the region of greatest significance is displayed in the plot. Thus, 
observation of the characteristics of C32, Fig. 5-2, should indicate that the human 
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controller need not reduce his to almost zero while maximizing his lead, thus 

operating in the neighborhood of the origin of Fig. 4-22, as was required in C21; 
but, that there is room for assuming parameter values in Q that yield a critical 
sub-region in the neighborhood of 3. 5 + jl, Figs. 5-2 and 4-22. If the mathe- 
matical modeling here represents in some small way the actual physical situation 
that obtains, it is to be expected that the human controller would feel less vitriolic 
about this system. 

For the system dynamics of the Hall report, C32, the "best-fitted” trans- 
fer function for the pilot was (Pilot II): 



H 

o 






(5-7) 



Again the value of gain actually assumed by the pilot was less than that used in the 
plot of Fig. 5-2; the influence of delay time T =0. 15 seconds also is not included; 
both facts must be incorporated in the final interpretations. The associated map- 
ping function for Q is 



O 3 (l + ^>5juj) 
(1+0. 77juu) 



(5-8) 



The form of equation (5-7) may seem to imply a 0 which would be physio- 

logical nonsense. Frequently, experiments connected with tracking studies and 
the servomechanism approach have incorporated some of the minor lags into the 
delay time wherever such a technique did not prejudice the data. The pertinent 
data for (5-8) is given in Table 5-3, 

TABLE 5-3 (C32) 



uu 


0 


. 5 


1 


2 


3 


7 


10 


n(juu) 


1 


1. 29+ j. 75 


o 

-1- 

00 


2. 59 + jl 


2.9+j. 82 


3. 18+j. 40 


3. 22 + j. 29 



Try 



McRuer, D. T. and E. S. Krendel, The Human Operator as a Servo System 
Element, Journ. of Franklin Instit, ,Vol. 267, Nos. 5 and 6, May and June, 1959. 
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The results of actual human control shown in Table 5-3, viewed in con- 
junction with Figs. 5-2 and 4-22, again bear out the previous implications of the 
analysis. It is seen that the human controller did indeed assume these values 
for his parameter characteristics which the two plots of Fig. 5-2 and 4-22 appeared 
to indicate. Since not as large a value of T^ was required to stabilize this sys- 
tem, the human operator did not assume any larger value than needed when the 
time lag is considered, and since experiment has shown that generating lead 
(anticipation or derivative error control) is regarded by the operator as a strain, 
the system C32 was regarded as superior to C21 and rated ’'acceptable'*. Again, 
it is obvious from Fig. 5-2, that reduced gain will tend to stabilize the systemvthis 
was recognized empirically by Pilot II in the experiment as demonstrated by the 

gain actually assumed, K =0. 5. 

h 

The plots for system configurations C30, C25, and C31 of the Hall study 
are shown in Figs. 5-3, 5-4, and 5-5 respectively. It is noted that the shape 
of the curve for the particular system being considered, as well as the relative 
location of its frequency markings, lead to estimates of pilot-derived handling 
evaluations with orderings equivalent to those appearing in Table 5-1 under the 
column Opinion Rating (average). A quick check of Figs. 5-3 through 5-5 and 
their relation to the mappings on the transparent overlays corroborates this point. 
Thus,C30 is regarded as a quite good "aircraft”, with C25 and C31 rated lower, 
in that order. The significance of these facts and their ramifications as discussed 
in detail for C21 and C32 follow directly. 

Any person conducting research in the field of control systems and the 
human operator soon encounters a preponderance of studies, published and un- 
published, using a nominal value for T^^ of approximately 2. 0; indeed, many 
experimenters feel that human-generated lead cannot exceed this figure. It is 
informative if the effects upon the stability hull of assigning such a bound to the 
range of Tj^ are considered in some detail. As in Chapter Four changing the 
upper bound on the range of T^^ has no effect upon the associated circles them- 
selves (in the complex plane), but only upon the extent of the arc segments. Thus, 
for a 2, one expects to see smaller arc lengths than for a T^^= 5. 3, which 

has been used up to now. It is possible to infer already the effect upon the sta- 
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Fig. 5-3 : Plot of Configuration C30 
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Fig. 5-4 : Plot of Configuration C25 
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Fig. 5-5 : Plot oi Configuration C31 




i 



126 



bill zing ability of the human operator without actually plotting the new specifications. 
One notes that points 5 6 8 7 corresponding to face 5687 in three space, Figs. 

4-18 and 4-16 respectively, do not move along their circular loci to the extent 
previously, when 5. 3. Thus, line 7 '- 8 ' makes a steeper angle with the 

real axis and there are many more systems with amplitude-frequency character- 
istics which leave the hull and tend toward instability. It is interesting to note 
that line 7 - 8 relates to a T^= 0 which, as mentioned above, is usually not 
possible (unless the effects of can be incorporated into an expanded t); 

.hence, for T^=0.7. the stability criterion becomes associated with line 5^-6^ 
and further difficulties with stabilization are implied. It is in this way that the 
designer is given considerable insight into the effects of human- operator para- 
meter variations upon overall system stability and even, perhaps, handling qual- 
ities. For example, with maximum lead in, attention may be confined to face 
5687 corresponding to the augmented triangular mapping 5 ' 6 ' S' 7 \ Clearly, the 
outer bound of this ’’sub-hull" most conducive to stabilizing a system is line 
7^ - 8^ i. e. , relating to max and min T^. It is plain to see now why all 

previous examples drawn from empirical results had stability determination 
dependent upon points clustered in a very reduced area of the large region Q , 
i. e. , this reduced area corresponded to those extreme values of possible para- 
meter variation most conducive to stability. The consequences of various modi- 
fications in human- operator capability are all quickly discerned with this tech- 
nique, i. e. , the region of criticality for stability, if T^= 0, is the mapping of 
face 1375 which is quickly found to have a hull 1^3^7'5^ (with two arc segments), 
and the ramifications of this state of affairs are rapidly seen. Of course, if the 
ranges of Tj and T^ must be altered, the hull in the complex plane is changed, 
but the new points may be found expeditiously from relations such as (4-12) in 
suitably modified form. The effects of varying Tj and Appendix B.2, and 

the characteristic that the T^ circles always pass through the origin, can be 



Appendix C contains an example of the use of the degenerate hull in a two para 
meter mapping problem as well as an instructive demonstration in stability 
interpretation. 
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used to good advantage in plotting any modified range on the parameters 

Tl- t,. t^. 

5. 4 The Smith Study 

In September, 1963 Ralph Smith reported results from an experimental 
investigation attempting to define the limits of manual control for unstable- 
second order controlled elements, ^ ^ There are but two results given! 

[l] C=-.70 '<■ =2.5 K = 300 (5-9) 

n p v-’ 7/ 

Cz] C = -. 35 0)^ = 7. 2 Kp = 1200 (5-10) 

Both of which are for the controlled element having the form 

K uu ^ 

G = E_J1 (5. 

P ? 7 

s s + 0 ) 

n n 

To plot Set (5-9), from (5-11) one obtains 



G 

P 



300(6, 25) ^ 1875 

+ (-1. 4) (2. 5) s+6. 25 s^-3. 5sf6. 25 



For s = and after some manipulation, it is found that 
, -11718. 75 +1875uu^-j3. 5(1) 

-L = ^ >! — j— 

(6. 25- i) ) + 12. 25u) 

and, in a similar manner for (5-10), but in a non- rationalized form 

-L = -62. 4 X 10^ 

" ( 51. 9 - (D^) - j 5 . 040) 



(5-12) 



(5-13) 



Qualitative plots of (5-12) and (5-13) appear in Figs, 5-6 and 5-7, re- 
spectively. The actual coordinate information for (5-12) and (5-13) appears. 



Smith, R. H. , An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, 
On the Limits of Manual Control, IEEE Trans, on Hum, Factors in 
Electronics , Sept. 1963. ) 
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Fig. 5-6 : Open Loop Unstable Plant (C 



I 

I 



I 




I 



I(cu) 



3 




129 



Fig. 5-7 : Open Loop Unstable Plant ((“ = 
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for selected frequencies, in Tables 5-4 and 5-5. A glance at the row for 
-L^ in these tables explains why Figs. 5-6 and 5-7 are presented qualitatively, 
i, e. , the scale is such that the universal overlays, given with this thesis for 
the D regions, cannot be used. 



TABLE 5-4 



uu 


. 5 


1 


2 


3 


4 


5 


-L^(juu) 


-288-j83. 6 


-246-J164 


-78-J243 


43. 8-jl67 


62. 2-J89. 2 


54. 8-j51. 2 



uu 


10 


15 


20 


30 


-L^(jLU) 


17. 6-j6, 6 


8. l-jZ.O 


4.61-j. 8 


2.1-j. 2 



TABLE 5-5 



UU 


1 


10 


20 




-1214 


617 


165 




-j647 


-j48 



In the face of an open-loop-unstable plant the plots must be interpreted 
with even greater care than before. The successful application of the variable 
parameter mapping technique to systems having open-loop instabilities may be 
of direct use in extending this study to the monitoring situation where the 
human operator is called upon to control a destabilized plant (or plant plus 
compensator). In the normal Nyquist plot, if a system is open-loop unstable, 
e. g. , has two right-half-plane poles, then for the system to be closed-loop 
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stable the net number of encirclements of the critical point must be exactly 
equal to the number of those right-half-plane poles and the encirclements 
must be in a CCW sense. This results directly from the number of encircle- 
ments yielding the difference between the number of zeros and poles of the 
do sed- loop- sy stem-func tion denominator (characteristic equation when set 
equal to zero). Here, zeros (CW rotations) are regarded + and poles (CCW 
rotations) are regarded -. For no zeros to be present in this example, the 
net number of encirclements N = Z - P must be N = 0 - 2 = -2. That is, for 
stability, the plot of the so-called "difference vector" as s traverses the 
entire right-half plane must undergo two complete counterclockwise revolutions 
about its origin (or, equally well, one CCW revolution for the positive fre- 
quency branch, since there is symmetry about the real axis). If N = 0, there 
are also two zeros in the right-half plane and the system is closed-loop 
unstable. What does a CCW encirclement appear as in the dual Nyquist plots? 

The single parameter variation as in Chapter Four and Figs. 4-8 through 4-12 
may be viewed to gain some insight into this question. Specifically, suppose 
the situation were as depicted in Fig. 5-8, where there is shown a hypothetical 

plot of some -L . Clearly, if iJU > oi , there is a CCW encirclement, \^ile 
^ n H ^1 

if uu - < 0) , there is no encirclement. If attention is now turned to Smith's 

fl VI 

results, it can be seen from the data points of both unstable plants that there 
are no encirclements, and since Smith used an unstable second order plant 
(two right-half-plane poles), it appears that this is an indication of two zeros 
in the right-half-plane and closed-loop instability. In helping to interpret the 
somewhat enormous plot of Fig, 5-7 (C = -0. 7), the following quote from (1) 
appears pertinent. "Controlled element gain had no apparent effect on the 
operator's ability to stabilize the system. Changes in K (K in this thesis) 
from 300 to 900 were barely noticeable by the operator at = 2. 5 rad/ sec. " 

For the C = -0. 35 case, the following comment appears. "The controlled 
element gain, K, (K in this thesis) was originally set at 200. The final 
value was 1200. Changes in gain produced no apparent change in the operator's 
technique. The only noticeable effect produced by gain changes was to decrease 
the maximum control deflection. " Since the Smith study addressed itself to 
the control problem presenting extreme difficulty to the human operator, i. e. , 

^^^Smith, R.H., An Experimental Determination of the Limits of Pilot Controlla- 
bility for Unstable, Oscillatory Second Order Systems , Thesis for degree of Aero- 
nautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper, On the 
Limits of Manual Control, IEEE Trans, on Hum. Factors in Electronics 
Sept. 1963. ) 
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Fig. 5-8 : Interpreting Open-Loop Instability Effects 
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controlling an unstable-open-loop system, one further quote from (1) is appro- 
priate. "Observations and operator comments during the conduction of the 
experiments for Q = -0. 35 indicate the appearance of a very important pheno- 
menon. In watching the operator at his control task, it seemed as though his 
reactions were intuitive. The changes in error were so fast and the error rates 
were so great that satisfactory control looked impossible; yet, somehow, suffi- 
cient stability was generated by the operator. After becoming accustomed to 
the control task the operator often remarked that he was simply outguessing 
the system. He made the statement that he was attempting to get his stick 
motion in phase with the error and trying to keep it there, i. e. , system 
stability was his foremost requirement -- keeping error minimum was not 
consistent with maintaining stability. " 

The task of interpreting Figs. 4-6 and 5-7 (actually the data of Tables 
5-4 and 5-5 in conjunction with Fig. 4-22), must now be considered; it is 
obvious that to stabilize the systems all the human gain must have been quite 
low. Unfortunately, the Smith study did not provide information relating to 

since Smith's stability criterion was taken to be positive phase margin for a 
corresponding amplitude ratio in terms of the calculation technique is not 
clear. The upshot, however, as recognized also by Smith, is that the general- 
ized transfer function for the human operator (2-1) is a model incapable of 
stabilizing the systems under consideration. Yet the human operators in the 
experiment were able to stabilize the system (5-11) for the conditions (5-9) and 
(5-10), though just barely. That is, these results are supposed to reflect the 
"limits" of manual control. The plant characteristics (5-9) and (5-10) define 
the boundary of marginality, i. e, that point at which the human operator can 
no longer control the system, for a plant having form (5-11). Greater confi- 
dence accrues to the predictive capability of this theoretical approach if it too 
displays such marginal conditions. In order to account for the experimental 
results. Smith hypothesized a second lead term for the human operator model 
so that stability results from the theoretical description, also. This tack is 

Smith, R. H. , An Experimental Determination of the Limits of Pilot Con- 
trollability for Unstable, Oscillatory Second Order Systems, Thesis for degree 
of Aeronautical Engineer, Univ. of Cincinnati, 1961 (also appeared as a paper. 

On the Limits of Manual Control, IEEE Trans, on Hum, Factors in Electronics, 
Sept. 1963. ) 
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not taken here. It is sufficient to demonstrate that the techniques of this 

thesis imply a marginal stability situation. Thus, for uu = 10 and K =01 (not 

h 

an infrequent value for human- operator gain) from Table 5-4, -L = 1. 8-j. 7, 

and it is conceivable that if K were slightly smaller than 0. 1, the -L locus 

^ n 

might intercept the small semicircle (which is not a function of frequency, 

Table 4-3) projecting below the real axis on all the 0 plots. Under such cir- 
cumstances, there is an enclosure and concomitant CCW encirclement indi- 
cating closed-loop-system stability. If it is felt that this discussion affords 
nothing more than an amateurish job of rationalization, it is also evident that 
the plots seem to indicate no clear cut conditions of one kind or the other. 

This substantiates that the mapping analysis adequately depicts the marginality 
of the situation. Application of the variable parameter mapping technique to 
the Smith data certainly demonstrates that the human operator is being pressed 
to the limit of his abilities. 

The possibility of applying the techniques of this thesis to the multi-aixis 
situation should not be discounted, i. e. , the study provides a theoretical method 
of approaching the multi-axis situation where a particular ’’nominal*' transfer 
function may not apply with invariance over all channels. The plots display the 
entire adaptive range of the human-operator transfer function and, since the 
bandwidth dealt with is so limited (0< JU< 20), the number of plots required to 
cover the necessary frequency range is pleasantly restricted. This may prove 
also, to be a more expeditious and fertile approach even for the single aixis 
problem, since the automatic part of the system and the human operator are 
displayed separately. This isolation permits observation of the control charac- 
teristics of the automatic part of the system from the same plots used with the 
human controller analysis. Indeed, the positioning of the locus relative 

to the (+l+j0) point gives precisely the information normally obtained in a stan- 
dard Nyquist analysis with respect to the (-1+jO) point. 

No experimental monitoring studies as such have been available to the 
author and none exist to his knowledge. Information vital to multi- axis moni- 
toring by human operators may devolve from the Elkind studies of human 
response to abrupt plant changes, i. e. , sudden changes in the gain and polarity 
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of a controlled element that was simply a gain; however, considerable experi- 
ment with more complex situations is needed to generate confidence in anv 
(1) (2) ^ 
inferences. The theory presented herein works remarkably well for the 

single- axis situation and an equivalent simplified monitoring topology with 
open-loop instability, i. e. , the theory predicts results consistent with the 
empirical evidence. Experimental results are needed to determine if the 
theory is totally applicable or, at the least, can continue to provide a valuable 
tool for design insight and "rule of thumb" criteria in multi-channel task 
situations. 



The difficulty with the experimental studies reviewed here is that all 
parallel compensators, G^, have been equivalently zero, Fig. 2-2. As indi- 
cated above there are no experimental results for monitoring. The question 
is: can theory predict that the human operator will have a less difficult task 

in a monitoring situation and how can this be determined? It would be instruc- 
tive to observe an controlling a with G^ = 0, then compare the human 
effort there with that required when an is paralleled by a G^ designed for 
good automatic control of the same G^. Such information is not available at 
present.* The problem is compounded by the valid question as to what role the 
human would play when the system is performing properly, i. e. , presumably 
the human enters only upon some catastrophic failure of one kind or other. 
Obviously, considerably more research is required in this area. 



5. 5 Summary 

The variable parameter mapping of the three time constants (which 
man can adapt to a given task) in the human-operator description, is used in 
conjunction with a slightly modified Nyquist plot of the plant characteristics to 



Young, Green, Elkind and Kelly, Adaptive Dynamic Response Characteris- 
tics of the Human Operator in Simple Manual Control, IEEE Trans, on Hum* 
Factors in Electronics , Sept. 1964. 

^^^Dander V. A. , Predicting Pilot Ratings of Multi-Axis Control Tasks from 
Single-Axis Data, IEEE Trans, on Hum. Factors in Electronics, Sept. 1963 . 
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demonstrate that stability of the closed-loop system (with a man in the loop) 
can be predicted. Experimental verification of the theory is obtained by 
analyzing certain systems which have been studied empirically in the past and 
for which results have been published in the literature. To this end, trans- 
parent overlays are provided which facilitate stability determination. The 
theoretical predictions parallel the experimental results remarkably well; 
one of the sets of data relates to a system which is open-loop unstable. The 
correlation possible between boundary lines of the mapping and the physical 
capacity of a human operator is cited and, in the cases of stable control, the 
margin of human capacity remaining (displayed on the same diagram) is 
indicated. 
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CHAPTER SIX 
PERFORMANCE 



6. 1 Introduction 

As in the discipline of inanimate systems engineering, several performance 
measures have been used for man-machine systems to provide a so-called figure 
of merit or index of quality and, as in the case of the all-hardware system, defi- 
ciencies have been found to exist with all these scalar indices. The difficulty in 
attempting to quantify a complex system's performance by virtue of a single num- 
ber is that satisfying the requirement on this unique scalar may still result in a 
system with other unsatisfactory characteristics. Of course, the designer of per- 
formance measures answers that all pertinent characteristics are duly weighted 
and included in the scalar index thus causing it to reflect those properties deemed 
important. Still, the final selection of those properties regarded as necessary for 
an acceptable system and assignment of relative weights are only two of many rea- 
sons for viewing single performance measures with some caution. Generally, the 
selection of performance index has often been predicated more on ease of calcula- 
tion and other conveniences than upon how accurately the truly efficacious char- 
acteristics of a system are reflected. In the psychological literature, time on 
target, an easily obtained number has been used extensively although studies have 
shown its deficiencies; it has even been demonstrated that there is considerable 
danger in inferring that system-measure changes imply system-component 
changes. ^ Other indices used by psychologists studying manual control have 

been average error, average-absolute error and root-mean- square error; attempts 



Bahrick, H. P. , P. M. Fitts, and G. E. Briggs, Learning Curves-Facts or 
Artifacts, Psych. Bull. , 54, 256-268, 1957. 

( 2 ) 



Taylor, F. V., and H. P. Birmingham, That Confounded System Performance 
Measure-A Demonstration, Psychological Review, 66, 178-182, 1959. 
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have been made to measure task difficulty but independent critics have dismissed 
them as measuring little about the specific requirements of a task. The engi- 
neering approach to the performance measure quandary has been to employ the 
seemingly dependable (and quite convenient) mean-square error, as well as cer- 
tain other control theory measures. The period of time over which the measures 
themselves are calculated serves to influence the result and thereby confuse the 
problem still further. Historically, pilot opinion has been used as an acceptance 
technique in the case of aircraft flight control systems. In this regard, the fol- 
lowing quote is considered appropriate. "It should be clear at this point that a 

single measure for a complex tracking situation cannot give information concerning 

_ 

both system performance and operator performance" ( underlining theirs). 

In this chapter, the various performance measures are discussed and a 
brief survey is given of some techniques available for "optimizing" a system so as 
to best satisfy these indices. The formulation of a "vector" performance index is 
considered as a means of circumventing certain of the disadvantages of scalar 
indices. The concept of a "best-overall" compensator (fixed compensation), for a 
plant with variable dynamics, is introduced and a technique for specifying the para- 
meters of this compensator is offered. This technique depends upon the formation 
of an "artificial" function which, when minimized to zero, identically satisfies a 
set of performance requirements. A "direct search" procedure, used in conjunc- 
tion with this "artificial" function, which is sequentially modified, serves to isolate 
satisfactory values for the compensator time constants. In this manner, the space 
of acceptable compensator time constants is defined and correlated with the space 
of required performance. The entire procedure is especially amenable to the 
electronic computer. An example of the search technique is given at the end of 
the chapter. 



Taylor, F. V. , Human Engineering and Psychology, in "Psychology: A Study 

of a Science", Koch, S. , (Ed. ), Vol. 5, McGraw-Hill, 1963. 

Summers, L. G., and K. Ziedman, A Study of Manual Control Methodology 
with Annotated Bibliography, NASA Contractor Rpt. CR-125, Nov., 1964. 
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6. 2 Varying Parameter Plants and "Best” Performance 

The question frequently arises how ’’best” to use the human operator in 
any given system; it is assumed that the stability criterion is satisfied. Given 
the requirement for a human operator one may ask, alternatively, how to design 
a satisfactory compensator for use in a system which includes this human con- 
troller? The design task may be complicated further by specifying a plant which 

itself contains some varying parameter. Generally, the controls engineer has 

✓ 

tended to avoid computer design for the more familiar Bode diagram, frequency 
plot or root locus technique. An approach to the performance problem, is pursued 
in this thesis which, it is hoped, can be suggestive of a standard technique, or 
algorithm, for handling far more severely complex situations than can be treated 
here. The method is particularly applicable to systems containing human con- 
trollers since the human adapts his characteristics to some unknown internally- 
generated criterion (which is frequently taken to be some function, or functional, 
of the error). The situation is quite fluid in the face of systems with varying- 
parameter plants since the error itself then changes due to this parameter varia- 
tion thereby causing further adaptation on the part of the human. It is felt intuitively 
that some design compromise on an adequate compensator should be possible so as 
to provide the ’’best” performance over some range of parameter variation. The 
question immediately arises regarding a proper criterion (performance index); 
this is a subject around which much controversy still rages within the engineering 
profession. Frequently, the physics of the situation suggests a performance 

index which leads to a mathematically intractable problem. The designer then 
must compromise on a choice of performance index which permits a meaningful 
solution, i. e. , does not preclude analysis, yet still provides a useful quality mea- 



Newton, G. C, , L. A. Gould, and J. F. Kaiser, "Analytical Design of Linear 
Feedback Controls", John Wiley ^ Sons, New York, 1957. 

Gibson, J.E., "Nonlinear Automatic Control", McGraw-Hill, 1963. 
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sure for the system. Usually, the indices yielding the greatest analytical ease in 
their calculation are the mean-square error and the integral- square error. For 
inputs having Gaussian characteristics and certain other properties which are 
frequently encountered (or at least frequently assumed) it can be shown that sev- 
eral error criteria offered as more appropriate lead to the same results as the 
mean square error criterion. ^ ^ Of course, the advent of high speed computation 
and concomitant expeditious numerical integration offers great latitude in the 
choice of a performance index. 

Some work has been done indicating that minimization of the mean- square 

error is a fair approximation to what the human controller attempts, i. e. , it 

(2) 

appears the ultimate effects are similar. Generally, the mean-square error 

is used as a criterion for stochastic inputs; indeed, this is the basis upon which 

( 3 ) 

the human-operator transfer function has been generated. For transient inputs, 
the analyst frequently uses the integral square error criterion; of course, the 
ISE is used only in the design of systems having zero final error since continual 
finite error leads to the non-convergence of the ISE index. With either measure, 
however, stability checks must be made at discreet (also discrete) moments to 
insure that the values of the adjustable parameters which minimize the performance 
index correspond to a stable solution. This is elaborated upon below. Thus, one 
finds that the MSE is used as an index of, loosely speaking, steady state perfor- 
mance while the ISE is used to measure transient performance characteristics. 
Certainly, the selection of a performance index is at the discretion of the designer 
(or the customer) and modern digital computation offers the possibility of numer- 
ically computing heretofore mathematically unattractive but physically useful 
response criteria. 



Sherman, S. , Non-Mean-Square Error Criteria, IRE Trans, on Information 
Theory, Sept. , 1958. 

/ 

' Roig, R. W. , A Comparison Between Human Operator and Optimum Linear 
Controller RMS-Error Performance, IRE Transactions on Human Factors in 
Electronics, Mar. , 196Z. 



(3) 



Graham, D. , and D. McRuer, ’’Analysis of Nonlinear Control Systems'*, John 
Wiley & Sons, 1961. 
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The remainder of this section is confined to the exposition of a possible 
approach to solving such a "best" performance problem. The attack is along 
conventional lines while that of the next section is considerably more original 
and exotic. The system studied here is somewhat simplified but it should be 
recognized that the simplification incurs no loss of generality. 

The system under investigation appears in Fig. 6-1. The configuration 

within the dotted lines represents a "choice" situation, i. e. , the human operator 

either chooses to operate upon the error with his dynamics or he does not. Again, 

such a situation where the human operator decides whether or not to enter the 

circuit may be considered as an alternate form of the monitoring situation. The 

representation here does not place the human operator in parallel with the original 

compensation but in cascade. It is recognized that the design of for good 

automatic control may be distinct from a good design with a human controller in 

the system; therefore, the human is cascaded always with some such that 

G = G^ G provides good compensation for the human controller-plant com- 
C c c 

bination. It will be observed that no loss of generality occurs if the situation is 
depicted as in Fig. 6-2 since the plant-human path is under study. The human- 
operator transfer function is considered in the simplified form 



H 

o 



-. Is 

e 

s + h 



( 6 - 1 ) 



where h can be varied adaptively over some range. The unknown compensation 
to be determined is cast in the form of a straight gain, a , and the plant has the 
form 

G = (6-2) 

^ s + 3 s + 1 

where 3 varies over some range. It should be noted as the study continues that 
any class of more complicated frequency- sensitive compensation, a plant with 
several varying parameters or a mathematical model for the human operator con- 
taining more adaptive variables, only makes the calculations more severe but does 
not change the general technique. Thus, Fig. 6-2 takes on the form of Fig. 6-3. 
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Fig. 6-1 




Fig. 6-2 




Fig. tJ-3 : Human Control of Va rying-Pa ramcto r Plant 
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Since the study is restricted to deterministic inputs (steps, ramps, etc. ) the 
performance index is chosen to be the ISE, 



ISE = = j e^(t) dt 

o 



( 6 - 3 ) 



where e(t) is the error signal. 

A design is desired, i. e. , a value for OC , such that the ISE is minimized. 
The standard analytic technique for solving this problem uses a rational transform 
approximation for e 



-Ts 

e 



1 - 



Ts 

“Z“ 




(6-4) 



❖ 

This permits evaluation of the ISE as a function of the free parameters. One may 
then attempt to minimize the integral by taking derivatives of ISE with respect to 
each of the free parameters setting them equal to zero and solving for the critical 
value of each parameter. Thus 



bl bl bl 



(6-5) 



The step above is often more easily described than accomplished, but the 
main difficulty is not merely that of explicit solution. In the case with a human in 
the system, one must assume there is an attempt to adapt operating procedure so 
as to minimize some function of the system error, irrespective of the value of a 
and 3. The 3 variation is not controlled and may represent many possible 



The evaluation is easily done by means of tables. 



( 1 ) 



U) 



Newton, G. C. , L. A. Gould and J.F. Kaiser, "Analytical Design of Linear 
Feedback Controls". John Wiley & Sons, New York, 1957. 
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phenomena; e. g. , change of response characteristics with altitude by a high- 
performance aircraft, drift in plant parameters, etc. The problem may be ap- 
proached as one wherein Nature is in direct opposition to satisfactory perform- 
ance and a design is sought such that performance is at least as good or better 
than some value of the index of merit, in this case, ISE. For tutorial purposes, 
consider CX and each having such range that they are adequately described by 
a grid of five (5) values, each grid covering the respective ranges of a and 
Thus, for each pair (oc,0) there is obtained that value of h^ which minimizes 
the general performance index S; this index is a function of the three variables 
a,3,h and is taken in ISE form. The minimization of 

00 

s(a,3,h) = J e^dt (6-6) 

o 



is then given notationally by 

CO 

1. min Ja^dt hVS) (6-7) 

h o 



i. e. , for each over the chosen grid 



(a. 3) = (0.1) h° 



01 



= (0,2) 



02 



25 h" 



= (i, j) 






This step represents the effect of assuming that the human operator always does 
the ’’best" he can, where "best" is the optimization of some error criterion (here 
the identical performance index, ISE, is retained throughout as the criterion but 
this is not mandatory, i. e. , distinct performance indices can be minimized for 
h, OC cind 3, respectively). 

In the next step, for each value CC can take on there is found that value of 
3 which maximizes the ISE. Here one observes the result of assuming the 
"cussedness" of the plant; this is, at times, called "worst case design". The 
disadvantage of this approach lies obviously with the extreme conservatism of the 
final solution. Thus, ^ 

2. max 

0 o 



e^dt 



3°(a) 



(6-8) 
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i. e. , for a = 1, run through g = 1 5 to get an optimal 3° 1 

* 1 5 3°'s 

a = 5, run through 3 = 1, .... 5 to get an optimal 3° I 

00 ^ I 

where the "optimal" 3 maximizes J e^ dt. 

o 

Finally, one looks at the five values of (6-6) for a = 1 5 and picks the 

minimum of these, thus yielding the "best" a. Hence, 



3 . 



00 



min J e dt 
a o 




(6-9) 



The above procedure insures that performance 

00 

S(a^, 3 , h^) < value J e^ dt = S(a^, 3^, h^) (6-10) 

o 

where 3^, it will be recalled, refers to the "worst" situation. For those familiar 
with Game Theory, the use of a game theoretic approach in the technique is notice- 
able, i. e. , a minimax solution is being obtained. Clearly, even in this extremely 
simplified example many calculations and evaluations are required and frequent 
comparisons between them must be made. The technique is especially suited for 
digital or other machine computation. 

Two additional points with regard to this approach should be mentioned. 

The first step (the adaptive action) yields an h^(CX, 3) from which one can proceed 
to evaluate S(a, 3) as indicated above, but care must be taken to insure that 
(a,3)^R^; that is, the pair (a, 3) must correspond to a region of stability. 
Otherwise, S(a, 3)“"® if (0',3)^Rg. The parameter space plot of Fig. 6-4 
should clarify this. If the region of stability has the form sketched in a 3 space 
and the range on a and 3 is indicated by lines (l)-(2) and (3)-(4), respectively, 

then it is clear that during the step max S(a, 3) for the value a = a the criterion 

3 

is maximized by picking 3=3^, since this combination corresponds to an unstable 
system. Either a manual check must be conducted to ensure that the optimal 
values determined correspond to a stable condition or the computer program must 
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Fig. 6-4 : Solutions Outside the Region of Stability 
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include a check for stability during the computational cycle. 

If the transform of e(t) is a ratio of polynomials then the evaluation of 
ISE is straightforward, by virtue of previously prepared tables, and only the 
minimization is carried out by computer. If e(t) transforms into an irrational 
form (as it does with a pure time delay) then evaluating the ISE is not so simple. 
The integral may be truncated, however, and evaluation carried out numerically. 
To follow such a procedure, however, reduces any feeling of commitment to ISE, 
or other mathematically attractive figures of merit, since their calculation would 
require as much work as that involved in other more physically appealing per- 
formance indices. Also, a coordinated computer program may be possible such 
that, while the performance index itself is being calculated, the information 
generated in that process is used to simultaneously minimize the index and pro- 
vide a stability check. 

6, 3 The Concept of a ^'Vector'* Performance Index 

A technique is now considered for determining the ’’best-overall" com- 
pensator for any given system with a human controller. The approach in this 
section may be used for a plant with variable parameters; the only restriction is 
the usual aggravation imposed by the "curse of dimensionality", i. e. , the larger 
the number of varying parameters the more iterations required, The solution 

to be presented depends intimately upon an initial stability study conducted in 

❖ 

accordance with the procedure previously outlined in this thesis. Thus, for the 
situation of Fig, 6-5, the result desired is values for K^, Cl and 0 which pro- 
vide good, overall performance, defined by specifications on the range of the 



^ Bellman, R. , '’Adaptive Control Processes-A Guided Tour", Princeton Univ. 
Press, Princeton, New Jersey, 1961. 

The author has been unable to find in the literature information relating to the 
step response of a human- operator-frequency- sensitive-plant combination, 
although there is considerable published on such responses for plants represented 
by pure gains. Thus, the study represents a theory which is yet to be checked 
empirically, although all previous experiments would seem to justify this 
theoretical extension. 
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triple (t^, o^, tg), in the face of a plant which has a varying C . The triple 

(t^, o^, t^) represents the values of rise time, overshoot and settling time, 

respectively, which follow in accordance with some arbitrary definition. Of 

course, the triple may be regarded as a vector quantity in three space . The 

stability analysis outlined in Chapters Four and Five is first accomplished; for 

the class of compensator being considered, one zero and one pole, it can be 

determined quickly what approximate ranges of values on Cl and 3 provide good 

stability. The value of is not difficult to determine since a nominal value of 

may be used, thus indicating quickly what the overall loop gain should be, 

from this a first approximation to may be calculated. Refinement of this 

value is made along with the search for final values of a and 3. Gain adaptation 

is the easiest act for the human operator to perform; indeed, his task is regarded 

as simplest when he acts merely as an amplifier. ^ This aspect of the variation 

therefore can be discounted in the remaining analysis. The competent designer 

desires a final empirical test anyway, if only by simulator, overall operation can 

be double checked at that time. In Section 6, 5 an example is given which includes 

consideration of K , 
c 

In addition to satisfying the general stability requirement the designer does, 
of course, attempt to use the freedom available in selecting Cl and 3 so that the 
plant's Nyquist curve passes through that part of the Q region corresponding to 
reduced task requirements on the human operator, (Chapter Five). The entire 
analysis so far has been conducted for a nominal value in the range of variable C» 
Generally, there are some ranges on a and 3 which seem to meet the conditions, 
or it is apparent they are impossible to meet under the constraints imposed, or a 



Quality measures other than the three suggested are, of course, possible; the 
technique permits considerable latitude in employing multidimensional per- 
formance measures. 

Birmingham, H.P. and F. V. Taylor, A Human Engineering Approach to the 
Design of Man-Operated Continuous Control Systems, NRL. Rep. 5333, 7 Apr. , 
1954 (also published as A Design Philosophy for Man-Machine Control Systems, 
Proc. IRE, Dec., 1954.) 
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different class of compensator may be required. With the ranges of a and 3 
that provide acceptable stability available, nominal values for the human controller 
can be read from the plot using the time-invariance asssumption. These values 
must also retain the high-frequency" end of the compensated Nyquist contour in a 
location^ satisfying the stability criterion, as reflected by the corresponding 0 
region. That is, the triple T^) must have value such that it lies to the 

left of the plant* s Nyquist curve and still be within the range of human adaptivity. 
Simultaneously, T^ and T^^ should be as small as possible; particularly 

since this parameter is almost a direct indication of "comfortable" adaptation. 
All these conditions must also be consistent with physiological realities. 

With the analysis above conducted for a nominal C » it is then repeated for 
a grid of C values covering the range of Q. The intersection’ ‘ of the various 
ranges on a and 0 yields those ranges of (X and 0 which satisfy stability over 
the expected range of variation in Q, These final ranges of acceptable Cl and 0 
must now provide the specific values of (X and 0 which yield the best performance. 
Again, a vector (rather than a scalar) definition of performance is being used. The 
requirement is to meet specifications on rise time, t^, overshoot, o^, and settling 
time, t^ . There is no problem in determining what values of t^, o^ and t^ 
result from any given set (a,0); the question is how to go about selecting sub- 
sequent trials of cx and 0 so that each iteration is at least as good as the last. 
What is sought is some "scientific method" for approaching the sub- space of 
acceptable performance within the space of acceptable stability. The following 
approach provides a routine technique, or algorithm which lends itself to com- 
puter solution. With the resolution of ranges on CX and 0 which provide satis- 
factory stability characteristics, there is conducted a search for the sub-range 



Inherent in this statement is the asumption that this is what the human does; 
i. e. , it is hoped this adequately reflects his unknown performance criterion. 

That range over cx and 0, respectively, whose end points are within all the 
other ranges of cx and 0 determined for various Q is called the intersection 
(greatest common range of cx and 0 over C). It should be noted that this 
step considerably reduces the experimental region remaining. 
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which yields good performance, Fig. 6-6(a) and (b). Quite possibly none may 
exist; this situation is discussed later. A grid search of Cl3 space could be 
performed to find acceptable (t^, o^, t^), but that is a strictly trial and error 
method. If the (a, 3) range is small, straightforward grid search could be the 
most efficient technique; if there are more than two parameters available for 
manipulation, however, the algorithm is probably more efficient. 

The function 

t = (6-11) 

is formed where 0<t <T,0<o <O,0<t <T ,i. e.,T,0, and T are 

r— r s— s s—s rs s 

the maximum acceptable bounds on t^, o^ and t^, respectively. It should be 

noted that the function f has the interesting property that when it is zero the 

specifications are met identically at their upper bounds, i. e. , cx and 3 have 

values such that the surface boundary of the acceptable performance region in 

(tr» o^, t^) space has been attained, Fig. 6-6(b). The function is an artifice 

created specifically to facilitate finding some point on this surface and, in this 

instance, is still a function of two independent variables Cl and 3, although its 

ultimate evaluation requires three numbers t , o ,t . Thus, it is desired to 

r s s 

minimize f over the space of (l3 which yields stability, recognizing that due to 

the sum of squares form of f it cannot be negative, while to meet specifications 

just barely, min f must equal zero. The form of the function therefore precludes 

the occurrence of saddle points when the min of zero is attained. The search 

techniques for "optimum" a and 3 must incorporate certain characteristics to 

keep the designer adequately informed. 

First, it is desired to check each of the terms t-T ,o-0,t-T 

r r s s s s 

algebraically; in the event that they are negative the program should STOP 
and PRINT OUT such information since this is a reflection of an adequate design 
whose response falls inside the acceptable performance region. Second, when 
the search pattern finds min f a PRINT OUT should again be required, since it 
must then be determined by the designer whether the results are sufficiently close 
to zero to warrant an interior search of the acceptable performance region. If 
they are not, this is an indication that the requirements on t^, o^, t^ cannot be 
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met and a different class of compensatoi' may have to be used or other compromise 
fashioned. The utility of such information should not be regarded lightly since it 
represents an analytical demonstration of inconsistent or impossible specifications. 
Too frequently there have been attempts at building systems to meet such impos* 
sible requirements and only successive futile results have led to their abandon- 
ment; any technique which uncovers such sets of impractical demands serves the 
important purpose of precluding the expenditure of funds and effort on an un- 
realizable quest. 

Having located the surface of the acceptable region of performance, one 
proceeds in the "same direction" until the opposite side of the region is found in 
which case the entire line between (and the associated values of a and 0) must 
be contained inside. Many times the "opposite side" is known due to constraints 
on the parameters, e. g. , CC and 3 must be 0. In this manner, the region can 
be crudely, or completely, defined by successive computation along lines per- 
pendicular to the one originally found. 

6. 4 The Gradient Method 

Frequently, a steepest ascent (descent) technique is employed in such pro- 
blems. This type of approach is described briefly, although it will not be used; 
the description allows comparison with the procedure to be given in detail some- 
what later. From (6-11) for an arbitrary (a^,3j) the function f (a^,P^) is deter- 
mined, i. e. , for the values (a^, 3^) the system response to the input function of 
interest (here a unit step) is calculated and °g’ found thereby permitting 
f to be enumerated. For some (6cx, A3) the functions f (cx^+Aot, 3j^) and 
f (a^,3j+A3) are determined similarly. The function below is formed 

f(a^ + Aa, 3j) - f(Oj,3j) 

Ka 



^ Wilde, D. J., "Optimum Seeking Methods", Prentice-Hall, New Jersey, 1964. 
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and also, 



f(a^,3^+A3)-f(a^,3^) 

A3 



In this way, estimates are obtained on 
But the gradient of f is 






_ , 



■5a 



5f-r- 

Wh 



5f 



respectively. 



where and are unit vectors in the positive oc and B directions, respectively. 
The gradient of f, of course, points in the direction of greatest increase of the 
function f. Since it is desired to minimize f the program proceeds in the opposite 
direction. By surrounding the point by a pattern of, say, two points 

oriented as in Fig. 6-7, the gradient at each of these two points may then be 
determined and iteration continue in an appropriate direction. Thus, the direction 
of incremental change in a and 3 yielding the greatest decrease in f (or, increase 
in minus f) would be followed to obtain the next trial point. The search procedure 
is then repeated as one follows one's nose to a minimum. When the optimal point 
(be it min or max) is apparently reached it is good practice to conduct a search in 
the local area surrounding this isolated point to ensure that it is not an apparent 
optimum. For this purpose, non-linear exploration of the adjacent region is in 
order; Box and Wilson have discussed the problems inherent in such studies. 



6. 5 Direct Search and Pattern Strategy. 

The upshot of the previous discussion is that any of several different type 

search techniques can be used to determine the "best" Ol and B. The approach 

used in this section, and which seems to offer the greatest advantage for the sys- 

(Z) 

terns under consideration, follows that of Hooke and Jeeves. This technique of 



Box, G. E. P. and K. B. Wilson, The Experimental Attainment of Optimal 
Conditions, J. Roy. Stat. Soc. , B13, 1(1951). 

Hooke, R. , and T. A. Jeeves, Direct Search Solution of Numerical and 
Statistical Problems, J. Assoc. Comp. Mach. ; 8, Z, Apr. , 1961. 
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Fig. 6-7 ; Perturbation Pattern for Gradient Determination 
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direct search with pattern strategy is especially suitable to minimization of a sum 
of squares, but its overriding virtues are its straightforward simplicity and com- 
patibility with the electronic computer. 

The starting techniques in pattern search are quite similar to those of the 
gradient or steepest ascent method, however, it appears at first glance that pat- 
tern search neglects some information that is readily determined and used in the 
gradient method. The effects of this neglect are more apparent than real. By 
"direct search" is meant a step-by-step examination of trial solutions and a 
method of comparing each to determine which is the "best" in some sense. Addi- 
tionally, there is some method offered for choosing the next trial point. This last 
aspect forms the strategy part of the search procedure; a "pattern*’ approach is 
used and from this is derived the strategy’s name. 

The example to be used is shown in Fig, 6-8, It is assumed that the 
ranges of a and & have been discovered which yield good stability; specific 
values within that range are now desired which meet the step-response require- 
ments on rise time, overshoot and settling time. The compensator gain K^, is 
included as a third parameter whose value must be determined. The ranges K^, 
a and 3 form' a space of points representing possible solutions. Again the func- 
tion, f, is formed 

f = (t -T )^+(o -O )^+(t -T (6-11) 

r r s s s s 

For economy of notation, let the compensator gain, denoted, K, then the 

point (K 2 ,a^,P 2 ) is said to be "better" than (Kpa^,g^)if 

f(K2,a^,32) < (6-12) 

The values K^, a^, B ^ are determined by adding some increment 6 to each value 
Kp 0.^, and 6 ^ ; the increment size used for each parameter need not be the same. 
This portion of the search procedure is called the exploratory phase. It is under- 
stood that by f (K.,a.,3.) is meant the value of (6-11) found from the substitution 



The inequality is given for a minimization problem. 
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Fig. 6-8 : Compensation of a Va rying-Parameter System 




C = O.OI 



Fig. 6-9 : Sketch of ” Plots 
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of the t^, and associated with the system step- response when the com- 
pensator has gain and time constants K. , a. and S., respectively. The initial 
base point is chosen arbitrarily or the investigator may exercise 

his hunches . There may be none, one, or more points such that a min f of 
zero is found. The ramifications of this have been discussed above. Under the 
assumption (6-12) holds, there seems little reason in the face of such success, 
to alter the direction of search"; therefore, exploration is continued by per- 
turbing the independent variables with increments whose proportions are based 
on the number of successes achieved along each dimension in the recent past. 

This strategy forms the so-called "pattern" aspect of the search procedure; of 
course, great strides are made with pattern moves in relation to progress 
resulting from exploratory moves. An example is offered below. In those situa- 
tions where a min f = 0 is obtained, the direct search method continues to offer a 
consistent philosophy. The designer unquestionably desires additional information 
with regard to the size of the acceptable performance region and does not wish to 
stop at the barely acceptable surface bound, e. g. , other satisfactory triples 
(K^,a^, 3.) may offer simpler design realizability. Direct search can forge across 
this region until it pierces the opposite surface; in this manner the performance 
region may be completely explored. 

Termination procedures are required when a pattern move failure is fol- 
lowed by consistent failure from all exploratory moves. A retrenching is occa- 
sioned by reducing the nominal step size of the increments and new explorations 
are conducted; this is continued until a new pattern is constructed or the step 
size falls below some pre-selected margin. The lower bound on this margin is 
usually imposed by the demands of the particular computation interval peculiar to 
the computer used. Termination does not imply an optimal solution so the wary 
investigator tests the solution generated in radial directions about the final base 
point distinct from the coordinate axes since pattern search restricts itself to 
motion along these axes. 

6.6 "Best-Overall" Compensation of a Varying Parameter System - An Example 

The system for which a best-overall compensator is to be designed is 
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shown in Fig. 6-8. It should be noted that the human-operator transfer function 
has been omitted and a simplified plant assumed; the simplification is imposed 
for clarity only. In this manner, the techniques to be demonstrated are not 
obscured by technical complexities; no loss of generality occurs. The plant*s 
damping ratio, C, is assumed to vary over the range 0. 01 < C £ 1. 0; the rise 
time, t^, is defined as the time required for the response to go from 10 to 90 
percent of its final value; the overshoot, o^, is the maximum departure from 
final value expressed as a percentage; the settling time, t^, is defined as that 
time after which the response is confined to a tolerance band of ± 5 percent 
about its final value. In order to demonstrate the search technique for more than 
the two parameters a and 3, the value of gain, K^, is included as an unknown in 
the search. 

Sketches of the negative Nyquist contours for the plant, for values of 

C = 0. 01, 0. 1 and 1. 0 are shown in Fig. 6-9. The limiting ranges on a and 3 
are taken to be 0 < a < 5. 3, 0^8^ 25. The variable-parameter mappings of 
Chapter Four can be used to aid in the design, since they are not required to 
represent the capabilities of a human operator, i. e. , has been assumed equal 

to unity in this example in order to illustrate the search method. Thus, the 
degenerate mapping shown in Fig. C-2, Appendix C, (for uu = 1.0) can be used to 
represent the available region of variation for the compensator's time constants, 
a and 3 . Similarly, the corresponding regions on the transparent overlays (in 
the pocket of the back cover) can be used to provide additional frequency effects 
other than for jj = 1.0. 

It is evident from the overlays, used in conjunction with the more accurate 
plots in Fig. 6-10 (or root locus considerations), that lag network configurations 
( a < 3 ) severely aggravate stability when Q is very small, i. e. , between 
0. 01 < Q < 0. 1. Appendix C contains further stability interpretation elaboration. 
From the overlays it is apparent that after introduction of a compensator, stability 
is possible for small Q (C»0. 01) only if the compensation is lead type. The over- 
lays also indicate that the furthest ’'to the left" the critical points governing sta- 
bility can be placed (for the curves representing 0. 01 C_£_ I- 0) is along the 

circular arc segment 3*-7\ since this arc is the image of the condition 
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Fig. (--10 ; Controlled Element Plots, for various 
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0 < a < 5. 3, [3=0. Also, in passing it may be noted that the condition 
0<a<5.3, 3=0. 7 is represented by arc T-S* (Fig. 4-18 indicates the definition 
of these points), i. e. , 3 would play the role of in the human-operator trans- 

fer function and the tiny hull formed by line 3' 1*, arc l*-5', line 5'7’ and arc 7*-3* 
would be of interest, (detail in Fig. 4-18). Since this region remains reasonably 
unchanged in the overlays for increasing frequency, it is possible to use the over- 
lay for 'S = 1, Z or 3, respectively, for an overall estimate of what is needed for 
a, 3 values. Furthermore, it is obvious that common values of a and 3 which 
provide stability over all C, can be found from those values of Ci and 3 yielding 
stability for the condition C = 0. 01, since it is this value which imposes the most 
severe stability requirement. Also, the plots suggest that increased gain might 
help the overall situation across the range of varying C, since the C = 1. 0 con- 
dition most assuredly benefits from increased gain, while the effects of increased 
gain upon the C = 0. 01 condition can hardly place that locus much closer to the 
real axis. Increased gain, however, may improve stability (for C = 0. 01) by 
shifting the frequency points further out along the locus and therefore further away 
from the region of criticality. Fig. 4-18 detail. Thus, values of ol ^ 5. 3 only 
place the critical points depressingly close to the plant’s Nyquist loci (all C)i while 
for values of 3 > 0. 7, a large portion of the circular segment for 0 < a < 5. 3 lies 
in a region of instability for much of the range of 0. 1. The largest overall 
degree of stability across the varying range of C appears to be given by 
3 < 0. 7, 0 < a < 5. 3. 

It is reasonable, therefore, to begin the search for the oc and 3 which 
meet specifications with starting values of oc = Z, 3 = 0. 5. The specifications are 
taken to be 

t = Z seconds = T 

r r 

o =15 percent = O 

s s 

t =4 seconds = T 

s s 



Since the problem includes a search for the proper gain setting, K^, this parameter 
is chosen to have a nominal value of unity, which represents the original conditions 
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for which this preliminary analysis has been conducted. It is noted that the plots 

suggest a final value of gain for which is greater than unity, I. e. , the Nyquist 

curves of Fig. 6-10 indicate that increased gain while not affecting the C = 0 . 01 

and C = 0 . 1 curves substantially, ought to improve the response for the C = 1. 0 

condition considerably. Thus, the starting values for the search are taken to 

be a = 2, 3 = 0. 5, K =1. 

c 

Also, 

f(K^,a,3) (6-13) 

where 

Tr = 2, 0^=15, T^=4 (6-14) 

The step sizes chosen for the intial exploratory move are as follows: 

For K^: take step 6 ^ = 0. 2 (6-15a) 

For a : take step 62 = 0.1 (6-15b) 

For 3 : take step 62 = 0.05 (6-15c) 

These exploratory increments have been chosen arbitrarily and, for explanatory 
purposes, they have been made distinct; obviously, they can also be set equal, 
if desired. The subscripts 1, 2 , 3 are assigned to represent step increments along 
the and 3 dimensional axes, respectively. Usually, the problem and the 

designer*s experience suggest basic step sizes for each independent variable. 
During the termination procedure, the nominal step sizes of the increments are 
reduced until no further reduction is possible due either to solution indifference 
or equipment limitations. 

Due to the pressure of time the search was conducted on analog- computer 
equipment rather than digital; the most satisfactory equipment for this type of 
problem is a hybrid computer containing both analog and digital capability, such 
equipment was not available. The computer was a Mid-Century Instrumatic 
(5068), a Mid-Century MC-300 six-channel recorder was used to obtain traces of 
system input signal, response, error, etc. The computation scheme is shown in 
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Fig. 6-11. If in Fig. 6-8, one lets 
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(6-16) 

(6-17) 

(6-18) 

(6-19) 

( 6 - 20 ) 

( 6 - 21 ) 



then Fig, 6-11 is a represenation of the above relations. 

It should be recognized that as a result of performing the search on analog 
equipment, the step-response characteristics (t^, o^, t^) of the system were 
determined visually by reading the pen-recorder traces for v(t). Once a near- 
optimal solution was obtained, it was not possible to reduce the increment sizes 
(6) any further since any differences in the pen- recordings, for smaller 6 values, 
were not distinguishable. 

Generally, the difficulties in meeting the specifications are worse at the 
C = 1. 0 end of plant variability, i. e. , worse for t^ and t^, therefore, the search 
is started at that end. The initial base point for the search is designated, b^, 

Fig. 6- 12, the horizontal bar signifies the vector (multi-dimensional) nature of 
this point, i. e. , 

b^ = (1,2,0. 5) 



( 6 - 22 ) 
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Fig, 6-11 : Analog Computer Setup for Direct Search Example 
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Fig. ()“12 : "Direct Search" with "Pattern Strategy" 

Only the first two iterations are shown 
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which is obtained from 



K^=l, a = 2 . e =0.5 (6-23) 

The step response of the system with a compensator which possesses such char- 
acteristics is 

(tr’Os'^s^ = (4. 8, 0,9) (6-24) 

Therefore, 

f(l, 2 , 0. 5) = [ 4 . 8-Z]^+ [0-ls]^+ [ 9 - 4 ]^ = 257.81 (6-25) 

This value for the artificial function, f, is now compared with the value obtained 
at the next observation. The next observation is performed at the point b^ , 
where 6^ is a vector with all components zero except the one representing the 
independent variable, K^. Thus, 

bj+'S’^ = (1.2, 2, 0.5) (6-26) 

since |"5jl = 6 = 0. 2, from (6- 15a). From these parameter values the response 

is found to be 

(t . o , t ) = (3. 5, 0. 6. 8) (6-27) 

r s s 



Therefore, 

f(l. 2, 2, 0. 5) = [ 3 . 5-2]^+ Lo-15]^ + L6. 8 - 4 ]^ = 235. 05 (6-28) 

In general, when there is little or no intuitive feel for a potential solution and one 
direction offers as much hope for success as another during the initial exploratory 
moves, the response for b^^ - 5^ is also obtained to determine if it yields better 
results than b^ + 6^ , or b^. In this case, since 

f(b^ +7^) < f(b^) (6-29) 

and, since it is recognized that a sluggish system (C = 1. 0) becomes even more 
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sluggish for reduced gain, the value of f produced by + 6^ is accepted as the 
smallest and b^ +6^ is designated the temporary base point, The subscripts 

indicate this is the first temporary head (in a vector sense) and only one variable 
(K^) has been perturbed. Thus, t^^, represents the temporary base point about 
which the next exploratory move for the second independent variable, CL, is made. 
To retain consistency of notation, let 



t 



10 




(6-30) 



To find tj 2 ’ T e. , the temporary head obtained from perturbing the second 
independent variable, the response for first determined. Since "5^^ 

represents an increment to a of 0. 1, the response for 

Si (6-31) 

must be found. The performance index for these parameter values yields 

(t , o , t ) = (3. 2, 0, 8) (6-32) 

r s s 

Therefore, 

f(l. 2, 2. 1,0. 5) = [ 3 . 2-2]^ + [0-15]^ + L8-4]^ = 242.44 (6-33) 



The response for 

Sl"'^2 " 1.9, 0.5) 

is now found. The performance index becomes 

(t , o , t ) = ( 3, 0, 6. 8) 
r s s 



Therefore, 

f(1.2, 1.9,0. 5) = L3-2 ]^ +Lo-15]^ + L6. 8 - 4 ]^ = 233.8 



Since 

f(Tll- 62 ) 



(6-34) 



(6-35) 



(6-36) 



(6-37) 



the next temporary head is taken to be 
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*12 * 11'^2 



To find, the third and last (in this example) independent variable, 6, is per- 

turbed by 6^. First, the response for 



+'5’^) = (1. 2, 1. 9, 0. 55) 

is found. 

The performance index yields 

(tr’ °s’*s^ " (3. 8, 0,5.8) 

Therefore, 

f(l. 2, 1. 9, 0. 55) = [ 3 . 8-2]^ + [ 0 - 15 ]^ + [ 5 . 8 - 4 ]^ = 231. 48 



The response for 



hz" *^3 " 



is now found. The performance index yields 

(t^, o , t ) = (3. 9, 0, 5. 5) 
r s s 

Therefore, 

f(l. 2 , 1 . 9 , 0 . 45) = [ 3 . 9 - 2 ]^ + Co- 15 ]^ + [ 5 . 5 - 4 ]^ = 230. 86 



Since, 



f(ti2--?3)^f(ti2-l63)< f(t^^) 



(6-38) 



(6-39) 

(6-40) 

(6-41) 

(6-42) 

(6-43) 

(6-44) 



the next temporary head is taken to be 



*13 *12‘^3‘ 



(6-45) 



All the independent variables have now been perturbed; the last temporary head, 
tj^^, is taken to be the new, or second, base point, b^ . 

With the establishment of the second base point, b^. the pattern strategy 
is brought to bear. The pattern move uses the information generated in the explora- 
tory moves which, initially, establish the pattern and later modify the pattern, as 



L 




i 

1 



169 



the example is to demonstrate. The pattern moves are formulated upon the phi- 
losophy of ’’nothing succeeds like success". Based upon the success in the past, 
of motion to from b^, the assumption is made that continued motion in the 

same direction remains the most desirable course of action. Thus, the new 
temporary head for the second set of exploratory moves is obtained by doubling 
the magnitude of (b^-b^) and proceeding in the same direction. Hence, 



Therefore, from (6-22) and (6-45) 


(6-46) 


r^o = 2b2-bi = 2(1. 2, 1. 9. 0. 45) -(1,2, 0. 5) 
= (2. 4, 3. 8, 0. 9) -(1,2, 0. 5) 


(6-47) 


= (1.4, 1. 8, 0. 4) 


(6-48) 


The performance index becomes 

(t , o , t ) = (4, 0, 4. 2) 
r s s 


(6-49) 


Therefore, 

f(1.4, 1.8, 0.4) = [4-2]^+L0-15]^ + [4. 2-4]^ = 229. 04 


(6-50) 


The response for 

= (1. 4+ 0. 2, 1. 8, 0. 4) = (1. 6, 1. 8, 0. 4) 
is now found. The performance index yields 


(6-51) 


(t , o , t ) = (3. 0, 4. 8) 
r s s 


(6-52) 


Therefore, 

f(1.6, 1.8, 0.4) = [3-2]^ + L0-15]^ + [4.8-4]^ = 226.4 


(6-53) 


It is found that 


(6-54) 



If one reasons as previously with regard to the gain, K^. the calculation for 

t - can be skipped and t,. can be taken to be 
20 1 ^ A 
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(6-55) 



*21 *20'* 


(6-55) 


The search for t^^ now begins. To this end, the response for ^2] ^^2 
where 


is found, 


The performance index yields 


(6-56) 


(t^^Os, ts) ={3.1,0, 4. 5) 


(6-57) 


Therefore, 

f(l. 6. 1. 9, 0. 4) = L3. 1-2]^+ [o-15J^+ [4. 5-4]^ = 226. 46 


(6-58) 


The response for 

4r"^2 =(1.6, 1. 7, 0. 4) 

yields 


(6-59) 


(t , o , t ) = (3. 0, 0, 4. 0) 
r s s 


(6-60) 


Therefore, 

f(l. 6. 1, 7, 0. 4) = [3-2j^+ L o-15J^+ [ 4-4]^ = 226 


(6-61) 


Again, since 

the next temporary head is taken to be 


(6-62) 


T = r -T 

2Z 21 2 


(6-63) 



To find, the last independent variable, S, is perturbed by 6^. First, 

the response is found for 



*22+ ^3 = n. 6, 1. 7, 0. 45) 
The performance index yields 


(6-64) 


(t ,0 , t ) = (3. 8. 0. 6. 0) 
r s s 


(6-65) 



Therefore, 






f 
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f(1.6, 1.7,0.45) = L3.8-2j^+L0-15]^4[6-4]^ = 232. 24 
The response is found for 

^22" S " 

The performance index yields 

(t^, o , t ) = (2. 8, 0, 4. 3) 
r s s 

Therefore, 



( 6 - 66 ) 



(6-67) 



( 6 - 68 ) 



f(l. 6. 1. 7, 0. 35) - [2. 8-2]^+ [0-15j^+[4. 3-4]^ = 225. 



73 



Since 



^^*^22'^3^ ^^‘22^ ^^‘22’'’ ^3^ 



(6-69) 

(6-70) 



the next temporary head is taken to be 



T = T -T - 

23 22 3 



(6-71) 



Since all the independent variables have again been perturbed the last temporary 
head, t^^. is taken as the new, or third, base point, b^. The process is continued 
from (6-45) and (6-71) is found 



t3o = 2b3-b2 (6-72) 

= 2(1. 6. 1. 7, 0. 35) - (1. 2, 1. 9, 0. 45) 

= (3. 2. 3. 4, 0. 7) - (1. 2, 1. 9, 0. 45) 

- (2. 0. 1. 5. 0. 25) (6-73) 

It is noted that repeated success in the same direction, for all the independent 
variables, has caused the pattern to grow. Thus, although only two iterations 
have been conducted, the gain value, K^, is at 2. 0 after starting at 1, 0 and being 
constrained to basic step sizes 5^ =0. 2. 

The performance index becomes for t 

(t , o , t ) = (2. 0, 0, 3. 6) 
r s s 



(6-74) 



I 



I. 





172 



The refore , 

f( 2 , 1. 5, 0. 25) = L2-2]^+Lo-15]^+ U. 6-4j^ = 225. 16 ( 6 - 75 ) 



It is noted that the set of values 

(K^,a,0) = (2. 0, 1. 5, 0. 25) (6-76) 

is the first triple for which the specifications are met by all three elements of 

the vector performance index (t , o , t ). The results of the search to this point 

r s s 

are encouraging. The artificial function f which has served the purpose of 
finding a surface of acceptable performance no longer provides a suitable com- 
parison test; the magnitudes of the individual differences 

It - T 1 , lo - O 1 , It - T 1 (6-77) 

r r s s s s 

are monitored in an attempt to distinguish "best” performance. The response for 



t^O +^1 25) (6-78) 

yields a performance index 

(t , o . t ) = ( 1 . 7, 0 , 3. 8 ) (6-79) 

r s s 

Therefore, 

f ( 2 . 2 , 1 . 5 , 0 . 25) = Cl. 7 - 2 ]^+ C0-15]^+ Ls. 8-4]^ = 225. 13 (6-80) 

The response is still without overshoot, so the exploratory move 
again omitted and 

The choice has been made, not because 

£(^30+ 5 1 ) < £(” 30 ) (6-82) 

although that is true, but because the response £or ^ 30 ^ is £urther inside 
speci£ications than that £or Tjq. This is discussed in detail below. 

The search £or is begun. First, the response £or 

T-. + T^ - ( 2 . 2 , 1 . 6 , 0. 25) 

^ i ^ 



(6-83) 
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is found. The performance index gives 



(tj.. O^. t^) = (2. 0, 4) 

Therefore, 

f(2. 2. 1. 6, 0. 25) = [ 2 - 2 ]^+ L 0-15]^+ [ 4 - 4 ]^ = 225 



F or 



Sr^2 = u. 2, 1.4, 0. 25) 



the performance index yields 

(tr» o^, t^) = (2, 0, 3. 5) 



(6-84) 



(6-85) 

( 6 - 86 ) 



(6-87) 



At this point, the investigator faces a quandary. With the boundary of 
acceptable performance pierced, and the artificial function f effectively dis- 
carded, how is the decision regarding the "best" temporary base point to be made? 
The three possibilities, t^^, 3,nd yield performance indices 

( 1. 7, 0, 3. 8), (2, 0, 4) and (2, 0, 3. 5), respectively. All are within specifications; 
which is best? While 6^ yields a set just meeting requirements, and on that 

basis could be discarded, how can a decision be made between t^^ and 
Obviously, what the one gains on rise time the other gains on settling time. Two 
possible techniques suggest themselves: a)a coin is tossed, i. e. , one invokes the 

Principle of Insufficient Reason, b)the specifications are reduced; a new artificial 
function, f^, is formed, and another comparison measure of which set (K^,a, 0) 
is ’’best" is thereby provided. 

To this end, a new set of requirements 

T = 1. 5, 0=15, T =3. 0 (6-88) 

r s s 

are specified; it should be noted that rise time and settling time have been 
retained at the same ratio as previously. Under these circumstances, 

= f,(2. 2, 1. 5. 0. 25) = [l. 7-1. 5 ]^+Lo- 15]^+ b. 8-3]^ = 225.68 

1 31 1 (6-89) 

^l^bl'^V = f, (2. 2, 1.6,0.25) = [2-1. 5]^+[0-15]^+[4-3]^ = 226.25 

1 (6-90) 

f.(t',,-'6^) = f,(2. 2, 1. 4, 0. 25) = [2-1. 5]^+ [0-15]^+ [3. 5-3]^ = 225.50 

r 31 2 1 (6_9D 
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Therefore, since 






(6-92) 

the next temporary head is taken to be 

b2"Sr^2 

To find, t^^, the variable, d, is augmented by 6 ^. First, the response 

(6-94) 

is found. The performance index yields 



for 



t^ 2 + 6 ^ = (2.2, 1.4, 0.30) 



(t ,o ,t ) = (2,0, 3.5) 
r s s 



Therefore, 



fj(2, 2. 1. 4. 0. 30) = [2-1. 5]^+L0-15]^+ [3. 5-3]^ = 225. 50 



The response for 



(2. 2, 1.4, 0. 2) 



is found. The performance index yields 



(t . o , t ) = (1. 9, 0, 3. 4) 
r s s 



(6-95) 

(6-96) 

(6-97) 

(6-98) 



Therefore, 



Since 



f (2. 2, 1. 4 , 0. 2) - ll. 9-1. 5]^+ L0-15]^+ [3. 4-3^ = 225. 32 
1 /A. 



f, (t33- 63 ). (. 3 ^+ 63 ) 

the next temporary head is taken, to be 



*^33 b 2 ‘ S ^^4 



(6-99) 

( 6 - 100 ) 

( 6 - 101 ) 



Since all the independent variables have again been perturbed, the last 
temporary head, T^y is taken as the new, or fourth, base point, b^. From (6-71) 
and ( 6 - 101 ) is calculated 
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^40 ‘^3 ( 6 - 102 ) 
= 2(2. 2, 1. 4, 0. 2) -(1.6, 1. 7. 0. 35) 

^ (4. 4. 2. 8, 0. 4) - (1. 6, 1. 7, 0. 35) 

= (2.8,1.1,0.05) (6-103) 

The response for, t^g> is determined and the characteristics yield 

(tr- Os> ts) = (1. 4, 0, 2. 75) , (6-104) 

1 herefore , 

fj(2. 8. 1. 1, 0. 05) = [l. 4-1. 5]^+ [0-15]^+ [2. 75-3]^ (6-105) 



The new function, f , , must now be discarded since the set of values for K ,a,3» 
1 c 

(6-103), provides a response which again is within the specifications on 

Thus, a second new comparison function is formed for the changed requirements 



T =1.0, 0=10, T=Z 
r s s 



Therefore, 

£^(2. 8, 1. 1, 0. 05) = Cl. 4-1]^+[0-10J^+L2. 75-2]^ = 100. 72 , 
The response for 

T^O+^l " (^'0. 1. 1,0.05) 



(6-106) 

(6-107) 

(6-108) 



is now found. The performance index yields 



(t , o , t ) = (1. 4, 0, 2. 3) (6-109) 

r s s 

Therefore, 

f (3. 0, 1. 1. 0. 05) = L 1. 4-1 J^+ L0-10]^+ [2. 3-2]^ = 100. 25 
^ ( 6 - 110 ) 

Since, the response is again without overshoot, the exploratory move (t 3 o"^i^ 

again omitted. Also, 



so 



'41 ^40"^ ^1 



( 6 - 112 ) 
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The search for, begun. First, the response is found for 



141+62 = (3. 1. 2. 0. 05) 


(6-113) 


The performance index is 

(tj., 0 ^, t^) = ( 2 , 0 , 2 . 9 ) 


(6-114) 


Therefore, 

f^O, 1 . 2 , 0 . 05) = L2-1]^+ Lo-10J^+C2. 9 - 2 ]^ = 101 . 81 

For 


(6-115) 


= (3. 0, 1. 0. 0. 05) 


(6-116) 


the performance index gives 

(t ,0 , t ) = (1. 7, 0, Z. 4) 
r s s 


(6-117) 


Therefore, 

^2^^4l'"^2^ = f^(3, 1, 0. 05) = [ 1 . 7 - 1 ]^+ [ 0 - 10 ]^+ [ 2 . 4 - 2 ]^ = 100. 65 


(6-118) 


Thus, since 

the next temporary head is taken to be 


(6-119) 


‘42 ^ ^41 


(6-120) 



To find, t^^ the variable, 3, is increased by 6^. First, the response is found 
for 



The performance index yields 


(6-121) 


(t ,0 , t ) = (1. 8, 0, 3. Z) 
r s s 

Therefore, 


(6-122) 



f^(3, 1. 1, 0. 1) = L 1. 8 -l]^+ L0-10]^+ [3. l-zf = 102. 08 ( 6 - 123 ) 



The response for 


(6-124) 
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cannot be used due to practical limitations on the analog computer, i. e. , due 
primarily to the layout used to solve the problem. An even more important 
justification for not reducing 6 to zero, however, is that it is not possible to 
realize such a compensation network, within the given class to which the com- 
pensator design is constrained. Furthermore, since the artificial function has 
been reduced twice already, the design is well within original specifications, for 
C = 1.0. Thus, since 



1 . e. , ~ ^2^^ 40"'^ ^2^*^42^ 



100. 25 < 102. 08 



(6-125) 

( 6 - 126 ) 



the next temporary head is taken to be 



t 



43 




(6-127) 



Since all the independent variables have again been perturbed, the last temporary 

head, T.,, is taken as the new, or fifth, basepoint, h^. From (6-97) and (6-108) 
43 ^ 

is calculated 



= 2(3, 1. 1. 0. 05) - (2. 2, 1. 4. 0. 2) 

- (6. 2. 2, 0. 1) - (2. 2. 1.4, 0. 2) 

= (3. 8, 0. 8, -0. 1) ^3. 8, 0. 8, 0. 05) 



(6-128) 



(6-129) 



That is, under the practical constraint that 8 not be less than 0. 05, 3 is set 
equal to 0.05 whenever the process indicates 8 < 0. 05. 

The response for 

(3. 8, 0. 8, 0. 05) (6-130) 

is determined. The performance index yields 

(t , o , t ) = (1. 4. 0, 2. 5) (6-131) 

r s s 






8.0.8,0.05) = [ 1. 4-l]^+ [0-10]^+ [2. 5-2]^ = 100.41 (6-132) 



Therefore , 
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The response is now found for 



t^ 0 + = (4. 0, 0. 8 , 0. 05) 



(6-133) 



The performance index yields 



(t , o , t ) = ( 1 . 2 , 0 , 2 . 0 ) 
r s s 



(6-134) 



Therefore, 



£ {4,0.8,0.05) = [ 1. 2-l]^+ L0-10]^+ L2-2]^ = 100. 04 (6-135) 



Again, the exploratory move, omitted, and since 






(6-136) 



sr- 

The search for, begun. First, the response is found for 



£ 51 + 5^ = (4. 0. 9, 0. 05) 



(6-137) 



The performance index gives 



(t ,0 , t ) = ( 1 . 1 , 0 , 2 . 0 ) 
r s s 



Therefore, 



(6-138) 



f.,(4, 0 . 9 , 0 . 05) = L 1 . 1-1 J^+ L 0 - 10 ]^+ [ 2 - 2 ]^ = 100 . 01 (6-139) 



For, 



£ 51 - 5^ = (4, 0. 7, 0. 05) 



(6-140) 



the performance index yields 



(t ,0 , t ) = (1. 3, 0, 2. 5) 
r s s 



(6-141) 



Therefore, 



f ( 4 , 0. 7, 0 . 05) =[ 1 . 3 - 1 ]+ [ 0 - 10 ] + [ 2 . 5-2T = 100 . 34 
2 (6-142) 



Thus, since 



r 
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f 2 (t 5 i+ 62 ) < < f 2 (tsrV (6-143) 

the next temporary head is taken to be 

I52 “ I5 ^2 (6-144) 

To find, t^^, the variable, 0, is increased by 6^. First, the response is found for 

"^52 + ^3 = (4. 0. 9, 0. 1) (6-145) 

The performance index yields 

(tr' ®s’ ^s^ " (6-146) 

Therefore, 

1 ^( 4 , 0. 9, 0. 1) = [1. 2-1]^+L0-10]^+ L2. 2-2]^ = 100. 08 (6-147) 

Again, the response for t>e found so the "^ 2"*^ ^3 response 

is matched against the t^^ response and it is found that 

f(F^2^ l’('i52'^^3^ . (6-148) 

Therefore, the next temporary head is taken to be 



t 



53 




(6-149) 



Since all the independent variables have again been perturbed, the last temporary 
head, 1^^* taken as the new, or sixth, base point, b^. From (6-127) and (6-149) 
is found 



t 



60 




(6-150) 



= 2(4, 0. 9, 0. 05) - (3. F 1, 0. 05) 

= (8, 1. 8, 0. 1) - (3, 1. 1,0. 05) 

= (5, 0. 7, 0. 05) (6-151) 



The response is determined for 

= (5, 0. 7, 0.05) 



(6-152) 



The performance index yields 



I 

i 

I 
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(tr- °s’ = n. 1. 0, 2. 1) 


(6-153) 


Therefore, 

1^(5, 0. 7, 0. 05) = [ 1. l-l]^+ [0-10j^+ [2, 1-2]^ = 100. 02 


(6-154) 


The response for 

= (5. 2, 0. 7. 0. 05) 

is now found. 

The performance index yields 


6-155) 


(t , o , t ) = (1. 1,0,2. 2) 
r s s 

Therefore, 


(6- 1 56) 



f^(5. 2, 0. 7, 0. 05) = Ll. 1 -1 ]^+ [0-10]^+ [2. 2-2]^ = 100. 05 (6-157) 



The response for 

r^O'^1 = 05) 

if found. The performance index yields 


(6-158) 


(t , o , t ) = (1. 1, 0, 2. 2) 
r s s 


(6-159) 


Therefore, 

f^(4. 8, 0. 7, 0. 05) = [ 1. l-l]^+ [0-10]^+ [2. 2-2^= 100. 05 

Since 


(6-160) 


^2^*60^ ' ^2^‘60''’ '1^ " ^2^'^60‘^1^ 
the next temporary head (basepoint) is chosen to be 


(6-161) 


^61 "" ^60 

The search is started for response is found for 


(6-162) 


r, , + = (5, 0. 8, 0. 05) 

b 1 c 

The performance index measures 


(6-163) 



(t , o , t ) = (1. 05, 0. 2) 

r s s 



(6-164) 
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Therefore, 



f2(5, 0 . 8, 0 . 05) = Ll. 05-1 ]^+ [ 0-1 0 ]^+ [ 2 - 2 ]^ = 100 . 0025 


(6-165) 


For, 

^61" ^2 ^ 

the performance index yields 


(6-166) 


(tr» o^, t^) = (1. 1, 0, 2. 0) 


(6-167) 


The refore, 

0. 6. 0. 05) = Ll. l-l]^+ L0-10J^+ [ 2 - 2 ]^ = 100. 01 


(6-168) 


Since, 

^2^*61''' ^2^ ‘ ^2^^61‘ ^2^ ^2^^61^ 
the temporary head is taken to be 


(6-169) 


^62 " ^61^ ^2 


(6-170) 



To find, the third independent variable, 3 , is increased by, 6^, and the 

response is found for 



^62“^ ^3 

The performance index yields 


( 6 - 171 ) 


(t , 0 , t ) = (1, 0. 2) 
r s s 


( 6 - 172 ) 


Therefore, 

12(5, 0. 8, 0. 1) = [ 1-1]^+ L 0-10]^+ L 2-2]^ = 100 


( 6 - 173 ) 



Again, the response for 6^ cannot be found so the ^3 response is 

matched against the t^^ response and it is found that 



^^*62’'’ ^^*^ 62 ^ 

Therefore, the next temporary head is taken to be 


( 6 - 174 ) 







i 




( 
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63 



^^ 62 -^^ 3 - 



(6-175) 



Since a.11 the independent variables have been perturbed, the last temporary head, 
*^63’ taken as the new, or seventh, basepoint, b,^. Thus, from (6-149) and 
(6-175) 



*70 ^** 7‘’’6 



(6-176) 



= 2(5, 0. 8, 0. 1) - (4, 0. 9, 0. 05) 

= (10, 1. 6, 0. 2) - (4, 0. 9, 0. 05) 

= (6,0.7,0.15) (6-177) 



The response for 

= (6, 0. 7, 0. 15) (6-178) 

can be deternnined, etc. The process can be continued, thus providing much 
additional information regarding the correlation between (K , a, 0) space and the 

C 

space of (t^, o^, t^). A solution providing satisfactory response across the range 
of varying C is (K^, a, 3) = (5, 1. 0, 0. 05). The values of the performance index 
for various C are displayed in Table 6-1. It is clear from Table 6-1 that, 

5, a— 1. 0. and 3 = 0. 05 form a satisfactory compensator- which meets the 
original specifications. 

TABLE 6-1 



K 

c 


a = 1. 0 


3 = 0. 05 


c 


t 

r 


O 

S 


t 

s 




1. 0 


1 


0 


2 




0. 9 


0. 95 


0 


1. 95 




0. 8 


0. 8 


0 


1. 6 




0. 7 


0. 5 


0 


1. 3 




0. 6 


0. 3 


0 


1. 0 




0. 5 


0. 25 


0 


0. 9 




0. 4 


0. 2 


0 


(T7? 




0. 3 


0. 17 


0 


OTTC 




0. 2 


0. 15 


2 


0. 17 




0. 1 


0. 1 


4 


0, 15 




0, 05 


0. 1 


7 


0. 5 




0. 01 


0, 1 


9 


0. 9 



k 
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- 2, = 15, = 4 (6-14) 

for all C in the range, 0. 01 < C < 1.0. Indeed, this compensator meets the more 
stringent requirements of 

T = 1, O = 10, T =2 (6-106) 

r s s 

across the entire range specified for varying C. 

A review of the example indicates how quickly the search converges upon 
a solution; a digital computation program can, of course, provide the answer 
obtained with considerably more speed and accuracy. The addition of the delay 
time (pure lag), which is the major change to the example, introduced by the 
human operator, in now way modifies the search procedure used. The delay 
interval has its effect upon the system response characteristics and may reduce 
the prospect of success in finding a "best-overall” compensator of a given class 
but the direct search method with pattern strategy retains its general applicability. 

It should be noted again that, in the example just offered, the human- 
operator transfer function was set equal to unity and the variable-parameter- 
mapping plots used to estimate the approximate values that a and 3 must have 
if compensation is to achieve the purpose prescribed. Thus, the three-variable 
mappings were allowed to degenerate to a two-variable representation and the 
values of a and 3 in this instance, had a direct effect upon the location of the 
critical points governing stability. Normally, when the Q mappings are used to 
represent the human- ope rator adaptive range the influence of a and 3 values 
upon system stability is not so easily discernible. In such cases, the selection of 
values for a and 3 influences the shape of the curve, which represents 

the system dynamics; the manner in which this curve passes throught the Q 
regions is, of course, what ultimately determines overall (with a man in the loop) 
system stability and performance. Thus, in an analysis with a human operator 
in the system, the effects of the a and 3 values are only indirectly observable 

before they are selected. Upon selection, however, if the compensated -L^ is 
plotted the effects are clearly visible. This situation is no different than that 
routinely faced by the control engineer when dealing with less animate systems. 
Experience is quickly gained concerning the results of introducing particular com- 
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pensator configurations (picking values for a and 8); to this end. with a definite 
goal for the shape of the modified curve the Bode diagram (more quickly- 

plotted than a Nyquist curve) can be used to advantage. 

6. 7 Summary 

A discussion of scalar performance indices and some of their limitations 
is followed by a laconic survey of methods for solving the "best” performance 
problem, especially for vary ing-parameter plants. The notion of a "vector" 
performance index is introduced and a technique for establishing one response as 
better ’ than another is presented. This technique rests intimately upon the 
formulation of an "artificial" function which when reduced to zero identically 
satisfies the upper bounds on system requirements. These requirements are set 
forth in the "vector" index. The design for a "be st- ove rail" compensator to 
meet specifications upon a varying-parameter plant is obtained by using a "direct 
search" method {especially suited to the electronic computer). This search 
method incorporates a "pattern strategy" to aid in finding the next set of trial 
values for the compensator; the strategy affords rapid convergence to a solution 
as is demonstrated by an example. 

Finally, the example itself, displays the versatility of the variable-para- 
meter mappings presented in Chapters Two-Five. By setting the human-operator 
transfer function equal to unity, it is shown that the mappings have use for 
"inanimate" system design. 

Empirical justification of the conclusions from this theory (with a human 
operator in the loop) must wait upon additional research and the publication of 
experimental results issuing from deterministic inputs driving controlled elements 
other than pure gains. It is recognized that the bulk of the evidence justifying the 
form of human-operator transfer function used in this study has been accumulated 
for systems driven by stochastic input signals. Still, the initial range selection 
for the parameters K^, , T^, T^, and T^ included the effects of calculations pub- 

lished for deterministic inputs and modifications, if needed, are easily incor- 
porated. With a technique provided for comparing one multi-dimensional index 
against another, the methods of this chapter offer the prospect for more realistic 
performance measures than the standard "figure of merit" of scalar variety affords. 
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APPENDIX A 



A. 1 Human Operator Zero Locus 

It is desired to establish (4-13) as the equation of a circle. The 
formation (A-1) is to be considered 

a+jb 

w - u + jv “ — 

14x+jy 

Of interest are the c onditions x 0, y = x’T , thus 

i— ( 

a + jb (l+x)-jy 

w = ^ 

(1 +x)+jy (l+x)-jy 

a( 1 +x)+by+j [b( 1 +x )-ay ] 

^ 2 2 

(1+x) +y 

The real and imaginary parts of each side of(A-2) are set equal 



Similarly, 



a ( 1 +x) +by 

" 2 2 

(1+x) +y 

b( 1 + x) -ay 

~ 2 2 

(1+x) +y 



1 +x+jy 



a + jb u - j V 
u+jv u-jv 



au+b v + j (bu-a v) 
1 .x+jy - 

u + V 



trans - 



(A-l) 



(A-2) 



(A-3) 



(A-4) 
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therefore 



+ x 



auibv 
u +v 



(A-5) 



, y = 



bu - av 

2 2 
u +v 



(A-6) 



For X - 0, (A-5) gives the locus in the space of the mapping and (A-6) permits 
recovery of T for any given t 



2 2 , 
u +v -au-bv = 0 



(A-7) 



If the square in (A-7) is completed, one obtains the familiar equation of a circle 




(A-8) 
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APPENDIX B 



B. 1 The Parabolic Effect 

A typical plot of a face from (T , T , T ) space, T constant, appears in 

E i N E 

Fig. B-1. Without loss of generality, a set of xy axes is shown so oriented as to 
simplify the following exposition. Of interest is the complete mapping corres- 
ponding to face ODCEO in Fig. B-2. This mapping can be found by allowing a 
typical element of this face such as a0 in Fig. B-2 to scan from line OD to line 
EC. The effect, as has been seen in Chapter Four, is a plot similar to that 
shown in Fig. B-l, except for translation of axes. When a is halfway down line 
DC, g is halfway down OE, Fig. B-2, and similarly in the mapping space. 

Fig. B-l, i.e. , linearity holds through the mapping as is proved later in this 
Appendix. Thus, is also halfway down the line D^C^ of Fig. B-l as is 
halfway along O^E^. Each of these points is given arbitrary coordinate 
designation in Fig. B-l and it is noted 




for this choice of axes. The mapping of any element ag is obtained by joining the 
points in Fig. B-l. As shown in Fig. B-3 the family of these lines extends 

outside the end point mapping of Fig, B-l. It is shown that the envelope of this 
family of lines is a parabolic arc as suggested by the many dotted lines of 

this family drawn in Fig. B-3. For the geometry of Fig. B-l 




y 

/ 

a 

x 



>^(y3-Y5) = 



X(x^-x^) 



= Xx 



2 



a 

y 



X { V - y ) + y 



where the subscripts indicate the x and y components of the moving points a and 
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y 




Fi^. B-1 : Typit al Map from T^, Space 
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Fig. B-2 : Scanning a Face in (T , T , T ) Space 
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Fig. B-i : The Parabolic Effect 
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8^ in the space of the mapping. Froni these designations may be written the slope 
of the moving line: 



^ _ 



B'-a' -^(y,-y ,)-y , 

V V _ 3 2 4 4 



dx 



6 -a 

X X' 



0-\x^ 



(B-1) 



Either by integration and evaluation of the constant of integration (when x = 0, 

Y = \y^) or through use of the slope intercept form y = mx+b the equation for the 
family of lines may be written 



y ^ 



y2-^3-y4 






From the theory of envelopes it is known that if 



( 1 ) 



(B-2) 



F(x, y. >.) - 0 



(B-3) 



represents a family of curves, then differentiating with respect to holding x 
and y constant and eciuating to zero, yields 



/ it 

^(x.y.X) 



(B-4) 



The envelope is obtained by eliminating ). bi-tween (B-3) and (B-4). 

y. 



F(x,y,)) y 



'^2-''3-y4 M , - 



and 



F^(x,y,X) 






X - V 0 



>'4^ ^ x'2>. y. 




(B-5) 



Ford, L.R., "Differential Equations", McGrawHill, New York, 1955. 
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K (B-4) is substituted into (B-2) one obtains, 



''z-H-U /4 f2’'3 

y - X + — / 

^2 ^4^ 



X + 



^4- 



2^3 



(B-6) 



Equation (B-6) upon simplification becomes 



^ _>' 2 -’' 3-''4 , jU'<i /'' 4''3 

y - :: X + / -7 — yx + / — ^ — yx 



(B-7) 



y = 






X + 2 



'^4^3 



(B-8) 



A change in notation yields 



a = 



y2-y3-^4 



and b = 2 



'y4Y3 



y - ax 



b^/jT 



(B-9) 



Equation (B-9) squared becomes 



2 2 2 

y - 2axy + a x - bx 



(B-10) 



Equation (B-10) is a quadratic form whose parabolic nature is not 

obvious. Standard formulas for transformation of coordinates to principal axes 

( 1 ) 

form are available. Briefly, a quadratic form 



2 2 

Ax +Bxy + Cy + Dx + Ey + F = 0 



(B-11) 



may have its axes rotated by an angle a so as to eliminate the cross-product 



( 1 ) 



Thomas, G. B. , Jr., "Calculus and Analytic Geometry", Addison- 
Wesley, 1956. 
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term. A general rotation applied to the axes of Eq. (B-11) results in a new 
quadratic form 



a' x' +B^x^y^+C^y^^ + D^x^+E^y^-I-F^ =0 



(B-12) 



with coefficient relations given by 

/ 2 2 
A = A cos a + B cos a sin a + C sin a 

/ 2 2 

B = B(cos a -sin a) + 2(C-A) sin a cos a 

/ 2 2 
C =Asin a-Bsinacos a + Ccos a 

- Dcos a + Esina 

E^ - -D sin a + E cos a 

F' = F 



(B-13) 



If (B-11) has coefficient B ^ 0 then an angle a can always be chosen such that 
B^ =0; this angle may be found from 



cot 2a = 



A>C 

B 



(B-14) 



With the above information and some trigonometric manipulation one 
finds that since 



then 



So, 



cot 2a 



1 -a 
2a 



cos 2a - 



1 -a 

a^l 



Sin a 



ttt : 



(B-15) 
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and 



cos a 



/a^+l 



(B-16) 



Use of (B-15) and (B-16) in Eqs. (B-13) produces 



,2 ,2 ab 

(a +l)y + y 



yr^+T y a ^+1 



x' = 0 



(B-17) 



which has the form of a parabola. However, let 



h = 






and k = 



ab 






then 



y ^ ^-hx ^ + ky ^ =0 



and completing the square 



(y' + |)^-^- hx' =0 



Let 



and 



then 



k ^ 



t . ^ 



y = y + 2 



// t , Is 

^ = X + ~ 

4h 



// ^ , / 
r = hx 



(B-18) 



Equation (B-18) is the simple form of the equation of a parabola. 

It should be noted that the area added to the end point mapping to form 
the complete hull is very small. In the light of this smallness, its effect upon 



195 



the gross hull over the range of T , and its location relative to the plots of 
most physical systems, there is sufficient justification to call the mapping 
"almost complete" in the spirit of Section 4.3. 

B. 2 An Exploration of T^., Space 

The prospect of defining exactly the mapping of the rectangular 
parallelipiped of Fig. B-2 onto the Q plane must be considered. Previous results 
were for the mapping of the boundaries of this volume only. As above for face 
ODCEO, the envelope of the region in Q space corresponding to the mapping of 
every single point within the volume shown in Fig. B-2 must be obtained if both 
necessary and sufficient conditions for stability are to be provided for the dual 
Nyquist plot. 

The situation (u = 1 given in Fig. B-4 (same as Fig. 4-18) is first 
considered. It should be noted that face 1243 of Fig. B-2 maps into the triangle 
plus short line segment 1^2^4^3^ and face 5687 maps into the triangle and short 
line segment 5^6^8^7^. Indeed, as face 1243 is moved along lines 1-5, 2-6, 3-7 
and 4-8 it eventually coincides with face 5687. Thus, every point of the volume 
in (T^, T^, T^) space is scanned as face 1243 is moved along for varying T^, 
i. e. , as T goes from zero to 5. 3. The effect in the complex plane plot is to 

Lj 

have the triangular shaped mapping move from the position 1 2 4 3 to that of 
5 '6 '8 ' 7 '. This motion occurs under the constraint that point Z‘ moves along arc 
2 '6', i. e. , the motion of point 2 along the line 2-6 in (T^, T^, T^) space corres- 
ponds to motion of 2 ^ along the circular arc 2 6 in the space of the mapping. 
Similar situations pertain to points 1^,3 and 4 in the complex plane as they 
move toward coincidence with 5^,7 and 8 , respectively. 

It is desired to show that lines 2^-4' and 4 '-3', as they move toward 
lines 6^-8^ and 7^-8^ never contain points which lie outside the bounds of the 
hull plotted in Fig. B-4. Later, the results of the parabolic effects displayed 



ssiii:'* 




TO 2‘(C*) 



3 





1 6 



Fig. B-4 : Effect of Varying T 
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in the first part of this Appendix are considered. 

For the situation of Fig. B-4, i.e., uu=l, the line E'C '{4'-2') has end 
points 1 + j25 and - 1 6, 5 + j25. 7, therefore its equation is 



y- 25 X- I 

25. 7-25 " -16.5-1 

9c; (•’ 7 ) , M 

y-25 

or 

y= ^^{x-l)+25 (B-19) 



i.e., the line has slope 



m 



1 



-. 7 
17. 5 



The end points of the radius arm of arc F E are 1 + j25 and the center 
ofarcF^E^, i.e., ^+jT2.5. The equation of the radius arm is therefore 



r-25 X- 1 



12.5-25 J_ ^ 
2 



- 1 2 5 

y-25 = (x-1) 

'2 



y = 25(x-l) + 25 



(B-20) 



i. e. , the line has slope 



Since 



= 25 



m.m. =7^(25) = = -1 



T 2 17. 5 



17. 5 






ynVU .. 

; ' ■ .. ■! : 
■ '■■■ ' 
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then line E'C' is perpendicular to the radius arm Opj,-E' or 

is tangent to arc F'e'( 4'8'). The above is true for uu = 1 . It must be shown 

true, in general, for any a\ 

For E'c': The end points are l+j25it' and 1 - 1 7. 5 u)^+ j25. 7 uj therefore 

line has equation 



y-25 uu _ x-1 

25. 7UU-250J “ , ^ r ^ 

- 1 7. 5 uu 

. 7 ID __ -.7 

■■ ■ -17.5.)^ ^ 

For radius arm Opg,-E' in Fig. B- 4 the end points are 1 +j25uu and the 
center Op^ =-g 4 ; the equation of Opp-E' is then 

y- 25 UU _ X - 1 
12. 5 uu-25ju " 1 T 



and 



2 



25uu 



m 1 



7 

17. 5 i' 



(35a) 



-17. 5 u 
17. 5 JU 



indicating, again, that perpendicular to line E^C^, and therefore line 

E^C^ is tangent to arc f'e\ 

A similar analysis answers the question, ”Is tangent to arc A'd'? 

For D ^ : The end points are 1 + j. 7 '»j and 1 - 17. 5uu^+25.7juu therefore 

the line has equation 



!t 



X-1 

-17. 5(ju^ 



y- . 7 JU 
25. 7 JU-. 7iV 



1 
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y - . 1 iv 
25uu 



X - 1 

-17. 5u)^ 



and m j 



25uu 

-17. 5u)^ 



(B-21) 



For the radius arm of arc A'd', i. e. , 0 ^q- 1' of Fig. B-4 the 

end points are 1 +j.7u) and the center giving an equation for the 

radius arm of 





Equations (B-21) and (B-22) give 



35uu 

J. 

" 2 



. 7 uu 



(B-22) 



m, m- = j {. 7 ju) = - 1 

^ -17. 5 x' 

and line - is tangent to arc 

Thus, it has been proved that initially line Z'- 4' is tangent to the cir- 
cular arc F'e'( 4'8') at point 4\ If in (Tj^, Tj, space, face 1243 moves 

some distance along the axis toward face 5687, then in the space of the map- 

ping, point Z' experiences some corresponding angular excursion along the arc 
2^6^ toward 6\ Point 4' experiences the same angular excursion along the 
arc 4^8^ toward 8^, etc. 

Actually, it is sufficient to show that for a specified increment along 
equal angular increments occur on all the arc segments, i. e. , it shall be proved 
that T =x T 0 < X < 1 (T =max T ) corresponds, as x varies, to arc seg- 
ment increments having equal angular measure. The proof is facilitated by the 
passing of all circles through the origin. Thus, it is only necessary to show that 
a particular increase along T^ produces equivalent phase change on all the cir- 
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cles. 



For arc B'c" (z' 6^), since 



then, for 0 



and 



^ _ 1-17. 5o) + 25. 7juj 
1 + TJOJ 



Q = (1 - 17. 5uu ) + jZ5. 7uu 
-1 Z5.70J 



cp = tan 



1 - 17. 5uu 






and 



^ (1 - 17. 5uu ) +Z5. 7joi) 
1 + X T^j aJ 



cp = tan ^ - tan ^ (xT^a) 



1 - 17. 5a 

Therefore, the extent of arc Z^6^ is given by 



Acp = tan ^ X T^a 



For arc.F^E^ (4^8^), since 



, 1 + ^ 5 j jL 
1+T, ja 



then, for = 0 



= 1 + j Z 5 a 



and 



cp = tan ^ Z5 a 



For T^=xT^ 



_ 1 +j25x 

\ i. ~ . rp • 



1 + X T^j a 



and 



= tan ^ Z5 a - tan ^ x T, a 



Therefore, the extent of arc 4^8' is given by 



(B-Z3) 




! 
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Acp = tan"^ xT^iu (B-24) 

and similarly, along arc G ^ O ^ (3^7^) and arc (1^5^), 

Since (B-23) equals (B-24) the statement is proved that corresponding 
angular motion occurs on all arc segments. Indeed, the fact follows from the 
properties of the bilinear transformation. ^ 

It can be mentioned at this point, also, that if is set equal to x T^, 
where T^=max and 0 < x < 1 , similar to the above exposition, then the effect 
again is to move correspondingly along the mapped line. For example, Fig. B-5 
shows a plot, 0) = 1 , of the situation where 0< T^<12. 5, rather than 0<Tj<25, 
as originally. The equations of lines D^-C\ O^-E^, A^-B^ and G^-F^ are 
the same as those in Fig. B-4, but the line segments cover only half of their 
previous extent. This property follows directly from the mapping 

(1 + j ijuT ){1 + j 0) T ) 

" = (l4j.Tj- 

For U)=(JUj^, T^= constant, T^^= constant 

n = (1 +j'ji|T|) (Cj+jCj) 

= (C,.C2.,T,)tj(C2+Ci«,T,) 

= K^Tj+K^ (B-25) 

where K.= some complex constant. Equation (B-25) demonstrates that, irrespec- 
tive of the value of Tj., the mapping lies along the same straight line. 

To be considered now is the case where movement has occurred corres- 
ponding to 30*^ of arc along 2 6 , shown as point 2 , and 30 of arc along 4 8 , 



Churchill, R. V. , "Complex Variables and Applications", 2nd Ed. , McGraw- 
Hill, New York, I960. 
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V 




Fig. B-5 : Effect of Reduced Range on 
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shown as 4 , Fig. B-4. The line joining these new points corresponds in 

space to an element of face 2486, Fig. B-2, which is parallel to line 
2-4, or equally well, line 6-8. In the complex plane the angle Z''4'z' equals 15^ 
(an inscribed angle equals one-half its intercepted arc). Also, the angle formed 
by the extended line 2^-4' tangent to arc 4^8' and the line 4'-4'' equals 15^ (the 
angle between a tangent and a chord equals one-half its intercepted arc). These 
two angles are therefore "vertical angles" and points Z'\ 4" can be connected 

by a straight line. Thus, in the space of the mapping, as 2 '-4^ moves toward 
6 -8 , the line connecting the moving points 2^-4^^ always passes through the 
original position of point 4'; thus, the line NEVER passes outside the "end 
point" hull. Of course, this state of affairs pertains only until point z!' reaches 
4^ but after that it is obvious that the line z!* is inside the hull as in the case 
when 2 -4^^ has finally become 6 -8 . Similar conditions prevail in the pro- 
gression of line 3 -4 to coincidence with 7 -8 , i. e. , the line connecting any 
points 3^^ -4^^ always passes through the original point 3^ Thus, except for the 
parabolic envelope effect mentioned at the beginning of this Appendix the end 
point mapping of the volume in Tj, T^^) space would be a complete one. It is 

granted that the area contributed by the parabolic fillet adds little to the complete 
hull that was not already in the end point hull. This can be seen by envisioning 
the triangular shaped region 1 ^ 2 ^4 ^ 3 ^ sweeping down to 5 ^ 6 ^ 8 ' 7 ' , as T^^ ranges 
from zero to 5. 3, only keeping in mind during the sweep the presence of the 
parabolic fillet. If at any time the stability of a physical system is dependent 
upon points in the affected region they then can be taken into account. It is seen 
in Chapter Five that in those systems considered in this study the additional area 
may be disregarded. Due to the presence of this small effect the mapping is 
classified "almost- complete. " 
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APPENDIX C 



C. 1 A Two Parameter System 

The use of the three-parameter Variable-Parameter Mapping of Chapter 
Four for a two varying-parameter problem is considered in this Appendix. An 
interesting stability demonstration occurs in the following system; also, the 
example serves a useful purpose in the performance considerations of Chapter Six. 
The situation arises in the quest for a best-overall compensator which is des- 
cribed in detail in Chapter Six; only that portion of the design task of tutorial 
pertinence is investigated. 

The system and compensator are as shown in Fig. C-1. The system is 
analyzed for values of C = 0. 1, JU^= 1. 57. As in Chapter Four stability obtains if 

1 + F(s) = 0 (C-1) 



has no right half plane zeros, where 



- 


r ^ 1 

i' 


sa+ 1 


n 


s3 + 1 


s( s + ZC ) 
n 



If (C-1). is multiplied by 



(C-2) 



s3 + 1 _ 1 

sa + 1 " ^ 

the fixed and varying portions for the plots are obtained. That is, 

Z 

L' 

s3 -f 1 ^ n _ Q 

sa + 1 s(s + Zt x^) 



(C-3) 



(C-4) 



and y 

n _ + 1 (C-5) 

V sa+ 1 

where V is again taken to be unity as in Section 4. 2, i. e. , a = 0 = 0 
n 

Also, ^ 

L - - 

n S( S + ZC ti;' ) 

n 



(C-6) 




o 
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r(t) 





sa + 1 










s^+ 1 




s(s+2^a;n ) 



Fig. C-1 : Second-Order System with Compensator 
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This results in the standard plotting structure 



V 

n _ 

V " 



-L 

n 



(C-7) 



C. Z The Degenerate Mapping 

The constraints on permissible a and 6 are taken to be 0 < a < 5. 3 

V 

^ ^ ^ apparent that the mapping is obtained from the image 

of face 3487 shown in Fig. 4-16. This mapping is actually shown in Fig. 4-ZO, 
for OJ = 1.0, i.e. , it is represented by the bounds formed by line 3 ^4 \ arc 4 ^-8 ^ 
line 8 7 and arc 7^-3 . The degenerate mapping just outlined is shown in Fig. 
C-2 and again in Fig. C-3, along with the plot of -L given by the negative of 
(C-6). 

The primary purpose of this example is to demonstrate that sufficient 
caution must be exercised in interpreting the plots to preclude succumbing to the 
temptation of making stability decisions based upon limited frequency information. 
One must guard against using the thumb rule of Section 4. Z superficially, by con- 
sidering only one or two frequency points. The remainder of the example clarifies 
this notion. Upon observation of Fig. C-3 one might infer that some set of CX, B , 
such as a = Z, 3 =Z5, would result in a stable system. That is, for 3 =Z5, as a 

V 

varies, the image point of — ^ runs along arc segment 4 -8 , and for a = Z if the 
image point is "to the left’' of (tracing - as frequency increases), such a 

condition seems to imply a stable situation. The point corresponding to a = Z, 
B=Z5 and, of course, i = 1.0 is shown as M(jl)inFig. C-3 where 



That is, 



M 



ja(Z5) + 1 
ji:(Z)+ 1 



V 

M = ^ (a = 2, e =25) 



(C-8) 



and it is observed that M(jl) does fall to the left of - implying stability. 
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Kcu) 




Fig. C~2 : The Degenerate Three-Parameter Mapping 

for Two- Paramete r Analysis 
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Kcu) 




-0. I, =1. 57 

n 
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An actual stability check for these specified circumstances can be obtained 
from the Routh test. From (C-1) and the characteristics given one obtains 



, (2. 46)(s2 + 1) 

(s25 + l)(s+. 314) ■ 



(C-10) 



From (C-10) there follows 

s^+. 354s^+ . 20956s + . 0985 = 0 



3 

s 

2 

s 

1 

s 



1 

. 354 

. 0742-. 0985 
. 354 



0 

s 



(-) 



(+) 



. 20956 
. 0985 



The two sign changes in the array above signify two right-half-plane roots 
in the characteristic equation and, consequently, instability. The apparent in- 
consistency between the graphical indications and the analytical results can be 
cleared away by a more detailed plotting of the M(j uu) locus at more frequencies, 
particularly those less than unity. 

The pertinent points of the locus for M(j uu) are shown in Fig. C-4 along 
with -L (io)). It is noted that the "difference vector" drawn from -L to M 
undergoes a clockwise rotation between M(0. 3) and M(0. 7) thus indicating in- 
stability and consistency with the Routh test. The difference vector actually maJces 
two complete clockwise rotations corresponding to the existence of the two right- 
half-plane roots in the characteristic equation. To observe the complete revolu- 
tions the negative frequency branches of ^ must be considered; also, 

the effects of s = 0 are handled by the standard technique of setting 

s = lim Re*^ ^ 

R-0 
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I(w) 




Fig. C-4 : Two- Parameter Variation for a=2, 8 = 25; also shows -L^ 



to 



where 0 



varies from 



and the origin of the s*plane is circumvented by 

>(i) 



a small semi-circle of radius R' 

Difficulties with "encirclement” interpretation can normally be resolved 
by sketching the entire locus of for positive and negative frequencies 

(similarly for the varying portion); the excursions of the difference vector between 
these two loci may then be more readily observed. 



ti) Del Toro, V. and S. R. Parker, Principles of Control Systems Engineering, 
329-330, McGraw-Hill, N. Y. , I960. 
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Fig. 4-20 : Variable Parameter Mapping 
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Fig. 4-22: Variable Parameter Mapping 



lP 




J' 









o 



c 

t 

c 

► 

f: 

V 



^ r 



O'. 

3 



.• 




c 

•-%■ 

o 



C/5 

n 

cr 

o 

o 



I(w) 




Fig. 4-23 : Variable Parameter Mapping 



thesS557 

Stability and performance of manned cont 




